AND PUBLISHED B 


The 


PHILOSOPHICAL 
MAGAZINE 


FIRST PUBLISHED IN 1798 


. 46 SEVENTH SERIES 


No. 374 


A Journal of ; 
Theoretical Eaperimental 


and Applied Physics 


EDITOR 
PROFESSOR N. F. MOTT, M.A., D.Sc., F.R.S. 


EDITORIAL BOARD 
SIR LAWRENCE BRAGG, 0.B.E., M.C., M.A., D.Se., F-B.8. 
SIR GEORGE THOMSON, M.A., D.Sc.. F.R.S. 
PROFESSOR A. M. TYNDALL, C.B.E.; D.Sc., F.R.S, 


PRICE 15s. Od. 
Annual Subscription £8 0s. 0d. payable in advance 


Y TAYLOR & FRANCIS LTD., RED LION COURT, FLEET ST., LONDON, E.C.4. 


New Interseience Books 


LECTURES ON PARTIAL 
DIFFERENTIAL EQUATIONS 


By I. G. PETROVSKY 
Translated from the Russian by A. SHENITZER 


CONTENTS: Introduction. Classification of Equations. Hyperbolic 
Equations: The Cauchy Problem for Non-analytic Functions. Vibrations 
of Bounded Bodies. Elliptic Equations. Parabolic Equations. 


1954 255 pages, 17 illus. 41/- 


HIGH-ENERGY ACCELERATORS 


By M. STANLEY LIVINGSTON (Massachusetts Institute of Technology) 
(Interscience Tracts on Physics and Astronomy, Volume II) 


CONTENTS: High-Energy Accelerators as Tools for Nuclear Research. 
Principles of Acceleration to High Energies. The Electron Synchrotron. 
The Synchrocyclotron. Linear Accelerators. The Proton Synchrotron. 
Alternating Gradient Focusing. References. Index. 


1954 5x8 165 pages, with many diagrams 26/- 


THE METHOD OF TRIGONOMETRIC SUMS IN 
THE THEORY OF NUMBERS 


By I. M. VINOGRADOV 


Translated from the Russian, revised and annotated by K. F. Rotu, B.A., Ph.D., 
and ANNE DAVENPORT, M.A. (University College, London) 


** Since 1934 the analytic theory of numbers has been largely 
transformed by the work of Vinogradoy. This work, which 
has led to remarkable new results, is characterized by its 
supreme ingenuity and great power....” 

from the preface. 


CONTENTS: Preface. Notation. Introduction. Note on Vinogradov’s 
Method by the Translators. I. General Lemmas. II. The Investigation on 
the Singular Series in Waring’s Problem. III. The Contribution of the 
Basic Intervals in Waring’s Problem. IV. An Estimate for G(n) in Waring’s 
Problem. V. Approximation by the Fractional Parts of the Values of a 
Polynomial. VI. Estimates for Weyl Sums. VII. The Asymptotic Formula 
in Waring’s Problem. VIII. The Distribution of the Fractional Parts of 
the Values of a Polynomial. IX. Estimates for the Simplest Trigonometrical 
Sums with Primes. X. Goldbach’s Problem. XI. The Distribution of the 
Fractional Parts of the Values of the Function. 

Each Chapter with Notes by the Translators. 


1954 6x9 190 pages 33/- 


INTERSCIENCE PUBLISHERS LTD 
88-90 CHANCERY LANE, LONDON, W.C.2 


XXIX. The Growth of Ice Crystals from the V. apour 


By D. Saw and B. J. Mason 
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[Received November 18, 1954] 


SUMMARY 

The growth and evaporation of single ice crystals have been studied in 
relation to the temperature and supersaturation of the environment. 
The sequence of crystal forms which occurs in clouds as the temperature 
falls from 0°c to —40°c has been reproduced in these experiments. 
Measurements of the growth rates of individual faces of crystals growing 
under fixed environmental conditions reveal that dR?/dt—const. where R 
is a linear crystal dimension, but the growth rates of basal and prism 
faces are different to a degree depending mainly on the temperature. At 
constant temperature, dR?/dt varied as the square of the supersaturation 
o over the range of o studied. A critical supersaturation for growth was 
found which varied in an apparently random manner from face to face, 
there being no systematic difference between the values for the prism 
- and basal faces. The crystal habit, defined by limiting values of the 
axis/diameter ratio, was found to be controlled very largely by the 
temperature, the supersaturation having a much smaller, non-systematic 
effect. 


§ 1. IvTRoDUCTION 
THE growth of ice crystals, besides being an interesting and difficult study 
in crystal physics, is of particular importance to the cloud physicist 
concerned with the development of the ice phase in natural clouds. Snow 
crystals occur in a wide variety of shapes and forms and many workers 
have tried to relate the occurrence of particular crystal forms with the 
meteorological conditions. That they have not found a clearly marked 
correlation between the predominant crystal habit and the temperature 
at the ground is, perhaps, not surprising, since only the conditions prevailing 
during the growth of the crystals are likely to be of major importance. 
It is only in recent years that crystals have been collected from different 
types of clouds having widely different conditions of temperature and 
water-vapour concentration, and that a correlation between the crystal 
habit and the environmental conditions has been established. The most 
comprehensive study of this kind has been made by Weickmann (1947) 
who took samples at different heights up to cirrus levels where the temper- 
ature was below —40°c. His observations are summarized in table 1. 


* Communicated by the Authors. 


S2 


250 D. Shaw and B. J. Mason on the 


The familiar needle-like crystals, which are not represented in Weick- 
mann’s observations, are observed at the ground when the temperature is 
only slightly below freezing, and according to Wall (1947) and to Gold 
and Power (1952), originate from clouds with temperatures between 
—3°c and —8°c. 


Table 1. Weickmann’s Observations of Predominant Crystal Forms in 
Different Cloud Types 


Level of Temp. range | Cloud types Crystal forms 
observation 


Lower 
troposphere | 0°c to —15°c | Nimbostratus |) Thin hexagonal plates — 
Stratocumulus| | Star-shaped crystals show- 
ing dendritic structure 


Stratus 
Middle —15°c to 
troposphere —30°c | Altostratus Thick hexagonal plates 
Altocumulus Prismatic columns—single 
prisms and twins 
Upper <—30°c Isolated cirrus| Clusters of prismatic co- 
troposphere lumns containing funnel- 
shaped cavities. Some 
single hollow prisms 
Cirrostratus Individual, complete prisms 


The existence of a number of basic forms of snow crystals, e.g. hexagonal 
plates, hexagonal prismatic columns, dendritic forms, and of an almost 
infinite number of variations on each of these main themes, suggests 
that their growth and development are complicated matters. As some 
of these variations may arise from changes in environmental conditions as 
the crystal falls through the atmosphere, it seems judicious to study the 
growth of crystals in the laboratory where the temperature and super- 
saturation of the environment can be controlled. 

The growth of crystals in artificially-produced supercooled water clouds 
has been studied by aufm Kampe, Weickmann and Kelly (1951) and by 
Mason (1953). The clouds are formed by introducing steam into thermo- 
statically-controlled, room-size, cold chambers and crystal formation is 
initiated by seeding the cloud with a minute quantity of solid carbon- 
dioxide. In the presence of a water cloud it can be assumed that the 
crystals are growing in an atmosphere saturated with respect to liquid 
water, and therefore supersaturated with respect to ice by an amount 
determined solely by the temperature. 

The observed changes of crystal habit with temperature were very 
similar in both sets of experiments and are summarized in table 2, which 
shows a marked similarity to table 1 giving the temperature classification 
of natural snow crystals. 


i 
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Table 2. Changes of Crystal Habit with Temperature in Artificially- 
produced Water Clouds (aufm Kampe et al., and Mason) 


Temperature range Crystal habit 
0°c to —5°o Single, clear, hexagonal plates with no surface 
; markings. Some trigonal plates 
—4°c to —9°o Prisms, some showing marked cavities and similarity 
to needles 
—10°c to —25°c Hexagonal plates showing ribs, surface markings 


and tendency to sprout at corners. Sector stars 
Dendritic stars, most prominent below —14°c 

—25°c to —40°c Single prisms, twins, and hollow prisms. Aggregates 
of prisms and irregular crystals 


The most striking feature of table 2 is the remarkable sequence of habit, 
plates—prisms—plates—prisms which occurs as the temperature is lowered 
from 0°c to about —25°c. Such marked changes in habit within this small 
temperature range, involving changes from preferred growth along the 
principal c axis for prisms, to preferred growth along the secondary a 
axes for plates, are apparently peculiar to ice. 

To investigate in more detail the factors which control the relative 
growth rates of different crystal faces, it was decided to make continuous 
observations on single crystals growing in an atmosphere where the 
temperature and supersaturation could be controlled and varied in- 
dependently. The crystals were grown on a metal surface, their growth 
being recorded by time-lapse photography so that the rates of advance 
of the different crystal faces could be measured. 


§ 2. THe EXPERIMENTAL ARRANGEMENT 


A diagram of the apparatus is shown in fig. 1. The clean, smooth, 
metal surface on which the crystals are grown is situated in the centre of 
a cylindrical metal chamber, internal diameter 3 inches and depth 2 inches, 
the hollow walls of which may be cooled down to a temperature of —55°c 
by circulating petrol chilled with solid CO, through the annular space. 
The temperature of the chamber walls is controlled by the temperature 
of a cold bath in which is immersed a long copper spiral forming part of 
the circuit for the cooling fluid. The base of the chamber is coated with 
ice, so that the air in the experimental space, which is stirred by a small 
electric fan, is uniformly saturated with respect to ice at a temperature 7’, 
indicated by a sensitive thermocouple embedded in the surface of the ice 
layer on the floor of the chamber. The temperature over the floor does not 
vary by more than 0-2°c, and the thermocouple is located in the coldest 
spot near the inlet for the cooling fluid. . | 

The experimental surface, a metal plate, is supported on a long copper 
rod which is insulated from the chamber and dips into a Dewar ilask 
containing liquid air. The plate is thus cooled to a temperature 7’, 


lower than that of the surrounding air, and which ean be controlled by 
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varying the current through a heating coil wound around the copper rod, 
The temperature 7’, of the growing crystals is recorded by a thermocouple 
soldered chinedintely below the experimental surface, and the saturation 
ratio of the air in the immediate vicinity of the crystals is altered by 
adjusting the temperature difference 7',—7',, and is measured by the 
ratio of the equilibrium vapour pressure of ice at the temperatures 
T , and 7’, respectively. 


Fig. 1 


Pointolite 


Source WD U7 
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diguid Air 


Diagram of the apparatus. 


The crystals were viewed through a metallurgical microscope fitted 
with a vertical illuminator and a 4 mm objective, and were photographed 
at one-minute intervals with a 16 mm camera using a fine-grain film. 
The growth rates of the crystal faces were determined from measurements 
made on the negatives with a micrometer eyepiece. 


§ 3. FoRMATION OF CRYSTALS ON THE METAL SURFACE 

In order to produce ice crystals on the test surface it was generally 
found necessary to cool it to the dew point of the surrounding air, so that 
condensation occurred on the plate. This confirmed the eathes observa- 
tions of Weickmann (1947) that the great majority of crystals form only 
when the air is saturated with respect to liquid water, and therefore 
probably arise by the freezing of the supercooled liquid. Only an 
occasional crystal formed slightly before the dew point was reached ; 
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in these cases the air was highly supersaturated relative to ice and con- 
densation of very small quantities of water may have occurred in surface 
cavities of solid nuclei or in fine scratches on the metal surface. To 
obtain a regular crop of well-formed crystals the metal plate was exposed 
to a small dose of silver-iodide smoke which provides active nuclei at 
temperatures below —5°c. 

Typical examples of crystals forming at different temperatures are 
shown in fig. 2 (Plate 3). Crystals which grew with their c axes normal to the 
surface permitted measurement of the rates of advance of three pairs of 
opposite prism faces. In the case of hexagonal prisms growing with their 
c axes parallel to the surface, measurements could be made of the growth 
rates of both the basal faces and the prism faces. Since it was possible 
to measure only the diameters and length of a crystal, the sum of the 
growth rates of opposite faces, rather than the growth rate of each face 
separately, was determined. 

It was necessary to work with a sparse deposit of crystals, their spacing 
being large compared with their linear dimensions, so that the growth 
of one crystal did not seriously interfere with that of its neighbours. 


§ 4. CORRESPONDENCE BETWEEN CrySTAL Forms APPEARING ON METAL 
PLATE AND IN CLOUDS 


The first experiments were designed to show whether the various crystal 
forms appearing in natural and laboratory-produced clouds as described 
in the Introduction, could be produced in the same temperature ranges 
on the metal plate. 

The plate temperature was held at various temperatures between 
—5°c and —40°c and crystals were grown at constant temperature, the 
supersaturation being adjusted to correspond to water-saturation. The 
crystal habit was found to alternate between prismatic columns and thin 
hexagonal plates as the temperature was lowered according to the 
following scheme : | 

==5°o'to —9°c Prisms, 

= Cito 20 GC Plates, 

below —25°c Prisms. 
As a very few crystals appeared at temperatures above —5°o, crystals 
were initially formed at a lower temperature but grown subsequently 
at temperatures above —5°c. Although they may have originated at a 
temperature normally associated with prismatic columns, when the 
temperature was raised, they assumed subsequently the plate-like habit 
which was typical of the temperature range 0°c to —5°c in the cloud 
experiments. 

Thus, the variation of crystal habit with temperature followed a very 
similar pattern to that found in both natural clouds and the model clouds 
of the laboratory. But, in addition, a new and interesting eryetal form 
quite often appeared on the plate at temperatures between —4°c and 
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—8°c and also below —22°c. This took the form of a prismatic column 
terminated at one end only by a pyramid (see fig. 2). The appearance of 
such hemimorphic forms, which have been found only very oceasionally 
in natural clouds (usually cirrostratus) but not at all in the small-scale 
laboratory clouds, is of considerable interest in connection with the 
possibility of ice possessing a polar lattice (see Owston (1951), Mason 
and Owston (1952)). 
Fig. 3 
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Plot of the square of the linear dimensions of the crystal against time. 


§ 5. GrowTH RaTEs or CrystaL Facks 


The photography at regular time intervals of crystals growing with 
their c axes either normal or parallel to the metal surface allowed the 
rates of advance of both the prism and basal faces to be determined. 
When the square of the linear dimensions ((2c)? or (2a)?) of a crystal 
were plotted against time, straight lines resulted as shown in fig. 3. 

The growth rate of a crystal under steady-state conditions is determined 
by the balance between the rate at which material diffuses towards the 
crystal and the rate at which it can be built into the crystal. If the 
diffusion rate is much greater than the rate of incorporation at the crystal 
face, the growth is said to be surface limited. If, however, the reverse 
is the case, the growth rate is determined almost entirely by the rate of 
diffusion, i.e. we have diffusion-limited growth. Under the latter 
circumstances, the rate of advance of a particular face should be inversely 
proportional to some linear dimension of the crystal. The fact that in 
our experiments the quantities da/dt and dc/dt were inversely proportional 
to a and ¢ respectively and not to c and a respectively (the ratio y=c/a 
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varied considerably during growth), enables us to describe the growth of 
the faces by the equations 
(2c)P?=Ut-+«; (2a)?=U,t+8 
where c and a are respectively the semi-principal and semi-secondary 


crystal axes, U,, U,, « and § are constants. It follows that the crystal 
grows towards a limiting habit which may be defined by 


P= [ela |in=V(UJ0.)=v| 5 (20/5, (2a |. 


The limiting habit, which can be determined experimentally from the 
plots of (2c)? and (2a)? against time, is the most suitable parameter to 
use in discussing the crystal habit because, unlike the actual c/a ratio 
of the crystal, it is independent of the period of growth. 

Fig. 4 
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The growth of three opposite pairs of prism faces of the same crystal, showing 
that different faces grow at different rates. 


Although the square law of growth held for all crystal faces indicating 
- that the growth was largely controlled by diffusion, the actual values of 
da®/dt and dc?/dt were somewhat different for different crystals growing 
under the same environmental conditions, while erystallographically 
similar faces of the same crystal sometimes grew at different rates ea 
ences of 50° were not uncommon) as shown in fig. 4. Furthermore, me 
growth rate of a particular crystal face sometimes.changed abruptly, even 


though the temperature, supersaturation and rate of stirring remained 
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constant. A typical case is illustrated in fig. 5. Very similar behaviour 
has been observed by Bunn (1949) and by Humphreys-Owen (1949) in 
the case of crystals growing from solution. These results indicate that 
the growth of crystal faces is not determined solely by external conditions, 
but is partly influenced by certain properties of the faces themselves. 
Thus the variations in habit may be partly attributed to different faces 
having their own individual characteristics which control the rate at 
which the available molecules can be built into the crystal structure. 


Fig. 5 
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Discontinuities in the growth rate of particular crystal faces. 


Detailed examination of all the available data on crystals growing at 
constant temperature revealed no correlation between the limiting habit 
of a crystal and its rate of mass increase which is a measure of the con- 
temporary total flux of vapour towards the crystal. This point is of 
importance in that Weickmann (1950) and also Marshall and Langleben 
(1954) have suggested that the crystal habit is essentially controlled by 
the flux of vapour towards the growing surfaces. 
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§ 6. VARIATION OF GROWTH AND EvaporaTION Rates with SATURATION 
RATIO 


_ To study how the rate of propagation of a crystal face varied with the 
supersaturation of the environment, observations were made on individual 
crystals growing at constant temperature, as the supersaturation, governed 
by (7,—T,), was changed. The first step was to establish the absolute 
value of the supersaturation for each experiment. To do this, graphs 
of ,/U, and ./U, were plotted against (7,—T',). The resulting straight 
lines cut the (7',—7',) axis at a distance + AT from the origin which was 
the same for the plots of \/U, and 1/U,. The intercept 47 was always 
small in magnitude (<1°c), varied in sign from experiment to experiment 
and showed no systematic correlation with the plate temperature 7’. 
These small fluctuations of the straight-line plots about the origin were 
probably due to slight changes of temperature distribution within the 
experimental chamber and also small changes in the air temperature within 
the cold box from experiment to experiment. These results indicated 
that ideally the supersaturation of o of the air in the immediate vicinity 
of the crystals was given by o=p(7';)/p(7',)—1 where p(7',) and p(T.) 
are the equilibrium vapour pressures of ice at temperatures 7’, and 7’, 
respectively ; the actual supersaturation prevailing in any one experiment 
was assumed to be [p(7, + 47)—p(7')|/p(Z'2). In this way it was possible 
to obtain plots of \/|U,| and ./|U,| against the supersaturation 
c, the temperature 7’, remaining constant, for a large number of crystals 
and for different values of 7',. Typical plots are shown in figs. 6 (a), (5), 
which represent observations made during both growth and evaporation 
of basal and prism faces. In general, it was found that the quantity U was 
proportional to the square of the supersaturation o (the values of o employed 
ranged from 0 to 40%), but with the higher supersaturations there were 
occasions on which the growth rate was smaller than that. which would 
obtained by extrapolation of the U —o? relation. 

The behaviour of crystals evaporating under large subsaturations was 
interesting ; under these conditions of rapid evaporation the reciprocity 
between evaporation and growth suggested by fig. 6 broke down. The 
erystal developed well rounded corners and eventually lost its external 
crystalline appearance, the rate of evaporation being much faster than 
the growth rate under a supersaturation of the same magnitude. It was 
suggested by Frank (1949) that at moderate subsaturations, dissolution 
of a crystal face should proceed in a manner closely equivalent to growth, 
i.e. by unbuilding at molecular steps ending on screw dislocations, but 
that under large subsaturations, evaporation may arise from attack at the 
edges and corners to produce curved faces ; evaporation then proceeds 
rapidly and the crystal loses its geometrical shape. 


§ 7. CRITICAL SATURATION Ratios FoR GRowTH AND EVAPORATION 


For some of the crystals studied, the saturation ratio «=e-+1 had to 
exceed or be less than a certain value before growth or evaporation of the 
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faces could be detected (see fig. 7). 
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The value of the critical saturation 


ratio was often appreciably different for different faces of the same — 
crystal, but there was no systematic difference between the initial resistance — 


Fig. 6 


| 
/|Ua| (10 ém sec 


3.0 

Crystal Temperature 
-17-.8-C 2.0 
1.0 


(2) 


-0:20 -0.10 


=%) 


(er 

t— + —+—__——_—_—_—_+—_> 

0.10 0-20 0.30 0.40 
(a) 


JiUel om sec) 
or lUcl (lOcm? sec") 


Crystal Temperature 


—21-8°C 
Hi 


(b) 


\Ucl 


ent 5 
ieee! 


s+ 
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of prisms and basal faces. This is obvious from fig. 8 in which are 
plotted the critical values of « for crystals growing at different tem- 
peratures. Again, there was evidence that growth and evaporation 
occurred more readily on some faces than on others which were crystallo- 
graphically similar. 


Fig. 7 
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Values of | U,| plotted against supersaturation showing the critical saturation 
ratios required for growth and evaporation of crystal faces. 


§ 8. VARIATION OF CrysTaL Habit WITH TEMPERATURE AND 
SUPERSATURATION 


That the curves of growth rate plotted agaist supersaturation 
(figs. 6 (a), (b)) for the basal and prism faces were, in general, of constant 
slope indicates that the crystal habit defined by P=1/(U,/U,) was 
generally independent of the supersaturation. This was confirmed by 
plotting I against the supersaturation for each of a large number of 
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crystals grown at the same temperature; for any given value of the 
supersaturation, there was a considerable spread in the values of I” for 
different crystals but there was no obvious dependence of J’ on the super- 
saturation. The maximum and minimum values of I’ obtained at various 
fixed temperatures when the supersaturation was varied are shown in 
fig. 9. Remembering that values of "<1 indicate a tendency towards 
plate-like growth and values of [>1, a tendency towards columnar 


Fig. 8 
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The critical saturation ratios for growth of basal and prism faces as a 
function of temperature. 


development, fig. 9 suggests that there is a correlation between the 
temperature and the crystal habit, with a suggestion of the sequence 
plates—prisms—plates—prisms described in §3. This indicates that the 
crystal habit is controlled much more by the temperature than by the 
supersaturation. The manner in which the crystal habit varied with 
temperature was brought out more clearly by confining attention to those 
crystals in each batch which had the largest growth rates (and so were 
probably least influenced by neighbouring crystals), and to regimes of 
supersaturation in which both 1/U, and 4/U, were proportional to o, 
so that J’ was independent of the supersaturation. Walues of I were 
obtained from the slopes of the .»/U—c plots (e.g. fig. 6), and are shown 


—— 
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as a function of temperature in fig. 10. This indicates a tendency towards 
plate-like growth between 0 and —5°o, for prismatic columns between 
—5°c and —9°c, for plates between —9°c and —26°c and for prisms again 
below —26°c. Thus, in reducing the scatter in the observations occasioned 
by the individual peculiarities of each crystal, we have shown that the 


Fig. 9 
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fixed temperatures when the supersaturation was varied. 
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XXX. Theory of Ionization Fluctuations 


By J. E. Moya* 
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SUMMARY 

The distributions of the loss of energy by ionization of a fast primary 
and of the numbers of ion pairs it produces are derived. It is shown that 
down to quite small values of the primary ionization, both can be 
represented by the same ‘ universal’ distribution if the variables are 
reduced by a proper choice of scale and origin, which accounts for the 
experimental fact that ion pair numbers are proportional to primary 
energy loss. These conclusions remain valid when one takes into account 
quantum resonance effects and the details of atomic structure of the 
absorber. 


§ 1. [INTRODUCTION 

THE main object of the present paper is to derive expressions for the 
distributions of (a) the loss of energy by ionization of a fast primary 
particle passing through an absorbing medium; (b) the numbers of ion- 
pairs produced by such a particle. We shall try in particular to explain 
the experimental fact that these two distributions are very approximately 
proportional to each other: i.e. that the energy lost by the primary 
particle per ion-pair produced is approximately constant (of the order 
of 35 ev). 

Landau (1944) has shown that under certain simplifying assumptions, 
the distribution of energy loss can be expressed as a universal curve in 
terms of certain reduced energy variables (depending on the charge, 
mass and velocity of the primary, the atomic properties and the thickness 
of the absorber). We shall first derive a closed analytic expression for 
Landau’s distribution. In view of discrepancies between this distribution 
and the results of recent experimental work (see West 1953) we shall 
examine possible departures from it due to: (a) small thicknesses of 
absorber; (b) the influence of the detailed atomic structure of the 
absorber; (c) the influence of quantum resonance effects in distant 
collisions. Finally, we shall derive a theoretical expression for the 
distribution of the numbers of ion pairs produced by the primary. This 
will allow us to decide whether these discrepancies can be due to the usual 
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assumption made in cloud chamber, ionization chamber and propor- 
tional counter work of a constant energy loss per ion-pair. The experi- 
mental evidence will be analysed in a forthcoming paper by Owen and 
Eyeions; a review of the subject by Price will appear in Reports on 
Progress in Physics. 


§ 2. GENERAL EXPRESSION FOR THE 1ONIZATION ENERGY 
DISTRIBUTION 


In deriving an expression for the ionization energy distribution, we 
shall make the following assumptions : 

(1) Successive ionizing collisions are statistically independent. 

(2) The total energy loss of the primary is very much smaller than its 
initial energy, and hence the decrease of primary energy may be neglected. 
Our theory will not therefore apply to slow primaries or great thicknesses 
of absorber. We also assume that energy losses due to radiation or 
nuclear interactions are negligible. 

(3) The absorbing medium has a homogeneous constitution. 

We shall develop the theory first for an arbitrary total cross section 
o(H#), where No(EH) dE dt is the probability of an energy loss between 
& and H+-dE in an absorber thickness dt, N being the number of absorber 
atoms per unit volume. The primary ionization rate is then 


fees , o(E) dE, 


and the ionization energy distribution per collision is 
No(H Pie 
d(h)= ae with | d(E) dH=1 
“0 
(ydt is the probability of an ionizing collision in the thickness dt, 
$(H) dH the probability of an energy loss between # and H+dE given 
that a collision has occurred). 
Let x(L , t) dE be the probability of an energy loss between # and H+-dE. 
in a thickness f. It is easily seen to follow from assumptions (2.1), 
(2.2) and (2.3) that 
Kh 
X(H#, ty +t.)= | x(H—W, t,)x(W,t)dW. . . . (2.1) 
“0 
For a small thickness df, the probability that no energy is lost is 
1—qédt-+-o(dt); the probability of an energy loss between HF and H+dK 
is gf(H) dHSt4-o(dt) ; hence 


x(L, 8t)=(1—gSt)8(Z)+9¢(E)St+o(8t); =... (2.2) 
in particular, y(H#, 0)=6(£).* 
We introduce the Laplace transform of y(Z, t) : 


r co 


M(a,t)= | exp(—aH)y(H,t)dH. . . . , (2.3) 


“0 


*'The Dirac 3-function 8(£) should not be confused with St. 
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Tt follows from (2.1) that 
i MM (Gxt ts) M (a, t,)M(x,t),- . « 2 . , (2:4) 
and hence from (2.2) and (2.4) (noting that M(«, 0)=1) that 


aM(a,t) 


at % 5 Lt ot)—1}M (a, t) 


ere face 
= Lim 54 (1a) + 48! | . exp (—al)4(£) 1B -+-0(3)—-1} M(, t) 
= {a [exp (a8) —1y9(H) al at, 0 bd Mas Fhe 


The solution of the differential eqn. (2.5) with initial condition M (z, O=1 
is 


Mia, exp {at | [exp (—2#)—114(H) d| exp [QR], 2.6) 


where Y=gqt is the mean number of collisions (i.e. the primary ionization) 
in the thickness #, and 


Ria) | jexps( 4h) = lit (B) og hina ei ne (2-3) 
The standard inversion formula for Laplace transforms then gives 
ili C+io ; 
x(£, Q)= sa) _ exp [OR (z)--2E] dz » ieee (28) 


(it is convenient to express y in terms of the mean collision number @ 
instead of the thickness 1). 

‘An asymptotic expansion for y may be found from (2.8) by the saddle 
point method (see e.g. Jeffreys 1946 and Daniels 1954). As is usual in 
applications of this method, we can accept as a sufficient approximation 
to x the first term of this expansion, namely, the expression 


| s 
x(L, Q) dE= F [27QR’(x)|-1? exp {Q[R(a)—aR'(«)]} d#, (2.9) 
where « is related to H by the expression 


E=—QR'(«)=Q | 2 exp (—al)¢(F) dk, ye (I 


> CO 


Ri(aj=| Hexp (—aB)g¢(H) dE, . . . . (2.11) 


{ 
| 
hs 
and ¢ is a normalization constant, chosen to make 
7 CO 
[ x(8, @) dB=1. 
0 


The most probable energy loss H, is by definition the value of E fcr 
which xy(#, Q) is maximum. Maximizing the right-hand side of (2.9) 


with respect to «, we see that H p= QR" (%), where «, is the solution of 
Ri (0) + 204) B" (tp) = 9. ‘ me. ee, (el) 


TZ 
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§ 3. Tae CiassicaL (THOMSON) CROSS SECTION 
Following Landau (loc. cit.) we now take for o(#) the classical 
(Thomson) cross section (see e.g. Bohr 1948) 
{ BZ dE/E? for EH ,<H<E,,, 
No(E) dE= 2 ; ee ieee vs 
| 0 3) (Reh, f On heb 
22 et 
pvr” 


(3.1) 


where B=2nN (3.2) 


N is the number of atoms per unit volume in the absorber, « the mass 
of the electron, « its charge, ze the charge of the primary, v its velocity, 
Z the ‘ effective ’ atomic number of the absorber; Hy and E,, represent 
average values for respectively the minimum and maximum transferable 
energy per electron in an ionizing collision (see §5). It follows that 


im dH rich 1 BZ uae 
y=BZ Ie mE == BZ a — r) = TF and hence Q= EB,’ 
dE 


This cross section gives the ‘ classical ’ expression for the average energy 
loss in thickness f 
Dike o} “KH 
EK=Nt\ Ko(F)dE=BZt ies = BZt eee 5 ee SE 
Jo Jz, & By 
known to be too small by approximately a factor of 2. The integral 
diverges if we make H,, ~ co. However. it is found that for fast primaries 
the value of H,,, has a negligible effect on the shape of the energy loss 
distribution curve ;* this renders legitimate the simplifying assumption 
that H,,, is infinite. 
It is convenient when substituting from (3.1) into the expressions of 
§2 to change to the energy variable «-=H/E,; we then find that 


fr oO 


R(«)=exp (—a«)—1—« | exp |—«e](de/e) ; 


1 
7 A, =O =O 
R(w)=—| exp (—ae) 25 R= RM = — (14 *) 
(3.5) 
Hence x(e, Q) de= - (x5) exp [$a 1+Q(e"—l1)] de, . . .(3.6) 
O\iT 2 da 
here > ce xn f—1 i, ae 
where ( JB) I exp |[—t«+Q(e 1)] We 
0 
=exp (—Q) 2 —-——— 912, gi ti ne 


nao V(22-+1) n! 


* j ‘Tavey p jata ‘ x ef 
This curve (see fig. 3) consists of a sharp peak centered on FH, (representing 
the statistical effect of frequent collisions with small energy transfer) .followed 
by a long tail of small ordinate (representing the effect of rare violent collisions) 
r . . . 5 
aking the finiteness of #’,, into account merely alters this tail in such a way as 


to make #< 00 without appreciably affecting the main part of the curve. 
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The integral above is taken from 0 to oo because « decreases from oo to 0 
as « increases oo 0 to co. This can be seen from relation (2.9) which 
connects « to « 

€ , fe ae are 

eee Gree } e-7 ie tt); Or big (O58) 
the exponential integral Ei(x) is tabulated (cf. Janke-Emde 1938 or 
British Association Tables 1951). For small values of a, however, 
(<1), it is more convenient for computation purposes to transform 
(3.8) as follows : 


S < dx: - dx (ce ax 
—— e- 7 —— =| es SE eae | ae al 
@ 1 Os i ) 4. ) og « 
= —C—log «— (ney 
0 
= —C— log « x (—1)” a (3.8 a) 
Re n.n! e 


where C=0-577 is Kuler’s constant. The numerical evaluation of x is 
discussed in the next section in terms of corrections for small Q to 
Landau’s universal distribution. 
Substituting from the above in (2.12), we find the value of «,, of « 
which maximizes (3.6) by solving 
2Q exp (—a,)=14 =, Coe ere Sahin) 
a 
graphically (fig. 1), and the corresponding most probable value of the 
energy loss <,, from (3.8); ¢, is shown as a function of @ in fig. 2. It is 
seen that the equation has no solution for Q<2-44, indicating that for 
such low values of Q, y(e,@) does not go through a maximum, but 
decreases more or less exponentially as e increases. However, the 
accuracy of the saddle-point approximation is poor for such low Q; 
for Q>5. we see that %,<0-13; hence we can approximate to (3.9) by 


the equation 1 
: =(142)=5-1 Spe ean (StL) 
2Q a 5 Oy 
whose solution is to a good approximation 
1 1 
es) =D : 3.11 
x =20) = 2Q—1 ( ) 


p 
Taking the first three terms in the right-hand side of (3.9) we obtain for 
the most probable energy loss the agen expression 


eo lon (2022 ee eee | 05.12 


which for Q=>5 departs by less than 1%, from the value obtained by the 


graphical solution of (3.9) (cf. fig. 2). 
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Fig. | 


LSet ie ia Sp LT AL oaks i Cea, | (Ee STs 
O00! O005 CO! OOS OF os ee) 
Solution of eqn. (3.9) : the values of «, are given by the abscissae of the ee 
sections of the curve 1--1/a with the curves 2Q e-* for values of Q=2-44, 
5, 10, 20, 50, 100; there is no intersection (and hence no maximum) 
for Q<2-44. 


Fig. 2 


100-0 S000 10000 


Oa 


The most probable energy loss #,, (in units Q/,) as a function of the primary 
ionization Q: the full curve is obtained by numerical solution of 
eqn. (3.9), the dotted and dashed curves represent: i Se at the 
approximations of eqn. (3.12), and (4.4), 


———————— 
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a ae ee oe § 4. Tor Lanpau APPROXIMATION 

One finds that for large Q, the major contribution to x(e, Y) in (3.6) 
comes from small values of « (as we have seen, «, is small for Q>5 
and the curve for y is sharply centred about a Hence, following 
Landau, we obtain an asymptotic expression for y valid for large Q by 
neglecting terms of order « compared to log «: thus, we take 


Hilcje=e( C= Ise lomnye’ sue. keen, 2(4,1) 
Substituting these values, (3.6) and (3.8 a) become 


x(e, Q) de= : J (ss) exp (—Qx) de, 


il ag Ode. 1 
where c= exp L (4.4 
Vm) XP) 7a = 578" a 
g =~ C—log Rees Wat ets se (43) 
Yo the same order of approximation, we may take for the most probable 
values 
r | Sy 2. a 
Coes 0 and Q =p C+log2Q. . . (4.4) 
Changing to the reduced energy variable 
pa aK Jib iy 
ee a P =—log2Qu, . .. . (4.5) 


pene OMe BLE 


we find the following explicit expression for Landau’s distribution 


1 
Xz(w) dw= Van) exp {—4[w+exp (—w)]}dw. . . (4.6) 


The half-width Aw of this distribution is easily found to be dw=3-58. 

We see thus that the asymptotic expression for y has a ‘ universal ’ 
form independent of Q when the energy is expressed in terms of the re- 
duced variable w. The accuracy of this expression has been assessed by 
computing x for small values of @ from the more accurate expression 
(3.6) and changing over to the reduced variable (4.5) : the results are shown 
in fig. 3. Surprisingly, no appreciable departure is found from Landau’s 
distribution for Q>20; the departure for Q=10 and Q=5 is indicated 
on the figure by dotted lines. We must mention that the expressions 
for Q and H, above are not the same as Landau’s: the relativistic log 
rise with increasing energy is absent, because we have neglected quantum 
resonance effects (see the discussion at the end of § 6). 


is 
§ 5, ErrEctTs OF ATomMIC STRUCTURE 


__ In order to assess the detailed effects of the atomic structure of the 
absorber on the distribution of ionization energy, we shall use a crude 
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classical model; we assume that the classical cross section is valid for 
the electrons in each shell of the atom, with a minimum transferable energy 


Fig. 3 
\ 
\ 
\ 
\ 
\ 
| gry o rae 
40 (are) 20 100 
> 


Distribution of primary energy loss: the full curve is Landau’s ‘ universal ’ 
distribution (4.6), the dashed curve is computed from eqns. (3.6) and 
(3.8), for Q=5 and 10. 


equal to the ionization potential uf the shell. Let 2; be the ionization 


potential of the jth electron, j7=1, 2,..., Z; then the total cross section 
is z 
\ Th 5 B/E* for ESI, 
o(H#) dBE= KH) dB, with Hy} ee 
(L) sos o;(H) ¢ with No,(#)= ‘i ; B<1,, (5.1) 
where B has the same definition as in eqn. (3.2). ‘Then 
ie pl Ea F4, 
@=Nt| o(f)dE=Bt 2 =Bi 2) == , 
I ) any: Bil >= (5.2) 
where 4 tie ‘ i 5.3 
Le ae ey 2) 
and if 5 
ee SNe 7 
ne d(H) dE= ZN te Gs( 20 LE ree 5 5 aati Oca) 


: 14g, FAL nl 
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Changing to the energy variable «=H/E,, and substituting in the 
expressions of §2, we find R(«) and hence x(e, MY). It turns out that an 
asymptotic expression for y similar to that derived in $4 is valid for 
values of Q>20. Restricting our attention to this case, we find that to 
the same order of approximation as in § 4, 


e)x| C1 tox (2) |, Wists wre nt (OcD) 


a 


Z 
‘where log #,=(1/Z) Slog /;. (Note that H=BZt log (H,,/E,); ie. in 
j=1 


m/ 


J 
the present classical model, H,, and #, are the average ionization potentials 
appropriate for the calculation of respectively the primary ionization Q 
and the average energy loss #.) We then find for the most probable 
values «,,’=1/2Q, and hence 


& _ #,' eas 2 lie 
GAB CE 
ie 2BAt ; 
ae (| E, (5.6) 
Introducing the reduced energy variable 
—<, LF’ 
<a Ae Bes -log 200. ae ieee ee Oe 


Cana BIE 


we see at once that the asymptotic distribution is exactly the same as 
_ Landau’s, namely 


1 ; 
X1(w) dw= (2m) exp {—4[w+exp(—w)]}dw. . . (5.8) 


Thus, the details of atomic structure, with the assumptions made above, 
have no effect on the asymptotic distribution* but only modify the 
expression for the average number of collisions Y and the most probable 
energy loss #,,; these conclusions are not likely to be modified by any 
more refined theory of these effects, though we may expect improvements 
in the expressions for @ and £,,. 


§ 6. Quantum RESONANCE HKFFECTS 


The classical cross section proportional to 1/H? is valid for values 
of E>I, the average ionization potential of the absorber. Its failure 
in case <2/fiv—c/137v <1 (where fi is 47 times Planck’s constant) may be 
ascribed to quantum resonance effects in distant collisions, which increase 
the cross section for values of # near J, and consequently approximately 
double the mean energy loss (see e.g. Bohr, /oc. cit., p. 89 for a discussion). 
Tn order to assess the effect of this resonant increase on the ionization 


ots SMES el ena en am ne 
* Blunk et al. (1950, 1951) appear to arrive at a contrary conclusion. 
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energy distribution, we shall postulate (since an exact expression is not 
available) a cross section 
BZ dE 
o(H#) dE= 


N (8#-1-r?’ 


of the usual form for quantum resonance effects, /” being the half-width 
of the resonance curve. We expect that <J; hence o is approximately 
equal to the classical cross section when H>J. In calculating the ioniza- 
tion energy distribution, we must set o(#)=0 for non-ionizing collisions 
where H<J. Hence ; 


(6.1) 


. re dE 2 pS oda eee BZE MM 
Q=Nt] o(H) d= BZ | Hae = BH | pr 2a 
- SOF * de [ 
ne Ree ee } 

d(H) dEh= 7 (HTP ET | 


Transforming to the energy variable «=(H—1)/J° and substituting m 
the expressions of § 2, we find that 


2752 dle 2 


R(a)= =a [exp (—ae)—1] [eae 1+ = (sin « Ci «—cos « si «), (6.3) 


where Ci «, si« are the cosine and sine integrals respectively (cf. Jahnke- 
Emde, loc. cit., p. 3). For the same reasons as in §5, we need only consider 
the asymptotic distribution. It is easily seen from the series expansions 
of Ci «, si« that to the same order of approximation as in § 4, 


2 
a 


R(a)= < (C—I-plogia)y,e vo) Se ee a Oras 
' 2#H,'’—I 2 4 
ty = 7 , and hence a Se BZA ie (a a -0) Beat (8 


Hence we see that as in §5, if we change over to the reduced energy 


variable iy ere 2 ales a es 
then the asymptotic energy distribution takes the ‘ universal ’ form 
| . 
X1(w) dwo= Vn) exp {—3lw+exp (—w)]}dw. . .° (6.7) 


Thus quantum resonance effects too modify only the expressions for @ and 
£,,, but not that for y,(w), and here again a more exact theory of these 
effects is unlikely to change these conclusions (though we may expect a 
modification of the expressions (6.2) and (6.7) for Q@ and #£,). An 
estimate of J” may be obtained by equating the primary ionization rate q 
with the value calculated by Bethe (1933) 


BZ 2.028 2B 
= — — —p2s =" _ 
ads ( r{log 204 B*} mae 


Ee 
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‘where B=v/c, [, is the outer shell ionization potential, r and s are constants 
depending on the absorber (for hydrogen, r=0-285 and s=3-04). With 
this estimate, Q and #, exhibit a relativistic log rise with increasing 
energy similar to the expressions given by Landau (loc. cit.). Alter- 
natively, one may estimate I" from the experimental values found for q. 


§7. Tae Distrisution or Numbers oF [on Patrs 


_ The production of ions by a fast primary particle constitutes a multi- 
plicative or branching stochastic process: the electron ejected in a primary 
ionizing collision may, if it has enough energy, ionize another atom; the 
secondary electrons may then ionize again, and so on, the process ending 
when all secondaries have become too slow to ionize further. We shall 
develop the theory under the assumption that the collection or recording 
of ions is delayed long enough for the process to terminate, so that the 
ion pair distribution per collision is the same for all collisions: this 
assumption is valid as a rule in ionization or cloud chamber work: post- 
expansion cloud chamber tracks form however an exception, since they 
record only the primary ionization. We shall also neglect extraneous 
effects such as increase in ionization due to acceleration of the electrons 
in electrode fields and to photoelectric effects in the gas or chamber walls, 
which become important in proportional counter work. 

Let then gq, be the probability that / ion pairs are produced in any 
given ionizing collision, with s=1, 2,... and 2q,—1; let p,(t) be the 
probability that a total number n of ion pairs is produced in an absorber 
thickness ¢, and write g as before for the primary ionizing collision rate. 
It then follows from assumptions (2.1), (2.2) and (2.3) that 


Paltitte= > Pil) Pulte), ee ays, Om pee he) 
= 
Pnlt +St)=(1—95t)8,,9+97,5t+0(5t). . . - ~ (7.2) 
Introducing the Laplace transform of the ion pair distribution p,,(t) 
Vo ya CH CU a(t), © Mile ee oe Cua) 
n=0 
we find as in § 2 that 
M (a, t,t.) =U (a, #1) . M(«, te), em oem S) 
CMe) a G z (o-™—1)t4} (at eee tO) 
ot n=1 
and hence that 
M(«, t)=exp {at 2 (e-"=1)4} =exp [QR(«)], . . (7.6) 
n=1 


where Q=gqt as before, and 
Waste Rees fea} 


R(x)= Blen™—Wfy ee + (17) 


n=l 
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The inversion of the Laplace transform (7.6) is given by the well-known 
formula 


af te 


"exp [QOR(z2)-+-nzZ| dz... ss wien 


L 
PAQ)= oe 
This integral can be evaluated by the saddle-point method : the first 
approximation to it (which may be assumed as previously sufficiently 


accurate for our purpose) is 


| / 
P,AQ)= Cc [27QR''(a,,)) se: exp {Q[R(a,)—«,R (~,)]}; . (7.9) 
where «,, is related to n by 
n= =O (4, )=Q 2 b exp (kas) Gg, 2 = = tee 
k=1 


ive) 
R'(a)= 2 k? exp (—ka) q,, 
k=] 
and ¢ is as before a normalization constant 


c= y [27QR'(«,)) 1% exp {Q[R(«,,)—a,,F’(a,,)]} 
r=) 


Uy 


ae 
me |” IQR (a) exp {Q1R(a)—a.B'(a)}} dex + Spelt). (7-11) 
V (277) Ja, 
The last approximate expression for c is obtained by the use of Euler’s 
summation formula, with x9, «,, corresponding respectively by (7.10) to 
n=0 and n=0o. These expressions are close analogues of those found 
in §2. Similarly, the most probable number of ion pairs is, as in $2, 
the number x, (which need no longer be an integer) corresponding by 
(7.10) to the solution o of 


Ri(o5) 2Qa, RY (o,)P=0. 92 eee eee 


§ 8. lon Pairk DisTRIBUTION WITH THE 
CLASSICAL Cross SECTION 

We shall now develop the theory of ion pair distribution for the 
following crude model: we assume the classical cross section for ionization 
by the secondary electrons as well as by the primary particle, and set 
this cross section equal to zero if the energy of the electron is less than the 
mean ionization potential / of the absorber. Let EL, be the energy of 
the primary particle, #,, /,,... those of the secondary electrons succes- 
sively ejected in the chain process resulting from a single primary ionizing 
collision. The total cross section o(£ ; H,) dE for an ionization energy 
loss between H and H+d# by the kth electron of energy FE, is 


Gal if E> 1 E<E 
= I Sel ean ad ie 
off; B,)= 4) \#i/ ; : 


hoe ifH,<I or E>E,: 


(8.1) 
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hence, its collision probability per unit thickness is 


By 5 kA Oa 
| oP: By) d= Ga 


i 20 if B<,I, 


while the probability that it should not ionize further in a thickness ¢ is 


| ) if ELSI, 
NE y)= 4 i (8.2) 


vexp{ qretN ( i ye ft El 

: x FT; 1 ae) 

(B,, t)=exp [—A(E,)t]= 4 Ae a ae ES ey 
] if, Hl 


The probability that the primary collision produce a first electron of 
energy between #, and H,--dEH, is IdH,/(H,4+-1)?: hence, the probability 
that just one ion pair be produced in a thickness ¢ is 

a IdE Te 

ih (By, ‘Op =4+ —gy, [l—exp (eA Y/21)]. . (8.4) 
The second term in the right-hand side of (8.4) is negligible for any 
reasonably large value of ¢ (if e.g. /=25 ev then for a monoatomic gas 
awe*N/I? = 28108 cm-!); hence the probability q, of exactly one ion 
pair per collision is just the probability that the first electron has an 
energy <J: i.e. q;=}. Similarly, we can simplify the calculation of 
dos 93,--- by setting «(H,t)=0 for H>TI and letting t— oo in the final 
result. 

The calculations get progressively more complicated for the succeeding 

qd, details for ¢, will be found in the Appendix. However, the distri- 
bution a 


i= ET? 


(e211 Os eS 


(where a is a normalization constant chosen to make 2q,=1) fits approxi- 
mately the first few q,, and has the right sort of behaviour for large k 
(analogous to the 1/H? law for the distribution of energy loss per collision). 
The value of a may be obtained from the formula (see e.g. Knopp 1928) 


cath a —4--1-077; hence a==0:929. (8.6) 


The distribution q, is independent of /: this gives us some confidence 
in its approximate validity in spite of the neglection of inner shell contri- 
butions. In fact, it appears from a rough estimate of the latter’s effect, 
using the model of § 5, that the ¢, are not appreciably modified, provided 
that the inner shell ionization potentials are much greater than that of 
the outer shell. 

Substitution of (8.9) in the expressions of §7 gives 


ar o pki Ied 
il 00 e- kx] if 2 ke heox i 5 ke 
— i= DD» 5 ve ao) Di, 3 —-R (a)= ond a ot nia : 
Ae (ee ee a per eel & ped 
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The series above converge very slowly for small values of «; for the 
purposes of computation, we transform the last two as follows: 


he = eho > oe x eis +log (1—e°%) | 
a () = * We+l) pn he en i )) os , (8.8) 
l RP ()ay E ots Dr ee ee 
Ph techs at meat pee em 1 


For values of «<0-04, it is more convenient to compute by numerical 
integration from the following transformed expression : 


FiRe)+ t= {ap 
=cos a4 E faa — fae} +sin « 
x {log (1—e-*) 4 Eaerh S (SS e. (8.9) 
~ BY(a)= = ia seed th —sin « 2 aa ae = as} | 
| 
| 


co fi r*1—cos 0 
So (eck re ———d@>. 
#00824 log( On ee gatry pmo") } 


— 


The substitution of these expressions in (7.9) and (7.10), taking 
@=BZt/I, yields the ion pair distribution p,(Q): in order to compare 
it with that of the ionization energy, it is important to evaluate its 
asymptotic expression for large Q. It is easily seen from the trans- 
formed expressions (8.9) that to the same order of approximation as 
in § 4, ; 

| n 


R(«)=ax(b—1+log a); «)’= 30a} Oa =—b+log 2Qa; . (8.10) 


where a is defined by (8.6) and b=21/k(k?4-1)=0-672. Substituting in 
(7.9) and changing to the reduced variable 
Cina (2 


hee acral? = —log 20aa, 0 soe ee ney 


we are led again, as in the case of ionization energy, to Landau’s 
‘universal ’ distribution: that is, we find that 


1 N—Ny’ 
P»(Q)= Qa XL (“a) > 


1 é 
where X1(w)= Van) oP {—t[w+exp (—w)]}. 
2} 


— 


(8.12) 


| 
] 
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The deviation from Landau’s distribution for small values of Q is 
somewhat more marked than is the case for the energy loss distribution : 
it is shown in fig. 4, where Qap,(Q), computed from expressions (8.8) 
and (8.9) for Q=20, is plotted against w. 


Fig. 4 


—— an Q =50 


a 


) 2 ir 


nas - — 


Distribution of numbers of ion pairs: the full curve is Landau’s distribution 
(8.12), curve (i) is computed with the classical cross section, §8 for 
Q=20, curves (ii) and (iii) with the quantum resonance cross section, 
§ 9, for Q=50 and respectively g=0-17 and 0-50. 


§ 9. Ion Parr DisrRiBUTION wiTH QUANTUM RESONANCE EFFECTS 

The effect of quantum resonance is to increase the primary energy 
loss cross section for values of the energy near the ionization potential J, 
and hence to increase the probability g, of producing a single ion pair 
per collision. Repeating the calculations of §8 with the cross section 
(6.1) for the primary, and the classical cross section for the (slow) ejected 
electrons, we find that 


PALES rH O} 2 / 
a= — | elo aces me ee (9.1) 
vin 0) Wye [P= TT e 
For relativistic primaries, where I/I’ is large, we have approximately 


(9.1 a) 


278 J. E. Moyal on the 


Continuing the calculations in the same way, one finds that the distri- 
bution q,, can be fitted approximately by the relation 
i 


a + NE 
= ———— = ————_. — = 1-190 — = eae ans 9.2 
ay (geri): where g n(ieayd 1-196 7 (eee ae )> say 


one can easily verify that 2,g,—1. This may be interpreted by saying 
that at every collision there is a chance 1—g of a resonant collision, and 
a chance g of a non-resonant one, given which the distribution of ion pairs 
is the same as in §8. 

Let us now write R,(«), R,’(%), ete. for the expressions found in §8. 
and R,(«), Ry’(«), etc. for those corresponding to the quantum resonance 
cross section. It then follows from the foregoing that 


Ro(w)=(1—g)(e*—1)+9R(a);  Ro'(x)=—(U—g)e*+gR,"(«); (9-3) 


and so on. To the same order of approximation as in § 4, we then have 


Ro(«)=—(1—g)a+-gax(b—1+-log «) ; ] 
sn) aes ae Oe 
Ro’ («)=—(1—g)+ga(b+-log «); ete. } 


and hence 


etl Ae ae | 
hy fo 2Qga ar 2.21Q T°? Ny =I —g+ga(log 2Qga—b)], : (9.5) 


where Q=7BbZt/2l. Substituting in (7.9) in terms of the reduced 
variable ’ 
peep log 2Qga« (9.6) 
CO Se “= Ky ° . . 5 . a. 
Qga & au 


we find as in $8 that 


l n—n,,’ 
)\ = ieee 
p, (Q) aga ( Qga i 


where X,(@)= 


! 
| 
(9.7) 
Van) exp{— }|w-+exp (—w)]}. | 
The deviation from Landau’s distribution is more pronounced than is 
the case with the ion pair distribution found in §8, but is still not very 
marked. This may be seen in fig. 4, where the curves of Qgap,(Q) as 
a function of w, computed from expressions (9.3), are shown for Q—50, 
y=0-5 and g=0-17. We may remark that if J” is evaluated by equating 
q=7BZ/2I" to the experimental values of primary ionization (ch s'G). 
one finds that near minimum ionization g is approximately equal to 4 for 
most gases. ; 
The tendency of the ion pair distribution to the ‘same (Landau) 
distribution as the primary energy loss accounts (at least in the case of 
fast primaries) for the experimental fact that the energy lost per ion pair 
is a constant, approximately independent of the primary energy : 
furthermore, we may presume that this result is not sensitive to the 
actual value of the cross section, since it holds for both the classical 
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(§8) and the quantum resonance ($9) cross sections. Comparing the 
expressions for w in §§8 and 9 with those respectively in §$4 and 6, we 
see that the energy per ion pair must be approximately equal to the 
_ average ionization potential. This is substantially less than the experi- 
mental values (of the order of 35 ev), presumably because the latter 
refers to total energy loss, including excitation of the atoms, whereas our 
calculations refer to energy lost by ionization only. Furthermore, the 
use of a single average ionization potential is undoubtedly a very rough 
approximation, especially at the higher atomic numbers, where inner shell 
ionization and excitation, Auger effects and the like must contribute 
substantially to ionization and primary energy loss. We cannot in fact 
expect the theory in the present section to have more than qualitative 
validity for elements heavier than helium. 


§ 10. CONCLUSIONS 
_ We conclude from the foregoing : 


(1) That Landau’s ‘ universal’ distribution appears to be valid for 
both primary energy loss and numbers of ion pairs, which explains their 
proportionality in the case of fast primaries. 

(2) That it remains valid down to unexpectedly small values of the 
primary ionization Q. 

(3) That none of the effects considered in this paper (atomic structure, 
quantum resonance, proportionality of numbers of ion pairs to primary 
energy loss) explain the experimentally found deviations from it. 
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APPENDIX 
Calculation of > 
In order to evaluate qo, we note first that the probability that electron 1, 
of energy #,>J, suffer an ionizing collision between thicknesses 7 and 
z+dr, thereby producing electrons 2 and 3 with energies respectively 
between H,, H,+dH, and E,, H,+dH;, where H,=H,—EH,—J, and that 
electrons 2 and 3 ionize no further between thicknesses 7 and f, is 
x(H,7) . No(H,; H,) dE, dr. «(H,, t—7) . x(H,—H,—I, t—7), (A1) 
where o and « are given respectively by eqns. (8.1) and (8.3). We then 
obtain g, by multiplying (A1) by the probability distribution Jd#,/(£,+)? — 
of electron 1, integrating with respect to H#, (from 0 to #,—J/), to Hy, 
(from I to 00), to 7 (from 0 to ¢), and finally making t— oo: thus 


j ro 6d t FI 
qe=lim | Ga | (Ba 7) dr) © o(Bly; By)«(E'y, (7) 
(eT tebe Ft ae 


The only positive contribution to g, in the above comes from the range 
of values of H, and H, which make #, and H,—H,—EH,—I both <I, 
and hence /</,<3/. If #,>3/, then one at least of electrons 2 and 3 
will have an energy >J, and will therefore ionize further. If /<#,<2/, 
then #,+H,;<J, and therefore both #, and #, are <J/, and will not 
ionize ; hence, we get from this range the contribution 
pe LI, 1 
i (Ei pL ar6: 
If 27 <E,<3/, then the requirement that #,<J and H,—H,—I<J means 
that the range of #, must be limited to H—21<H,<J. The last two 
. integrals in the right-hand side of (A2) yield in this range 
i ea e. me§N dH, E,(3I—£,) 
J, ePIME ede} E Eyt IPO Bo 
integrating with respect to #,, we obtain the contribution 
[oe TO ot in 3 
Wer DUH yt (Eee a ee Oe 
Adding these two contributions, we finally obtain 


1 3 


Note added in proof.—Contrary to the conclusions stated above, it now 
appears from a critical analysis by Dr. E. P. George of the available experi- 
mental data that agreement with the theory is substantially improved, at 
least for gases, by the results of § 6 (owing to the broadening of the scale for 
the reduced distribution by a factor 7/2: cf. eqn. (6.6)) ; this was also pointed 
out to me in a letter by Dr. B. T. Price. It now seems possible that, with the 
introduction of a more exact cross section, the theory should account for all 
the facts. 

Dr. U. Fano has pointed out that the approximate distribution (4.6) does 
not agree with Landau’s numerical evaluation for large w owing to the break- 
down of the saddle-point method in this range. This is not very important 
as regards comparison with experiment, because the frequency of events with 
energy loss greater than twice the probable value is small (of the order of 5%) 
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ABSTRACT 


The scattering of a parallel beam of unpolarized sunlight in a plane- 
stratified atmosphere is considered, taking into account the polarization 
of scattered light and allowing for the optical anisotropy of the molecules 
of air. Allowance is also made for the attenuation of the ultra-violet 
light due to absorption by atmospheric ozone. Expressions are derived 
for the intensities of the primary, secondary, and tertiary scattered light 
emitted and received in any direction and at any point in the atmosphere. 


§1. [yrropuctIon 


In this paper we consider part of the problem of the ultra-violet ‘ illumin- 
ation ’ of the earth’s atmosphere and surface by scattered sunlight. Of 
special interest is the ultra-violet radiation which is received at the earth’s 
surface from various directions and for various positions of the sun in the 
sky. 

We confine our attention to wavelengths for which the absorption 
coefficients of ozone may be assumed to be independent of atmospheric 
pressure and temperature. The problem is also restricted to a plane- 
stratified atmosphere, assumed to rest on a flat surface. Reflection at 
the ground and refraction are neglected, as also are the scattering effects 
of dust and haze. It is assumed that energy is not transformed from one 
wavelength to another by absorption and re-emission. ; . 

A similar restricted problem, but referring to the scattering of visible 
light in an ozone-free atmosphere, has been treated by Hammad (1948) 
and Chandrasekhar (1948) on different bases, but both take account of the 
polarization of the scattered light. Hammad modified an earlier step- 
by-step treatment (Hammad and Chapman 1939), considering only primary 
and secondary scattering, but his adopted law of scattering allowed for 
the anisotropy of the scattering molecules. 

By Chandrasekhar’s method, an equation of transfer is formulated. 
The solution of the associated integral equations yields an expression for 
the intensity of the scattered light, which includes contributions from all 
orders of scattering. Difficulties arise in this method when allowance 
is made for the anisotropy of the molecules (Chandrasekhar 1950 d). 


* Communicated by the Author. 
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Since such allowance is desirable (Hammad 1953), and since a compari- 
son of the contributions due to various orders of scattering is of interest 
in work on the vertical distribution of ozone (Chapman 1935, Dobson 1951), 
we here proceed in steps. However, we incorporate Chandrasekhar’s 
treatment of polarization, which facilitates the derivation of the secondary 
and tertiary expressions. 


§2. NoraTiIon 


Let M, and M signify the total mass of ozone and air respectively 
contained in a vertical column of unit cross-sectional area extending 
throughout the atmosphere ; m, and m the fraction of ozone and air 
respectively contained in this column above a height h. Then 


Momy=|_ po(h)dh and Mm= | p(h\dhs;os 0. Cea 


h 
where p,(h) and p(h) denote the density of ozone and air respectively at 
the height A. At the top of the atmosphere my=m=0; at the bottom 
My=m=1. The variable m will be called the level, and the vertical 
distribution of ozone may be regarded as characterized by a function 
Mo(m). 

The mass absorption coefficient of ozone and the mass scattering co- 
efficient of air, for light or wavelength A, will be denoted by «(A) and 
a(A) respectively, or, briefly by « and o. We also write 


a=aM, and s=oM. Nase SAR Fit 4 


Consider any point O in the atmosphere and construct rectangular 
Cartesian axes OX, OY, and OZ such that the positive direction of the 
Z-axis is vertically downwards, and such that the sun (regarded as a 
point) lies in the negative quadrant of the plane YOZ. Any direction 
may be specified by the polar angles (6, ¢) relative to OZ and the plane 
YOZ. We let dw, dw’, and dws respectively denote small solid angles 
about this direction, any other direction (6’, 4’), and the direction (65, 0) 
of the solar beam, 6) signifying the zenith distance of the sun. We shall 
also write 

p=cos 8, p’=cos 6’, and pp=Cos Oy. 

The intensity of monochromatic light which is received at height A and 
is travelling in the direction (6, 4), when the solar zenith distance is 05, 
will be denoted by R(h, 0, 4,4). The letters d and wu will follow the symbol 
R when we wish to distinguish between light travelling downwards and 
upwards respectively. We let H(h, 0, 4, 0)) signify the rate of emission 
of monochromatic light-energy per unit mass of air at height h, per unit 
solid angle about the direction (0, 6), when the solar zenith distance is 4). 
The parts of these functions referring to direct sunlight, primary, second- 
ary, and tertiary scattered light will be distinguished by the addition of 
the suffixes 0, 1, 2, and 3 respectively. The variables h, 0, and 4, specified 


in these functions will often be replaced by the corresponding variables 
m, p, and po. 
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§ 3. A REPRESENTATION OF PaRTIALLY PLANE:PoLARIZED Licut 


It may be shown (Chandrasekhar 1946) that the intensity of partially 
plane-polarized light (which is the most general type of light we consider 
in this problem) may be represented by three parameters, equivalent to 
three of the parameters introduced by Stokes (1852). For instance, if R, 
and £, denote the components of the light-intensity in the direction J and 
in the direction r at right angles to it, the intensity R(is) resolved in a 
direction inclined at an angle y to the direction of / is given by 


R(})= Rk, cos? +R, sin? 4+R, cosysing, . . . (3.1) 
where #,, denotes the third Stokes parameter. In the particular case 
of unpolarized (natural) light, R(x) must be independent of %s, whence 

= fea fe, and tt. ==() arecleee ected.) 
where F# denotes the total intensity. 

In the following work the intensity of received light is represented by a 
vector R, or {#,, R,, R,,}. the symbol R being reserved for the total intensity 
(2,+ #,). Similarly, the rate of emission of energy, per unit mass of air 
per unit solid angle, is represented by a vector E, or {H,, H,, H,,}, and # 
denotes (H,+ E,). 


§4. Tae ATTENUATION OF DiRECT SUNLIGHT 


The intensity of the light received at any height and in any direction 
will depend on the amount of attenuation, due to absorption by mole- 
cules of ozone and scattering by molecules of air, which the light has 


undergone. This dependence may be expressed by the relation 
hr 


R(h, 0, 4, 0) =R(h’, 6, #, 69) exp [—| (avec pines 6 dren emn rel) 


h 
In the case of direct sunlight, which we assume is unpolarized, the 
relations (3.2) are satisfied. By (3.2) and (4.1) we have 


et SS) ed 


0 
where Rd,(/, 95, 0, 09) represents the intensity of direct sunlight received 
at height h and travelling downwards in the directions confined to the 
solid angle about the direction (4), 0); and J,, briefly denotes Rd,(00, 6, 0,8) 
and is the total solar monochromatic energy received outside the atmo- 
sphere on a unit area placed normal to this direction per unit time. 
From (2.1), (2.2), and (4.2) we get 


Rd, (h, 8, 9. Ber ie exp[— | (apptop) secOdh, . (4.2) 
Yh 


if 
2 


Rdy(m, Mo; 9, eA 1 17, exp[—(amp+sm)/u], . + (4.3) 
0 | 


where Rd,(m, j49, 0, fo) replaces Rdo(h, 4%, 9, 85) (see § 2). Also | 
Rd,(m, 14, 9, 49) =0 when u-~Lo, and Ru,(m, u, &; Hy)=V- - (4.4) 
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§5. THe Law or ScaTTERING By ANISOTROPIC MOLECULES 


The Rayleigh law of scattering, generalized to include the case of 
scattering of arbitrarily polarized light characterized by a set of Stokes 
parameters (Chandrasekhar 1950 a) may be expressed as follows : 


E(m, , 6, bo) =(/477) P(u, d, pw’, 6’) R (m, pw’, 6’, ao) dw’. . (5.1) 


Here R(m, 2’, 4’, uo) represents the intensity of the light received at level 
m and travelling in the directions confined to the solid angle dw’ about 
the direction (cos~! yu’, 4’); E(m, uw, J, 49) represents the rate of emission 
of scattered energy, at the same point per unit mass, per unit solid angle 
about the direction (cos-!yu, 4); and P(u, 4, u’, d’) denotes the phase- 
matrix for Rayleigh scattering. 

The explicit form of P(p, d, »’. 6’) has been derived by Chandrasekhar 
(1950 b), who has also shown (1950 c) that when the anisotropy of the 
scattering molecules is included, P(p, 4, 1’, 6’) in the case of partially 
plane-polarized light should be replaced by 


: 12180 
(Te) rc (u, d, nw’, 6’) TES ray] 128 Cab iso) 
000 


where y=6/(2—6), and 8 is the depolarization factor. 
In the present work it is convenient to denote (5.2) by 


3 Q(t, db, bK’, ¢’)/(2+8). 


The law of scattering by anisotropic molecules then assumes the following 


form : 
30 
Elm, Bs $s Ho) = S(T gs) Me bo M's PYR(m, w’, P' wo)dw’. (5.3) 


The explicit form of Q(u, 4, u’, 6’) may be written as follows : 
Qu, $1’, 6')= QO (1, pw’) + (1—3)(1—p2) 12 (1—p’2)8? @M(, 6, pW’, 6’) 
' +(1—8) Q® (un do's!) = 7. ee ees 
where 
(L—p?—p!2+ 3 2"?) du oi E: 1 ae 
Q (nu, wu’) =(1—8) Au’? t+ O; +406); 1 1 Of; 

0 0 aa | | 
shay cos & —¢) 0 wsin ashi 
Q™Y (u,b, bw’, b= 0 0 
| aur ain ig’ 2p) oe acer 


and Q® (un, 4, nu’, 6’) 
Hn cos 2(f’—d) —du2 cos 2(¢’—4)  4y2u' sin 2 e4 
= Shae te (¢'—¢) + Cos 2(4’—¢) “in (5.5) 
—pp” sin 2(4'—d) je sin 2(4’—d) bye’ COS 2 aa 
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$6. THE Emission anp Recuprion or Primary Scarrerep Lieut 


The rate of emission of primary scattered light, per unit mass of air at 
level m, per unit solid angle about the direction (cos~tu, 4) will be 
represented by E,(m, u, 4, 49). It follows from (5.3) that 

Ei Gy) O(ubs us 0) Ryl mr, fig, 0. fh) d 
Oi aes eo) 87 (1+ 48) > P, Ko VY 0 > o> ), Mo) Wo- (6.1) 


Substituting for Ro(m, 9, 0, uo) from (4.3) and (4.4), we have 


oo. > 
E,(m, 1, b; bo) = 16x (1-45) Al, d; fo) K(m, po), - - (6.2) 
where 1 | 
Alu, d, Po) = Q(u, d, Ho> 0) ee igtte sere (6.3) 
0 
and* K(m, »)=exp|—(amy+sm))/9]. teeta ett 0-4) 


Putting ¢’=0 and substituting pu, for uw’ in (5.4), and carrying out the 
matrix-multiplication in (6.3), we get 


Ali, p, Ho) = AM (He, Mo) + 2(1—8)(1 —p?) (1 pg) "? rg A (12,6) 


+3(1—8)(1—py)A (x, $), poh ere a Cs) ACN) 
where 2— 1? — 2g? +32 U9? 1 ] 
A (tu, Ho) = 3(1—0) 1+ p49" eo |i, 
0 0 | 
[L COS & —p? cos 26 
AY(u, d)= 0 ,and A®)(u, d)= + cos 26 |. (6.6) 
sin —p sin 26 


From (6.2) we derive the following expression, which refers to the total 
intensity (see § 3) : 
Blom, ts d= TET Dy Alls Ba) Kons a)» «= (6-7) 
where, by adding the first and second rows of A(, 4, 49), we find that 
A (jy $5 49) =$(1 —8)[3—p2 pug? + 37 tu9? + 4pepey (Lp?) 1?(L p19) Cos 
L (1—p2)(1—pp?) cos 2f)+- 28. © Suess ae OrO) 
The reception of primary scattered light is treated by the method of 
Hammad and Chapman (1939), allowance being now made for the ab- 
sorption by ozone. Corresponding to eqn. (1) of § 12 of the paper citedt 
we have, in our notation, 


R,(m, 1, $, Ho) 
= z | E, (m’, 1, ; Ho) eXP[— (amy sm am, —sm')iu| dm, (6.9) 
pL 


* For the sake of brevity, the variables 4, s, and m, will not be specified in 

the function K, nor in several functions which are introduced later. 

The author's usage of m differs from that of Hammad and Chapman, who 
denoted this quantity by ™. 
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where 1,/ denotes the value of m, at level m’, and the limits of integration 
depend on whether the light is travelling downwards or upwards. We 
find that, for light travelling downwards (0<p <1) 


Rd, (m, wu, d; to) 
M Fi / / / 4 
=—| E,(m’, 1, 6, Wp) exp [—(amy+-sm—aim, —sm')/u| dm’, (6.10) 
MJo 


and, for light travelling upwards (—1<p <0), 
M i / / / 
Rw4(7, Lh, Ps 4g) = a E, (m’,u, 6,49) exp [—(amy)+sm—amy —sm )/w|dm 
(6.11) 


Substituting for E, (a’. u, ¢, 49) in accordance with (6.2), we have 
osie 


Rd, (m, 1, b; Ho) = Tén(14-45) A(1, b, Ho) Kd (m, 4, Ho), - (6.12) 
and 
—3slI, a 
Ruy (m, 1, bs bo) = Teo Uy A (1, $, fy) Ku (m, fb, Mo), - (6.13) 
where Kd(m, p, bo) =K (m, b) K (m’,—p) K (m', po)dm’, =. (6.14) 
FB; Bo 4 
1 
and Kulm, ph, bo) =K (m, »)| K(m’,—p) K (m’, 9) dm’... (6.15) 
0 


~ m 


Expressions for the total intensities Rd, and Ru, may be derived from 
(6.12) and (6.13). These expressions contain (6.8), and will not be given 
here. 


§7. THE Emission AND RECEPTION OF SECONDARY SCATTERED LIGHT 


The rate of emission of secondary scattered light, per unit mass, at 
level m, per unit solid angle about the direction (cos~1u, 4) will be repre- 
sented by E, (m, u, 4,49). It follows from (5.3) that the contribution to E, 
from primary scattered light received in the solid angle dw’ about the 
direction (cos~1y’, 4’) is 


30 
8(1+48) Q (u,b, 1’, f') Rd, (m, pw’, $'; Mo) dw’, when 0<p’ <1, (7.1) 
30 2 yt Faeys ’ , 
and 8m (1-+-40) Q(u, 4, x) )Ruy(m, po x) » bg) dw A when —l<p <0. (7.2) 


Forming the sum of the contributions from all directions, and replacing 
pe’ in (7.2) by —p’, we get 
30 a fe / , ‘ ‘ 
E, (m, Bb, d, Mo)= 8m (1-+45) | |, {Q (tu, p; > p )Rd,(m, bs d » Mg) 


0 


+ Q(x, 6,—pB', 6)Ruy(m,—p', $', Mo) fap’ dp’. . (7.3) 
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Substituting for Rd,(m, u’, 6’, wo) and Ru,(m,—p’, $', uo) in accordance 
with (6.12) and (6.13) gives 
9681, an rd ae ere . 
E(m.19; bo) = [De 45) | [ | eu $: MP )A(u’, $', uy) Kd (m, 11’, U9) 
+ Q(n,6,—w Pb) A(—w'. 6’, M9) Kulm, —p', Mo) $46! de’. 
(7.4) 
Our next step is to carry out the integration with respect to 4’. To 
facilitate this step, we express Q(y, ¢, 1’, 6’) and A(w’, ¢’, zo) in the forms 
(5.4) and (6.5) respectively, and we note the values of the following in- 
tegrals, in which p and q are integers : 


Se 0 if pA~q 
| cos p(¢’—¢). cos gd’. d6’= { 7m cos pd if p=q4~0 
0 


an (0 if pq 

ald’ =d). si br == 
I. cos p(¢'—¢). sin q¢'. dd lasinpd if p=g, 
an . , y} Va f 0 if pq 
I. sin p(}’—¢). cos gp’. do’ = | —7 sin pd if p=q, 
(0 if p~q 


Qn 
| sin p(¢’—¢). sin qd’. d’= ee cos pd ifp=q4~0. . (7.5) 


Pa 


We find that non-zero contributions only arise from the integration of the 
direct products Q° A®, QAM, and Q?) A®. After some reduction 


we get 
20 3 oa 
-| Ql, bw! BVA (HE Bs Ho) A= 2 BiH boo) WO. (7) 
TJ oO = 


Neglecting terms in 8”, we find that B;(, 4,49) may be expressed as follows: 
B,( 11, 5 to) = B; © (145449) +2 (1—28) (1p?) 1? (Lp?) 7 Bs (tes) 


at (1—28)(1—py”) B;) (Lu, dé), (@=1, 2, 3), 3 (7.7) 
where 
aie Baa | 7— 2p? — pg? 
0 1 1 2 1 45 5— py" 
By (u, bo) =2 (1—28) +Ho 2 0 
0 
— 34-4? + 5g? — 8" po” pee ae 
By! (u, po) =(1—28) | 1-H" +38{ 11 3y0° 
Lo 
2—3y? | 
B, (pu, Jig) = $(1 —28)(1 — 3p4" ) os | 
0 


BL) (u, d= B, (u, 6)=—3Bs Du, b= Nae b). 


and B,) (1, 6) =$Ba™(125 6) = Bs, =A (Usd); [see (6.8)]- (7.8) 


288 G. F. Walton on Atmospheric Scattering and 


Changing the sign of y’ in (7.6), we also have 


1 2% F 3 cite he 
=| Os Hs BAH $s Ho) Ab’ = E BiH, $, wow’. (7.9) 
0 ce 
From (7.4), (7.6) and (7.9) we then have 
9osT,, ~ L 710 
E.(m, bs $, bo) = 1287(1-+ 48)? me B. (1, b, bo) L a1 (™, Mo) }: - (7.10) 


where 
v5 . , i 
Ly; 4(m, 1o)=| [Kd (m, pw’, 49) + Ku(m,—p’, Wo) ]u' du’, (7.11) 
0 


the functions Kd and Kw being given by (6.14) and (6.15) respectively. 
The following expression, referring to the total intensity (see § 3), may 
be derived from (7-10) : 


Ey (m, bw, b, bo) = a Ed y B; La;-1(m, bo), - (7-12) 
1287 (1+ 48)? j21 
where B; briefly denotes B(y, 4, 49). We find that 
By = 3(1—28)[5— 3p? — 3 p19? + 5p? pg? + 4ppig (1p?) 1? (1 —pg”) "cos 
+ 3(1—p?)(1 pup”) cos 2¢]-+8 (6—p? po”), 
By=(1—28)[— 2+ 4p? + 419? — 6719? + 2pupey (1 —) 1? (1— py”) ¥? cos 6 
+3(1—p?)(1— p19?) cos 24 }+-8(—2+3p2+ 3u9”), 


B3=}(1—28) [ (1—3p?)(1 39?) — 8g (1 —p2?) 8? (1 pup?) 1? cos 
+43(1—p?)(1—j15?) cos 24). st ee SO cee ee 
The treatment of the reception of secondary scattered light follows that 
outlined in § 6 for the reception of primary scattered light. Thus we have 


9s? J, 3 
Rd» (m, Bb, d; fio) = 12877( l Ate 18 )2u me B; (L, d; Lo) Ldy;_4 (m, BK, Ho)}: 
2 i= 
(7.14) 
and 
SCY ee 3 
Rv, (m, LL, dp, Lo) = | 287(1 os 18) = { B, (K, Ps [o) Lug; 4(m, BK; Ho) 
2 v= 
\ (7.15) 
where 
P 7m 
Liddy; 4(m, 1, wy) =K (m, 1)| K (m', —p) Lo 1 (m'; po) dm’, (7.16) 
/ 0 
and 
“1 


Lg; (Mm, b, fo) =K (m, p) | K (m', —p) Lg; 4(m’, ty) dm’. (7,279 


m 


*The author has verified that the corresponding expressions of Hammad 
(1948), which refer to the directions normal and parallel to the plane of scatter- 
ing, combine to give (7.13), apart from two discrepancies, which have been 
traced to small errors in Hammad’s expressions for B, and C4. 
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§8. THe Emisston anp Reception oF Terttary Scarrerep LIGHT 


The treatment of tertiary scattering follows that outlined in the 
preceding sections, the number of the azimuthal expressions and the 


structure of the multiple integrals being greater in this case. Corresponding 
to (7.3), we have 


3c 27 
E.(m, Bb, d; Lo= ae i {Q(u, ps kK’, ¢') Rd.(m, [oa ¢', Ho) 
+ OH$, —p', b') Rus ( gah aeato) sles ae (8.1) 


Substituting for Rd,(m, »’, 6’, uw) and elie —', $’, Wo) i accordance 
with (7.14) and (7.15) we get 


Qos*I,, 3, (241 was Heed 
E3(m, p, ?: Ho)= Top4A(1-L IS), ls J [Q(u, dw’, $')B (u's $', Ho) 
X Lda, (m, M's bo) + Q(u;6,—1'. 6) B (1, P's Ho) Litas a (m, b's Mo) | 
Pen nace wreteal Me eS 48.3) 


We effect the integration with respect to ¢’ by the method of §7 
- making use of the results (7.5). In this case we write 


1 2% ‘ 3 Ra : 

a Q(x, db, ie ?') Bt’, P's o)dd = 2 C.(u, Ps [éo) Ne ) @=1, 2,3), (8.3) 
GTR EO) jail 

the same result being obtained when the sign of y’ is changed. Neglecting 
terms in 6?, we find that C,(u, 4, 4») may be expressed as follows : 


C(u, bs My) = Sj (4, Mo) +2 (1 — 38) (1 —.2)2(T— pg?) 8? fg Cys" (He, ) 
+$(1— 38) (1—py?) )C;; red (Lu, ), ( ), (@=1, 2,3; 7=1, 2,3), (8.4) 


where 
Cc, (LL; Ho) 
8— Tp? — 89? ia | 19 — 9? —5ypg? "p29" 
= ey sia +48} 10—4y,? 
0 (Shae 
C12 (u, Ho) i 
T+ 9? + Mr? — L179" 21 + 25? + Tpg"— 3% ULg 
=(1—38)} 2 2p" +45 4} 4g? 
0 0 
6— 9p? — 10g? + 15 wg” 2— 32 
C3 (hu, (9) =3 (1—38) —3+ Buy" +6 at ; 
) ) 
0 


— 6-4 Tp? 10pig?— pig? 
¢C,, (Le [ég)=(1—38) 1—pL)" | 
0 a 
— 144 11p2+ 14y,2— 376? |] 
oo = eee 
0 
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14—17 p?— 22 wy? + 27? 1,” 


Coo! (14, fg) =(1—38) |. —3-+-5y19? 
0 
30— 27 p2—18 pip? +9 p22 py? 
+46 3—9uy” > 
0 
— 84122 + 12,2182 149? —6+ 9p? 
Cas(4,449)—=(1—38) | 4— 6g? oe (lags 
0 0 
4*~ 5 2 | l 
C51 (wo) = 3 (1—35)(1—3 p19?) oy +35(L—3p9") | 1], 
; | 
—5+6p2 —3 
Coo") (449) = (1—35)(1—3 yup?) | +28 (1—3p9") | —3 J, 
0 0 
6— 9p? 
C53 (4, Wy) = 4 (1 38)(1—3149”) =e ? 


0 
= = = = ease = 
C31) (u,b) =Cyo ¥ Cys = CC.) V9=C, Y= —} 23°90 = —3 Cy) 
= — $3 Cg = 3 C55'V =A (1,4), 


and 


Cai?) (us, ’)=3 


@) 
an 
te 
I 
(@) 
a 
3S 
I 
bol 
i) 
bo 
~ 
z 
I 
Hale 
fe) 
no 
rn 
re 
| 


4C,,(2)— C,,2)—1 C,,) 
= C3, = A® (u, 4), nit el ee ee 


the forms of AM (yu, 6) and A® (u, 6) being given by (6.6). 
From (8.2) and (8.3) we have 


27 o8*1 Le 
E,(m, 1, b; Mo) = 024m (1-45), Father its $: bo) Mo; 1, 23-1(™s Ho) }, (8.6) 
where 


1 ; 
Mai -1, 95-1 (I, Hy) = i [Ldy,_1(m,p', hy) +-Lug,_4(m, —p’, Mo) | ei )dy’ 
(8.7) 
The following expression, referring to the total intensity, may be 
derived from (8.6) 


2Tos*I, 3 3 


E(m, fb, $, uo) = 1024n(1--45)8-, ‘a Cj Masa a5-1(™, My),  - (8.8) 
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where C’;; briefly denotes C,;(u,¢,449). We find that 
© C11 = 3 (1—38)[9— Tp? — Tyg? + 9p? p19? + 4 prpug( L — 2)? (1 —peg2) 1? cos 
+4(1—p?)(1—p1g) cos 26]-+45(29—9p.?— 9192+ 2152), 
= (1—36)[— 5+ 9p? + Tg? — 11 p29? + 24 9(1 —?) 1? (1 —e92) 1? cos 
| Hi p2)(1—pog2) cos 26] +-48(—17-++25p.2-+ 11 ig? —3y2 p92), 
= }(1—38)[ (1—32)(3—B jig?) — 82pzq (1—p2)¥2 (1 pug?) co8 ¢ 
Sai )(1—p1g2) cos 26] +3 (1-342), 
Egy (1L—38)[— 5+ Tp? + Opp? — 11 pr? pg? + 2p fly (1 —?) 1? (1—pug?) 4? cos 
+$(1—p2)(1—pig2) cos 26] +43 (—17 + 11p2-+25 4g? —3 p24162), 
O'g9=(1—38)[11—1 7? — 17 pg? + 27 7 9? + 2pupig( 1 —p?)¥?(1 p19)? cos p 
+4(1—p2)(1—pig?) 008 26 ]-+45 (38—27 2 —27 ug? +- 90292), 
C yg (1—38)[—2 (1—3p2)(2—3 ig?) —4pupty (1 —p2) #2 (1 pug?) 12 0s 
+4(1—p?)(1—pig?) cos 24 ]—38 (1—3p?), 
C5, =4(1—38) (3—5p.2) (1—3 jg?) — 8 prq(1 —p?) 2 (1 tug?) "2? cos 
+-4(1—p2)(1—pu9?) 008 26 +8 (1—3 442), 
C'39=(1—38)[— 2 (2—3p?)(1 — 3p?) — 4 pty (1 —p:?) 1? (1 py?) 1? cos 
+-H(1—p2)(1—pq2) 008 26 J—35 (1—3 p42), 
C'gg=}(1—38)[3 (1—32)(1— 3 4g2)-+ 16 1p2g (12) 2 (1p?) #2 C08 4 
1 1(1—?)(1—py?) cos 2¢ J. ee res eee a 
The treatment of the reception of tertiary scattered light follows that 


outlined in § 6 and § 7 for the reception of primary and secondary scattered 
light respectively. Thus we write 


Rds (m, 1, b, Ho) 


= Wael Ts x >» {Cc oo) ) Md )} (8 10) 
(o> Ps Ho) Mex 1, aj-1 (Ms Ms Mo) f> ; 
10247(1+ 46)? Tea See) (He 0 27— 15. 25-1 0 


and 


Rus (m, 1; b; Ho) 
1S oh 
~ 10247 (1+ 48) 8x ; 


where 


3 3 ; . 
EZ {Cy i(w, b, Ho) Mo, 1.95-1 (Ms Ms Mo) fs (8.11) 
=1 j=1 


Md; 1, 2j-1 (m, /; Mo) 
cin 
= K(m, p) | K(m’, —p) Moj-1 95-1 (1's Mo) dm’, 


~ 0 


—~ 
mn : 
= 
i 
— 


and 
Muy, 1, 25- ee Ls Ho) 
=K(m, w)| K pe reat) Moj-1, 25-1 (Ms bo) EM’. 


m 
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§ 9. CONCLUSION 


The preceding analysis completes the formal solution of the problem 
considered. We now have expressions for the primary, secondary, and 
tertiary scattered light-intensities emitted and received in any direction 
and at any point in the atmosphere. For each order of scattering, the 
vector representing the light-intensity is given in terms of the product 
of an azimuthal vector and a function generally involving a multiple 
integral. Similar products would evidently be derived by extending the 
treatment to higher orders of scattering. 

We have had in mind the scattering of ultra-violet sunlight, but the 
particular results which would obtain by reducing the absorption co- 
efficient of ozone to zero, might be utilized in the study of the ilumination 
of the earth’s atmosphere by visible light. Numerical results for an 
atmosphere containing an absorbing layer of ozone will be given in a 
forthcoming paper. 
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THE existence of the group of contact transformations, or that of the 
continuous Schrédinger group whose elements are commutative with 
Hamiltonian is essentially related not only to the prolem of degeneracy, 
as shown by Jauch and Hill (1940), but also to the problem of quantization. 
In the present paper it is shown that on the basis of the representations of 
the semi-simple Lie group generated by the integral operators F, 
(A=1, 2,...,7), which can be made by the method of classical analogy 
from a set of y independent integrals of the equations of motion in 
Newtonian mechanics, an eigenvalue problem can be constructed in 
terms of a complete set of the generalized Casimir operators 


CRS ER iG) ea ra le emer awl So io, ke Pee a mel) 
where / means the rank of the semi-simple group, and all G, constructed 


from the invariants of the adjoint group are commutative with any /;. 
An explicit formula of G;, as pointed out by Racah (1950, 1951), is 


DR ee RE Te oh See gO Mae earl eae 2) 
Ansan an As Pccoialy v 
and metric fundamental tensors are given by 
Gears Devi Crenice Ges fone hana arn) 


ial oth o 
The structure constants cj, of the Lie group may be at most operators 
depending on the Hamiltonian alone, and in principle they are not always 
identical with classical ones. The quantization procedure by the use of 
the representations of these operators is generally applicable to any 
quantum-mechanical system, as far as it is possible to find the above 
mentioned Schrédinger group for the system, although this restriction 
is considerably rigorous. 

Let us consider the present problem about a typical example of the 
hydrogen like atom. The six integral operators commutative with the 
non-relativistic Hamiltonian 


Heap 2m Let ene, «= (4) 
of this system are, as well-known, the Lenz—Pauli vector integral 
A=(1/2mZe2)(Mxp—pXM)+r/r . . ws (5) 


in addition to the orbital angular momentum operator M=r xp, where 


* Communicated by the Author. 
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charge, respectively ; 7, the distance between a given nucleus and electron, 
Ze, the nuclear charge. For the negative energy state the Balmer formula 
was formerly derived by Hulthén (1933) using the commutation relation 
of the above integral (5) and it was already shown by Fock (1936) and 
Bargmann (1936) that the Schrédinger group of this system is isomorphic 
with the four-dimensional rotation group, while from our point of view 
the set composed of the above six operators might be regarded as a 
semi-simple group of the rank 2, therefore there must exist a set of two 
generalized Casimir operators. But in this case it can be found that only 
one ordinary Casimir operator (1931) of these has the following functional 
relation with the Hamiltonian : 
+ 
G>=—}42 Ff=—(Zetm/27¥7H)—1, . . .-. (6) 


= 


A=1 
where # is the Planck constant h divided by 27. Each highest weight 
characterizing the irreducible representations of the Lie group is (0, 0), 
(1, 0),..., (k, 0),..., and then the degrees of their representations are 
1?, 2?,..., (k+1)?,..., respectively. The eigenvalue of the Casimir 
operator determined by a given highest weight (k, 0) is k(k-+-2), which is 
also connected with the eigenvalue of angular momentum in the four- 
dimensional Euclid space. Consequently by substituting » for k+1 the 
relation (6) can be reduced to 


n*?— 1 = — (Z*e4m/2h7H,,')—1, 
or 
Ae — Zretmi2hn* > |\(n=12, 9.5); cas See 
which is the well-known Balmer formula accompanied with the degree 
of degeneracy (k-+-1)?=n? as mentioned above. 

On the other hand it is known that the Schrédinger group of the positive 
energy state of this system is isomorphic with the Lorentz group. In the 
case of the continuous eigenvalue we (1953) found an infinite continuous 
(unitary) representation of a non-compact Lie group, which will be 
presented elsewhere, besides some applications of this procedure to other 
examples and the perturbation theory had also been formulated. 
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ABSTRACT 


The difference between the total cross section of deuterium and 
hydrogen has been obtained at effective energies of 109, 117, 132, 149 
and 169 Mev from attenuation measurements on light and heavy water. 
Using published values of the hydrogen cross section, the oxygen total 
cross section was also obtained at these energies. 


§ 1. [YTRODUCTION 


PREVIOUS measurements of the difference between the neutron total 
cross sections of deuterium and hydrogen have been made in the high 
energy region at 42 Mev (Hildebrand and Leith 1950), 85 mev (Cook, 
McMillan, Peterson and Sewell 1949), 95 mev (DeJuren and Knable 
1950), 270 Mev (DeJuren 1950), 280 mev (Fox, Leith, Wouters and 
MacKenzie 1950) and 400 Mev (Nedzel 1954). Neutrons produced by 
the 230 cm Uppsala synchro-cyclotron were used in the present work, and 
the experimental method was similar to that used in determinations of 
neutron total cross sections at Harwell, England (Taylor and Wood 1953). 
The results at 169 Mev have been communicated earlier (Taylor 1953). 


§ 2. EXPERIMENTAL 
(a) Geometrical Layout 

The circulating proton beam in the cyclotron was intercepted by a 
beryllium target (30 mm thick), the radial position of which was changed 
in order to obtain the different neutron energies. The neutrons emerging 
from the target passed successively through a thin aluminium window 
in the tank wall, a steel collimator, the attenuator, another steel collimator 
situated in the concrete shielding wall, and then struck the polythene 
The target, collimator holes, attenuator holder, and 
polythene scatterer were carefully aligned. The distance between the 
target and polythene scatterer was 14 m, between target and attenuator 
8-5 m, and the diameter of the steel collimator holes was 39 mm. 


scatterer (fig. 1). 


* Now at A.E.B.E., Harwell, Berks., England. ; 
+ Communicated by Professor The Svedberg, F.R.5. 
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(b) Neutron Detector 

Recoil protons emerging from the polythene scatterer (17 mm thick 
and 50 mm diameter) at an angle of 20° to the direction of the neutron 
beam were detected by a triple coincidence counter telescope using 
scintillation counters. The polythene scatterer was perpendicular to 
the direction of the telescope. Each scintillation counter consisted of a 
polystyrene crystal (7 mm thick and 25 mm square) mounted on top of a 
photomultiplier tube (EMI 5311) in a light tight aluminium case. An 
aluminium absorber was placed between the last two counters in order 
to set a lower energy limit to the recoil protons detected (see table). 
The coincidence circuit was of the bridge type using germanium diodes, 
and had a resolving time of 10-* sec. A block diagram of the electronic 
equipment is shown in fig. 2. 

Fig. 1 


NEUTRON EXPERIMENT 
OLOME TRA LAYOUT 


Zz - NEUTRON BEAM 
Z-< 
Z| a ee 


COLIMATOR sae ALUMINIUM 


CYCLOTRON 


In order to find the best operating conditions of the coincidence circuit 
the counting rate was recorded for various counter and discriminator 
voltages until a plateau of counting rate against voltage was obtained. 


(c) Neutron Monitor 


The neutron output of the cyclotron was monitored by a BF, counter 
placed in a cylinder of paraffin wax, situated behind a concrete shielding 
in the cyclotron hall. The operating conditions were determined as for 
the counter telescope. 

(d) Attenuators 

The ordinary and heavy water samples were contained in copper tubes 
(52 mm diameter) with copper end windows (0-5 mm thick). A remotely 
controlled sample changer held the required sample in position, and also 
had a position for an unattenuated beam. The tubes of ordinary and 
heavy water were of such lengths that they contained the same number of 
oxygen atoms when filled. The heavy water contained 99-7 5% of D,O. 
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Effective neutron 

energy Mev 169+3 | 149 
Energy band width 

Mev 19 
Contribution to the 

counting rate from 

carbon content of 

scatterer % 10 12 12 12 
Aluminium g cm~? 12-97 9-19 6-48 3:78 
Deuterium minus hy- 

drogen cross sec- 

tion 10-2” em? 23-1+.2-0 | 24-8 +2-0 | 25-3 42-4 | 27-2+42-4 | 29-1436 
Oxygen cross section 

10-77 cm? 431+6 | 47346 | 518+6 | 598+6 | 656+6 
Hydrogen cross sec- 

tion used in deriving 

oxygen cross sec- 

tion 10-?? cm? 48 +2 50+ 5442 59+2 64+ 


Fig. 2 


jo} 


NEUTRON EXPERIMENT” 
BLOCK BIAGRAM 


TRIPLE 
SCALER 
1009 A 


COUNTER 
TELESCOPE 


SCINTRL ATION 
COUNTER aa 
SCALER 
10098 


kength ofcobtes from 
Cyctotron Holl to 
Control Room LOmerres 


WEUTRON 
MONITOR 


§ 3. Errective Neutron ENERGY 


The counting rates were recorded for various thicknesses of absorber 
in the telescope and these thicknesses were converted to neutron enaneies 
using the range-energy relation and the Pokniaiseies ae sconeeonng 
recoil proton energy to incident neutron energy. The able aoe 
energy was taken as the mean value of the neutron Broigies sae 
by the corresponding number of recoil protons. This methud assumes 


x2 
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that the total cross sction does not vary appreciably over the energy 
band in question. Further, if an energy band greater than 13 Mev is 
used there will be a contribution to the detected protons from the carbon 
content of the polythene and this will have a different effective neutron 
energy. To estimate this effect the counting rate was measured for a 
carbon scatterer of the same stopping power as the polythene scatterer. 
This contribution from the carbon content in the polythene scatterer, 
together with the energy band used, are recorded in the table. Since the 
radiator and the last crystal have a finite thickness which cannot be 
neglected and since the solid angle is not negligibly small, there is an 
uncertainty in the calculation of the effective neutron energy. An 
estimate of this uncertainty is included in the effective neutron energy. 


§ 4. EXPERIMENTAL PROCEDURE 


A total counting time of about 40 hours was required to obtain an 
accuracy of 10% in the deuterium minus hydrogen cross section. In 
order to reduce possible systematic errors due to long-term shifts in 
electronic equipment the experiment was divided into shorter runs of 
about one hour. A complete run consisted of successive measurements of 
the ratio of triple coincidence counts to monitor counts with each of the 
two attenuators in position and without an attenuator. The overall 
spread of the readings was no greater than expected from the total counts 
taken. 

By inserting delay cables in the third channel of the triple coincidence 
counter it was found that the background of random triple coincidences 
was negligible. 


§ 5. RESULTS 


The cross sections calculated from the measured attenuations are listed 
in the table. The light water content in the heavy water sample was 
found to be negligible when calculating the deuterium minus hydrogen 
cross sections, but small corrections had to be applied for the copper 
end windows when calculating the oxygen cross section from the attenu- 
ation in light water. The values of the hydrogen cross section used in 
the latter calculation were obtained by interpolation from published 
values and are listed in the table. 

The values obtained in this experiment agree with those obtained by 
others, at higher and lower energies. 
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ABSTRACT 


A technique for the direct counter detection of heavy unstable cosmic 
ray particles stopping in a cloud chamber is described. It is shown that 
the method is about three times more efficient than methods previously 
used when detecting particles with lifetimes greater than 5x 10~° sec. 
The observation of eight, well-identified, unstable particles is reported ; 
they all have properties consistent with the K,,-particle. 


§ 1. INTRODUCTION 


Since the initial discovery of charged heavy mesons cloud chamber 
observations of these particles have been obtained by photographing 
penetrating showers. Nuclear emulsion experiments (Fowler e¢ al. 1951) 
indicate that about 3 x 10-* heavy mesons stop per gram per day beneath 
30 cm of lead at an altitude of 3000 m. However, multiplate cloud cham- 
bers, when selecting energetic penetrating showers in a conventional 
manner, identify only 0-6 10-® S-particles per gram per day at this 
altitude (Bridge et al. 1953). It was therefore thought worthwhile to 
attempt the direct selection of charged heavy mesons. This paper 
describes an experiment in which a multiplate cloud chamber is triggered 
by a counter system which responds directly to heavy unstable charged 
particles. About 1-4 10~* heavy mesons are identified per gram per day 
in this experiment. 


§ 2. PRINCIPLE OF EXPERIMENT 


The selection system and cloud chamber are shown in fig. 1. The 
counter system above the cloud chamber consists of two liquid scintil- 
lation counters S, and 8,, a Cerenkov counter C,, and a double row of 
Geiger counters G. There is a layer of 1-5 cm of lead between S, and G. 
The cloud chamber beneath these counters contains six lead plates and 
a Cerenkov counter Cy. 

The Cerenkov counters give signals only when traversed by fast particles. 
They are lined with diffusely reflecting surfaces, and respond to the total 


amount of Cerenkov light emitted without making use of its directional 
properties. 


* Communicated by G. D. Rochester 
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The counter system above the chamber selects heavy particles by means 
of coincidences S,+8,+G—C,, which will be referred to as H coincidences. 
This selection implies that a particle with velocity less than a critical 
velocity v,, and range greater than R, g cm-2 (the amount of material 
between C, and G) has traversed the system. Thus the particle must have 
a mass greater than M,(v,,, Ry). This system of selecting heavy particles 
has been described previously by Duerden and Hyams (1952), who 
showed that M)c?>6-8 R, mev when »,, is the critical velocity for water. 
In this experiment R, is chosen so that M,> 300 electron masses (m,), 
in order to exclude light mesons (7 and yp). 


Fig. 1 


LEAD 


PM PHOTOMULTIPLIER 


Front view scale drawing of counter selection system, and cloud chamber. 


The cloud chamber is expanded when coincidences H+C, occur within 
a resolving time of 5x10-* sec. In general no single downward-moving 
particle can actuate the system, since it would have to be moving slowly 

idly at Co. 

i ees oectorn is Eaained to select slow particles with mass greater 
than 300 m, which emit fast decay or annihilation products. 4 
- The C, signal amplitude for H+C, coincidences is recorded on a cathode 
ray oscillograph to enable the velocity of heavy-particle decay products 
re | eee unstable particles it is necessary to place a layer 
of material above the chamber. However, a slow heavy meson will a 
actuate the H system if accompanied by a fast particle giving ap es i- 
coincidence signal in the C, counter. For this reason the ee ete 
of the production layer was placed 80 cm above the apeater ee 0 a ee 
the slow heavy mesons to separate from any fast shower particles asso 
ciated with their production. 
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In the first half of the experiment a production layer of 80 cm depth 
of paraffin was used. It was thought that this would give few fast shower 


particles accompanying the heavy mesons, and that y-rays from any — 


7-mesons produced would give only small electron showers in the 
paraffin. In the latter half of the experiment a layer of 14-5 cm of lead 
was used instead of the paraffin, in the belief that heavy mesons may be 
produced mainly by ~-meson nucleon interactions, in which case large 
nuclei should be better targets for heavy meson production than small 
ones. It was hoped that any electron cascades associated with the 
heavy meson production would be absorbed by the lead. The area of 
each production layer was 150 x 150 cm?. 

There was no significant difference between the rates of observing un- 
stable particles beneath lead and paraffin. Beneath paraffin the apparatus 
was triggered much more often by air showers and it was necessary to limit 
the counting rate by requiring coincidences from an additional set of 
Geiger counters placed above §,. 

The time required for a particle produced in the paraffin or lead to 
come to rest in the cloud chamber is about 5 10~° sec. Therefore the 
system can respond efficiently only to particles with lifetimes between 
5x10°® and 5x10-® sec. Antiparticles, stable in flight but with very 
rapid annihilation, can also be detected. 


§ 3. DESCRIPTION OF APPARATUS 


The scintillation counter §, is filled with terphenyl in xylene, and 8, 
with terphenyl in phenylcyclohexane. Each of these counters has an 
active area of about 20 x 20 cm?. 

The Cerenkov counter C,, filled with distilled water, contains two photo- 
multipliers to. give high efficiency. It has an active area of 22x 22 cm? 
and a depth of 6-6 cm. Under normal running conditions pulses greater 
than about one-tenth the mean signal from relativistic particles give an 
anticoincidence signal. 

The counter C, is made of } inch brass and has a sealed dry compartment 
at one side for the base and wiring of the photomultiplier. The active 
area of the counter is 21x27 cm? and the depth of the liquid is 6-3 cm. 
A 19-core cable brings out connections to each dynode of the photomulti- 
plier and also two small tubes to maintain the dry and the liquid-filled 
compartments at atmospheric pressure. The dynode resistor chain is 
mounted externally to avoid having a source of heat in the cloud chamber. 
The signal from the collector electrode is brought out through a separate 
coaxial cable. Both cables are sealed to the counter and the cloud chamber 
by means of rubber O-rings. In over six months of operation no difficulties 
have arisen from running this counter inside the cloud chamber. 

The counter ©, is filled with a mixture of equal volumes of glycerine 
and water in order to lower the critical velocity and thus increase the 
efficiency of detecting decay particles. The critical velocity, v,, for this 
mixture is 0-71 c, which is still sufficiently high so that a slow proton which 


en 6 many tg 


Apne petite 


Study of Charged Heavy Unstable Cosmic Ray Particles 303 


just fails to give a signal in C, has a very small probability of giving 
asignalin C,. Signals in C, greater than about one-quarter of the mean 
signal from fast vertical particles give coincidence pulses. 

All the photomultipliers are E.M.I. type 6262 14 stage tubes. 

The cloud chamber, which has a photographed region 35 cm high, 
28 cm wide, and 35 cm deep, is illuminated from the rear, the expansion 
being obtained with two rubber membranes at the sides of the chamber. 

Each of the six lead plates has a thickness of about 1 cm, with a weight 
of 11:0 g cm-?. The total vertical thickness of C, is equivalent, with 
regard to ionization losses, to about 20 g cm-? of lead. 


§ 4. OBSERVATIONS 


The observations described in this paper were made at the Observatory 
of the Pic du Midi, at an altitude of 2850m. The equipment was run 
for a period of about four and a half months during which it operated 
successfully for about 2000 hours. 

Under 14-5 cm of lead 2520 H coincidences were registered per day. 
From a run of three days, photographing such events without requiring 
a coincident C, pulse, it was concluded that of these events 529 per day 
were due to protons stopping in the part of the chamber which was 
photographed. 


Spurious Hvents 


The unstable particle selection system (H-+C, coincidences) registered 
67 counts per day. Since the chamber recycling time was 5 minutes, 
56 pictures per day were photographed. Analysis of these pictures 
indicated that 19 per day showed evidence of an electron cascade in the 
chamber. Most of these were probably events in which different particles 
from electron showers traversed the counters 8,, S., G, and C,, without 
passing through C,. 

Three photographs per day showed a slow proton traversing the 
photocathode or first dynode of the photomultiplier in the counter Cp). 
This rate exceeded the accidental rate of such events by a factor of at least 
12. These events can only be attributed to pulses in the photomultiplier 
initiated by a charged particle causing fluorescence or direct secondary 
electron emission in the tube, rather than any effect due to Cerenkov 
radiation. 

Ten pictures a day showed single particles traversing the cloud chamber, 
with trajectories passing through the selection system. Most of these 
may be attributed to fast protons which, due to statistical ea eons 
in the signals, failed to be detected at C, but were detected at Cy. , 

Twenty-one pictures a day showed no charged particles traversing os 
counter C,, in spite of the fact that its whole active area was pnosgrayned : 
These can be attributed to y-rays from diffuse air showers, and the rate 


consistent with this assumption. 
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Three pictures a day showed a particle coming to rest in the cloud 
chamber lead plates, with no visible charged secondary particle, and no 
visible charged particle traversing Cy. This rate was eight times greater 
than could be accounted for by random coincidences between stopping 
protons and random counts (100 per sec) in the counter C,. Some of 
these events were probably due to slow protons giving rise to nuclear 
excitation with the detection of subsequent nuclear y-rays by Cg. 


Unstable Particles 


About 0-1 pictures a day showed particles appearing from ionization 
to stop in one of the lead plates and having minimum ionization tracks 
which begin within 1 cm of the point of stopping and enter Cy. Plate 
4 shows two such events, which will be referred to as S-events. In 
all, 8 well-identified S-events have been photographed. In three of these 
events the range of the secondary particles is more than 60 g cm of 
lead. Of the 8 events 5 were observed under the paraffin production 
layer. 

None of the secondary particles show evidence of multiplication or 
large angles of scattering as might be expected of electrons, although 
four of the eight traverse at least 1 cm of lead before leaving the chamber, 
and one of the others traverses C,. The remaining three enter C, with 
steeply inclined trajectories which go out through the sides without 
reappearing in the visible region. 

The fact that few of these could be due to stopping light mesons is 
evident even without reference to the properties of the H-selection system 
for all the S-events had lifetimes less than 4 x 10-8 sec (see following paper). 
Since the resolving time of the counters is 5 x 10~* sec the probability is 
2x 10-* that one of these events is due to a 7- or y-meson, and 1-6 x 10-7 
that four of the events are due to 7- or p-mesons. 

Taking these observations into account we consider that it is unlikely 
that even one of the S-events could be attributed to a stopping z- or 
#-meson, and that they are all due to the stopping and decay of heavy 
mesons. 

In none of the S-pictures do we observe any evidence for y-rays 
associated with the stopping particles, either in the direction opposed 
to that of the charged secondary, or in any other direction. 

In principle, the track length of the secondary particle through the 
counter C,, and the C, pulse amplitude, determine the velocity of the 
particle as it traverses the counter. In practice, the fluctuation of the 
signals makes any one velocity determination very inaccurate. However, 
if it is assumed that all the secondary particles are monoenergetic j.-mesons 
and allowance is made for the slowing down before reaching C,, the mean 
velocity of the eight particles may be determined. It is found to be 
(0-93*$'03)¢ which corresponds to a mean initial momentum of (275273) 


Mev/c, the errors being calculated from the r.m.s. deviation of the mean 
C, pulse amplitude 
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§ 5. Discussion 


Mt is consistent with present knowledge to identify all our S-events 
with the K,-particle identified by Gregory et al. (1954). In no case did 
we bring the charged secondary particle to the end of its range, since the 
chamber was too small for this to be probable. Thus we may place only 
a lower limit on the range of the decay particles. To give a coincidence 
in C, the momentum of the decay particles, if assumed to be /e-mMesons, 
must exceed 110 mev/c when they traverse the counter. The mean 
initial velocity observed for the eight particles using the C, pulse 
amplitudes corresponded to an initial momentum of (275418) mev/c 
which is compatible with that of the K,-meson. The absence of y-ray 
decay products, and the long distance of flight before decay of these 
heavy mesons is also consistent with the properties of the K,-particle. 

For every S-particle photographed 1750 protons traversed the selection 
system and stopped in the corresponding plates of the cloud chamber. 
If allowance is made for the geometrical inefficiency for detecting the 
S-particle decay products (assumed to be p-mesons with momentum 
220 mev/c) we find a ratio of 300 protons to one S-particle stopping in the 
chamber. 

The inefficiency of counting S-particles due to the C, anticoincidence 
may be estimated in various ways. We may compare the observed S- 
particle counting rate with that predicted by the emulsion data. We may 
also compare the observed ratio of stopping S-particles to protons with that 
obtained by Gregory (1954) and note that we obtain the same proton 
flux as Miller e¢ al. (1954). From these arguments we conclude that our 
efficiency for S-particle detection is certainly not reduced to less than 
30% due to the C, anticoincidence. 

Stopping z-mesons can give rise to relativistic decay particles, either 
through the usual 7 decay, or through electron pairs from the probable 
alternative 7°~y decay. In the former case the detection efficiency for the 
electrons is not very high. In either case, positive identification from the 
cloud chamber pictures is difficult. Although some events compatible 
with z-mesons have been observed, no satisfactory conclusions can be 
drawn from these. 

The selection system would respond efficiently to stopping anti-protons, 
but although during the running time of the experiment 5 x 10* protons 
have traversed the counters and come to rest in the chamber, no evidence 


for such particles has been obtained. 
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ABSTRACT 


A new counter method for measuring the lifetime of unstable mesons 
is described. Lifetime measurements on 13 charged heavy mesons are 
reported. It is concluded that if there exists only one heavy meson giving 
fast decay products, with a lifetime greater than 5 10~° sec, its lifetime 
is (11-0 *$4)x 10-® sec. No evidence for any longer-lived heavy mesons 
has been found. 


§ 1. IyTRODUCTION 


For some years it has been thought that there exists at least one type of 
charged heavy meson with a lifetime longer than 210° sec. The 
measurement of such a lifetime (see Rochester and Butler (1953) for a 
survey of early work), has eluded previous cloud chamber observations. 
This is because the lifetime is long compared with the time of flight across 
a chamber, and in a multiplate chamber, as in nuclear emulsions, 
practically all the particles come to rest before decay. Using the particle 
selection system described in the previous paper (referred to as I), a 
scintillation and a Cerenkov counter enable the lifetime of single heavy 
mesons to be measured directly. The particles are identified from cloud 


chamber photographs. 


§ 2. MetHop oF MEASUREMENT 


The counter system described in I (see fig. 1 of that paper) was used in 
conjunction with a cathode ray oscillograph to measure the time elapsing 
between a heavy meson traversing counter S, and a decay particle 
traversing counter Co. 

The time scale is calibrated by measuring the time base velocity and 
the separation of pulses from fast cosmic ray particles which traverse 
counters S,, S5, and Cy. 

These observations, in conjunction with the cloud chamber photographs, 
enable one to calculate the lifetimes of the heavy mesons in their rest 
system. In fact the contributions to the lifetime due to times of flight 


were in all cases less than 5 x 10~1° sec. 


* Communicated by G. D. Rochester. 
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§ 3. TECHNIQUES 


A block diagram of the counters and fast electronic selection and timing 
system is shown in fig. 1. The auxiliary circuits used in conjunction with 
the cloud chamber to give a 5X10-® sec resolving time are not shown. 
The amplifiers, discriminators, and coincidence circuits are all of con- 
ventional design. The S,+8, coincidence circuit introduces 50 psec of 
inoperative time after each coincidence so that photomultiplier after- 
pulses from the 10° fast particles traversing the system each day cannot 
give delayed 8,+S8,+C, coincidences, which would evade the C, anti- 
coincidence pulses. 


A.D.- AMPLIFIER & DISCRIMINATOR. 
M - MIXER. L- LENGTHENER. 


Block diagram of fast electronic selection system, and timing-pulse display. 


The C, signals are delayed by 5x 10-8 sec with respect to the 8, signals, 
and then added to them. The §S, signal which is a step, and the C, signal, 
a pulse of opposite polarity, are amplified by the same distributed amplifier 
and displayed on a cathode ray oscillograph. Plate 5 shows some typical 
oscillograph records. The time base velocity is 5107 cm sec, and 

the distributed amplifier rise time is about 7x 10~° sec. 

The 8, light pulses are larger and more uniform in amplitude than those 
in C,. Displaying the S, signal as a step ensures that S, photomultiplier 
after-pulses are not added to the C, signal and do not introduce timing 
errors. A correlation between the timing calibration and the C, pulse 
amplitude was observed, and all timing measurements were corrected 
for this effect. 

A photograph of the oscillograph is taken for every coincidence 
S,+8,—C,+C,. Ifa cloud chamber picture is taken at the same time, 
the signals from C, are temporarily stopped and one oscillograph cali- 
bration picture is taken at the next S,+8,+C, coincidence. Since most 
of the S,+S,+C, coincidences are due to single fast mesons, these 
oscillograph pictures are used to determine the mean pulse separation for 
fast particles and to give an estimate of the error in this measurement. 
The time distribution of a set of 126 8,+S,+C, coincidences is shown in 
fig. 2(a). The r.m.s. deviation of a single timing observation in such a 


es 
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set Js about 310° sec. Unexplained instrumental changes over a 
period of hours sometimes increased the error on a single time determina- 


tion to as much as 4 10~® sec. 


OBSERVATIONS 


OF 


NUMBER 


Fig. 2 
(A) 
16 16 
12 12 
8 8 
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(a) Distribution of delays from 126 calibration signals from fast mesons. 
(6) Distribution of delays from 13 heavy mesons. The 5 events not unam- 
biguously identified are marked with arrows. 


The observed lifetime ¢, of an S-particle in its rest frame is given 
precisely by 


ax, dx, dita, L 

=D es lira 
where t,, ¢,, are the separation between the 8, and C, pulses for the 
S-particles and calibration particles respectively; 2%, 01, y1, are the 
position, velocity, and time dilation factor, for the heavy meson ; x2, V2 
the position and velocity of the decay particle ; L, c, are the mean track 
length and velocity for calibration particles. In fact the terms in the 
brackets contribute in no case more than 0:5X10-° sec to the interval 
t,, which \is much less than the standard deviation on a single timing 
measurement. 

Figure 2(b) shows the distribution in lifetimes of the 8 S-particles 
described in paperI. If we assume that all these particles are of the same 
nature, we find for the reciprocal of the mean lifetime, and its standard 
deviation 1/7=(0.063-- 0-022) x 10° sec™#, giving 7= (15-8744) X 10° sec. 
The timing inaccuracy contributes little to the error, which is mainly due 
to the small number of events observed. 

There are also five events in which for various reasons the identification 
as an S-particle is not unambiguous. We estimate, however, that not 
more than one of these is due to any other type of event. The delays 
associated with these five events are shown in fig. 2 (b). We att ribute the 
small mean delay of these five observations to a statistical fluctuation. 
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Thus we consider that the lifetime of (11-073j) 10° sec obtained by 
combining these five with the eight events previously referred to gives 
the best estimate of the lifetime. 


§ 4. Discussion 


Because this experiment demands a long distance of flight for unstable 
particles, it has a strong bias against detecting particles with lifetimes 
less than 5x 10-® sec. A small contamination of short-lived particles in 
the sample could not be detected by timing measurements, due to poor 
statistics, and probably would not be identified from the cloud chamber 
pictures. 

The lifetime deduced is consistent with the order of magnitude estimates 
obtained by other cloud chamber workers. The result is not inconsistent 
with the conclusion of Keuffel and Mezzetti (1954) that the lifetime they 
measure precisely (8-7 x 10~® sec) is that of the K,-particle. 

A notable observation of the experiment is that no timing observation 
on an S-event has yielded a lifetime longer than 40 x 10-* sec. This implies 
that if there are heavy unstable particles with mean lifetimes in the range 
10°-* to 10 ® sec they are not very numerous, since decaying particles 
with lifetimes in this range would be observed efficiently. 
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ABSTRACT 


; Radiative corrections to the magnetic moment of an outer electron 
in alkali and hydrogenic atoms due to the atomic electrostatic field have 
been calculated to second order. The correction is found to be 

dn 26 <V) 
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where{V) is the expectation value of the electron’s potential energy. 
For hydrogen this is 4/u=2-2 x 10-7’ —(2/7)«? and will be of importance 
in the evaluation of future experiments of improved accuracy. For the 
alkalis, however, the effect is too small to account for the observed 
Z-dependence of the magnetic moment. 


§ 1. INTRODUCTION 


ONE of the most fruitful methods for checking quantum electrodynamics 
against experiment has been the measurement of the magnetic moment 
of the electron. The recent experimental achievements (Koenig, Prodell 
and Kusch 1952, Beringer and Heald 1954) in demonstrating the fourth 
order correction to the magnetic moment of a free electron as calculated 
by Karplus and Kroll (1950) appear to confirm, in one respect at least, 
the verity of quantum electrodynamics. However, experiments (Franken 
and Koenigt 1952) on the magnetic moments of the alkalis reveal a dis- 
crepancy between the moment of a free electron and the corrected apparent 
moment of a bound electron. Although estimates of the effects of 
interconfigurational mixing (Phillips 1952) have suggested a possible 
source for the discrepancy, the uncertain nature of the results made it 
seem worthwhile to calculate the radiative corrections to the electron 
moment due to the atomic electric field. Such a calculation resembles 
closely that done by Kroll and Pollock (1952) on the effect of the Coulomb 
field of a hydrogen like atom on the splitting of hyperfine levels. 
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Since the problems of radiative corrections to the magnetic moment and 
to hyperfine splitting differ only in the form of the external potential, 
the method employed by Kroll and Pollock (1952) in arranging their 
calculation may be copied here. The notation of Kroll and Pollock will 
be retained and the reader is referred to that paper for details. 


§ 2. THE RADIATIVE CORRECTION 


In terms of Feynman diagrams, the radiative corrections to be 
calculated are as in the figure. 


vo os Sm 


Feynman diagrams for the second order radiative corrections to a bound 
electron. Double lines represent bound states. 


The external potential to be considered is 


Afr) =iV ([r])3,4—Hy8,1 Jo re ae 


ust 
where V(|r|) is the average isotropic electrostatic potential of the alkali 
atom excluding the outer electron, and —Hy6,, is the non-symmetric 


form of the vector potential of a uniform magnetic field, H, in the 
z-direction. Thus, taking the Fourier transform 


d 


e Pes. < 0. . ? 
jo bw DV (la))3.a— 7 HB ()8'(My)8(G2)8a ss (2) 


where 6(x) is the ordinary Dirac 5-function and 8’(x) is its derivative, 
le. [ 0’(x) f(x) du=—(df/dz),. 

The simplest integral, and one which will illustrate the main points 
in the evaluation of all other integrals, is the first one considered by Kroll 


* Boldface type will denote 3-dimensional vectors while ordinary type will 
denote four vectors. p?=p,?+ p,?+ p,2+ p,?. 

( V(|q|) is the Fourier transform of V(|r|). It is trivially truethat a function 
which depends only on |r| has a Fourier transform depending only on | q]. 


to the Magnetic Moment of a Bound Electron 313 


and Pollock (ef. Kroll and Pollock, p. 882). It arises from the polarization 
energy 


AH’ ,=2hea(27)"| $(p,)G(py, Pa)6(P:) dpidp. . . (3) 
where 


1 up\* (? 2 
G(P1, Ps) = z [dw 2v 3) dq 
y-A¢P2—Pi—Q)ity - (a +9) —K ly - Aq) 


x het) fob cr pa oe (4) 


Now, by observing that ¢(p) is large only for small values of | p |, Kroll 
and Pollock argue that a sufficiently accurate answer will be obtained 
merely by evaluating G'(0,0) instead of the full G(p,, p,). However, it is 
easy to see that for the case considered here G(0, 0)= co because one of 
the A’ in the integrand must be a 5-function while the other must be 
1/q?. Thus we must calculate G(p,, p,) near the origin. Nevertheless 
we can replace p in the bracketed expression by p=(0, 0, 0, ix), in which 
case G(P;, Py) becomes a function of only p,—p,. Using the notation 
Po—Pi=p; p—q=D we obtain 


e\2 (1 AD) y.q—e(l+y)}y-4@) 
ip)=( jg) [de 208000 [aa a 


Considering only terms linear in H, we have 
iy. A(D)y.qy.A(q)=ty4{ 44(D)[q. A(q) +20. q x A(q)] 
+A,(q)iq. A(D)—is.qx A(D)]}.. . (6) 
Neglecting for the moment the term in y. A(D)(1+y,)y.A(q) and letting 
F(q)=5(11,)8'(q,)3(q.) we obtain 
| dv 2v3(1—4v") [aq 


TEN 
ate) = (75) nit] | 


V(|p—a|)liq. F(a) —2 - (0,9,4(9)), = 9F(4)] 
<| te FP") 
Vijq|q.A(D) +e. (0, ¢.F(B)), —4,F(B)] \. 6 
| Ay? + q?(1—v?) 
Observing that the choice of H in the z-direction was arbitrary, we get 
f en pe” HV a 1+ terms in ne Pubs fe, (8) 
cp=—F5(je) ne-HV UPD a ae 
to lowest order in «2p2. The succeeding integrals, all of which are similar 
to this one, will also contain terms in p,? and p,p.. It is clear that these 
terms must cancel because they single out directions in the xy plane 
whereas we know the physical problem to be isotropic in the pepene. 
Similarly, it will be seen that the term iny. A(D)(1+ %) y ANG ) will ti 
rise only to terms proportional to p,/ and hence es poet ee 
Of course the terms in question involve y-matrices and hence factors m 
p,, but no combination of them can give isotropic terms, 


Y2 
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Thus, if j is the electron moment, we have 
Ap. OAL, | eh 4 
Ly pened =—=—-—- Sees ee Bm eee, oe ee 5 V . . . 9 
& ). oH ae ifm ) (a 
where (V) is the expectation value of the potential energy of the electron. 
This result follows at once from the observation that if V(p) is the Fourier 
transform of V(r) then 


” 


[8(P2)eV(P2—Pa)$(P1) IP: AP2=(1/27)9<V ). (10) 

Similar results hold for the rest of the integrals, the remainder of the 
calculation being briefly outlined in the appendix. In all cases one gets 
G(p)~V(p) which according to (10) implies du~<V ). 


§ 3. ResuLts AND DISCUSSION 


Thus, to second order the radiative correction to the magnetic moment 
of a valence electron in an alkali atom is* 
Ape. 20 eC iD 


: ia 150 me (11) 
For hydrogen this becomes 
A 26 i 
(=) kane rei a3=—2-2x10-* (for hydrogen), eral $e 


As regards the alkalis, a sufficiently accurate answer is obtained by 
using the approximate formula of Lamb (1941), (V)=W,)—3W where W 
is the ionization potential and W, is the ionization potential of a hydrogen 
atom in the corresponding n-state. It is then found that 


th 5 (=) (all alkalis), n{) oa, Semen a 
hoe ENG (ie lal 

This result is too small to account for the observed dependence upon 
atomic number of the electron moment in alkalis. Moreover, it would 
be rather difficult to observe in alkalis because the corrections due to 
relativistic effects, inter-configurational mixing and internal diamagnetism, 
all of which are much larger, and which must themselves account for the 
discrepancy, cannot be conveniently estimated with sufficient accuracy. 

The situation is completely different for hydrogen where the afore- 
mentioned corrections either vanish or are known with great accuracy. 
Nevertheless there is still an effect which has not yet been calculated, 
namely the sixth order correction to the moment of a free electron—an 
extremely lengthy calculation. However, it is not unlikely that its mag- 
nitude is less than the effect calculated here, and consequently a crude 
estimate, if such could be found, would enable the bound state correction 
to be observed experimentally. 

To measure the electron moment two facts are needed, namely the 
ratio of the proton moment to the bound electron moment, and the ratio 
of the proton moment to the Bohr magneton. The first is almost known 


* cf. the Appendix to this paper, 
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to sufficient accuracy (Beringer and Heald 1954), an improvement of four 
times in accuracy being sufficient. The second (Gardner 1951) however 
_ requires an improvement of 102. 

It is probable that the next few years will see a decided improvement 
in the experimental data, in which case the effect discussed here will 
have to be included in a discussion of the results. 
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APPENDIX 


Remainder of the Calculation 


The next term to be considered is eqn. (43) of Kroll and Pollock. 
Although for hyperfine splitting this contribution to the fluctuation 
energy vanishes, it is not so here The procedure is the same as for the 
polarization energy 


AB = —Gr*hicx | $(P2)Gn'(Pr: Pe) (P2) AP: AP | 


Gp®(Py, P2)~Ge(P) (A.1) 
we\? (2 Qui -y.qtK(1—2v—2vy,)] 
Therefore, discarding non-isotropic terms and considering only the lowest 
order in p?, ae 
(+) pees Ji 2 i ee) 
|p 7 mc 


We come now to the first term of eqn. (36) 


AH O=2n%hica.| $(P2)G x (P1> P2)#(P1) EP: Pe, oe 
, 2 «i < ; 
Gp(P 4; p)=(F) A,,°(Po—P1)K,(P 2: P1)> 


where K,,(p., 1) 18 displayed in eqn. (33). The second term oa K (Pe 1) 

may be split into two parts as done by Kroll and Pollock. The first part 

gives the anomalous moment of Schwinger while the second may be shown 

a, 

to vanish. The other terms in K,, (2, P1) may be expanded using eqn.(23) 

and after some manipulation one can show that 
ee. (4.4) 
pw Je am mc 
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It should be pointed out that in obtaining this result a term involving 
(A,,°)? was dropped. This is certainly justified because the momentum 
is small compared to the rest mass. 

We now come to the second order parts, Lp) and Y. The justification 
in Kroll and Pollock for using only the first approximation for S,,° involves 
the nature of the magnetic potential. Such an argument fails in the case 
of a 6-function, but we can ignore higher order terms in S,,° on the grounds 
that terms in (A,°)? make only small contributions. It can then be 
shown that L, vanishes because of the symmetry properties of the 
integrand. 

Finally we must evaluate Q. 


Q=2inex| $(P2)Go(P)$(P) dP; dP 


(pa (2 ee AD ee eae mae ( 
cale)—(F) pe uD, Db) oq) bp ge ewe Pk) day (AS) 
: | 
ve \2 rs es | 
-(z) (I-+4II-+II1), | 
where I, I, III are the integrals indicated on p. 886 of Kroll and Pollock 
with A‘(—q) replaced by A‘*(D). As it turns out, there is no need for an 
infra-red cut-off in these integrals because they converge independently. 
After a considerable amount of manipulation, one can show that 
k,?-++ Pe ah 
k2(k®+-2ink,)> 
Thus Y=0 and on adding up all the terms we arrive at the final answer 
Au 26 <V) 
pal eh aes UES 


pe 15a me? ° 


(I+II-+II1)=8V (|p |)o. H dae 
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XXXVIT. An Exact Solution of the Spherical Blast Wave Problem 


By J. Lockwoop Taytor 
Norges Tekniske Hogskole, Trondheim, N orway* 


[Revised MS. received September 7, 1954] 


Various authors have investigated the problem of an explosion wave in 
air, assumed to start at a point, during the stage when the spherical 
shock wave is still strong, and some numerical results and at least one 
approximate analytical solution are known, but the present is, so far as 
is known, the first exact mathematical solution to be published. It is 
obtained by straightforward deduction from the known general equations 
of motion. 
The equations of motion can be written 


1 A(pua®) Op _ 


ee SO ee ee 
Ou ou ON =: 
1 (dp OMY By ($ Op 
a(#+ ne) =2 Poe cee (3) 


expressing respectively (1) continuity of flow, (2) radial force=rate of 
change of momentum, (3) adiabatic pressure—-density relation for a 
particle. By combining these three equations, or otherwise, we get the 


known energy relation 


a hes £ (aul) =0. oa Lee ene) 
The solution we are concerned with is of the form 
p=f i (a/t?/>), Le ie BRS > a een 
i tele wire), eu 1 S(O) 
aba) °F, (g)t""*), Se oot ek Se (7) 
with R=kt?/® defining the shock radius at time ¢. 


; Nw 
Thus the shock velocity U= ar 
It can be seen that (5), (6) and (7) make both the kinetic energy 


“R 
der | tourx* dx 
0 


5) 
kt-8/9, 


and the internal or intrinsic energy 

Lua als 
— [ pu da 
y—lio 


within the sphere constant at all times, also satisfying the shock} conditions 


* Communicated by the Author. 
+ See Appendix. 
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The specified dependence on ¢ makes the problem in effect one of a 
single independent variable (x). Thus 


Op __ 2 x Op _ 2 7% 
ot }08=—sSiB é Ox R~ ox’ 
Op... 2: Op =) Oras © OP 37 
Ot, a web orueGut BA Oral 


Equation (4) together with the relation 
oH 6H .2x0H 


a=+tse-tere =O... 2.2. 8 
Of 5 2et cB hoe (8) 
immediately following from the form H=t-*f(x/t?/*), gives 
CLL Sa IE ee ean Soe ee 
Be ie nat 
Fi ae Ye le, ew Bye a ee 
fe URI So 
which we rewrite as 
Tey P 
Uu =x ————— }. ot et ee ee ( ) 
ae ts 
2° * (v—l)p j 
Solving, 
a daca x a te ue sia), (12) 
p 2 yu'—x yu'—x 


ldp y—l ae 1 Dee ae y—lodp lou 2 (13) 
p 0x p 0x up at tp Ct WM Or. ae ree 
On substituting for dp/dt and 0p/ét and dropping the partial differentiation 
notation 


a  u\(ldp y—ldp 3 ldu 24 14 
Heo) @ae 2 at pon oa 


or 
1 du 
ld —1 Se 
LOD, a een Eee Ce ae fe STEN 
p dx p dx a x/R—uU } 
Integration gives 
p , , 
ora or —wu’)-ly'-2, (16) 
By alternately eliminating p, p from (12) and (16) we get 
p=C,[xu(x—u)]-?/2-) (yu—a) V2 iS ee oe 
p=C gt 2/2 yA (qq— yA (yy — J P- DI) (18) 


where x, « are now conveniently written for 2’, w’ (equivalent to taking, 
for the present U=1, R=1 which does not involve any loss of generality). 


} 
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Differentiating log p from (17) with respect to x and equating to 


Ledpe du 2u 
pdx” (x—u)\dx ' a 


(19) 


as can be deduced from the equation of continuity (1), gives a differential 
equation for uw in terms of « which need not be written out, but whose 
solution is 


2 log u=a log «+6 log (x 2 u) +h log (yw—ax)+logC (20) 
where 
ey 1) 
Ee 
b= (— 1377+ Ty—12)/(8y—1)(2y +1), 


yy fae 
h=5 (5) , 


The result (20) together with (17) and (18) express u, p, p in terms of 
x and therefore give the required solution. The constants C, C3, C, 
are to be found by substituting the known values of the variables at the 
shock (x or x’=1, u or u’=2/(y+1)). 

The exact solution as described allows for instance the accurate 
calculation of the value p/p, at the centre of the sphere, otherwise difficult 
to obtain precisely. The calculated value is 0-303 for y=5/3 and 0-365 
for y=1-4, as compared with Kynch’s value of 0-31 for y=5/3 and Sir 
Geoffrey Taylor’s of 0-374 for y=1-4. 


a= 


Symbols 
x—distance from centre of sphere. 
p=air density. 
u=air velocity (radially outwards). 
p=air pressure. 
¢=time from explosion. 
y=ratio of specific heats. 


R=shock radius=kt2”°. 
Suffix s denotes shock. Suffix 0 denotes free-air. 


l p 
E/p=total energy of air per unit mass= {5 u?+ (=1)e ' : 


u YP 
I/p=total heat of air per unit mass= {5 u?- = wai : 
2 : 
U=shock velocity = 5 ht ea a 
eel te aU. 
M=Mach No.=U/[¢o, (=, | — 
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C", Cy, Cy, C5, Cy, k, a, 6, h, constants. 
The values of C,, Cy, C3, Cy, got by substituting for the known 
pressure density and velocity at the shock are : 
1 Nae i | 
Co= 29? "(y—D[(y +1} 0?, 
C= 22/27 9 (yf 1)-V+D/2-M (y— 1)0-D/2), 
CO p= 2?!2-V og (y+ 1)-PtD/2-1) (y= 1) VAN G2, 


BP PN Dax 
Conditions at Shock 
When the shock is sufficiently strong, the usual shock equations give : 


==] ae 
M = if = aa U—u,)/¢,, C.= VPs Ps 


Yack Vail. 
») 2 
RS ap) 72 
Ps 41 Po if eee } 
Numerical Values for y= 1-4 

a/R 0 0-6 0-8 0-9 0-95 1-0 
u/U 0 0-429 0-584 0-685 0-7513 0-8333 
D/Dg 0-365 0-366 0-4035 0-510 0-651 1-0 
ol, 0 0:00675 0:0654  0:2056° 00-4106. ——s-0 
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XXXVI. The Measurement of Photoneutron Yields with a Sodium 
Lodide Crystal 


By K. G. McNemn 
Department of Natural Philosophy, The University, Glasgow + 


[Received December 8, 1954] 


ABSTRACT 


It is shown that a Nal (T1) crystal will act as a detector of photoneutrons 
by virtue of the production of radioactive 1281 within the crystal. This 
detecting technique has been applied to the comparison of the photo- 
neutron yields obtained from copper, cadmium, mercury and lead when 
these elements are irradiated with a bremsstrahlung x-ray beam of 
maximum energy of 22 Mey. The yields from these four elements are 
found to be in the ratios Cu: Cd: Hg: Ph as 1: 4:1: 9:5: 8-4, 


§ 1. INTRODUCTION 


THE importance of photodisintegration has given rise to the development 
of many different techniques for investigating this process. When the 
source of X-rays is a pulsed machine such as a beta- or synchro-tron, the 
instantaneous high intensity of the x-ray flux militates against the 
immediate detection of the product particles. Asa result, most of the work 
in this field has been done using techniques in which the measurements 
are taken after the end of the beam pulse. An example of sucha technique 
is that based on the production of radioactive isotopes (Johns, Katz, 
Douglas and Haslam 1950). A more recent method (Halpern, Nathans 
and Mann 1952) uses the fact that thermal neutrons have on the average a 
free life of a few milliseconds before capture by nuclei, and that therefore 
the photoneutrons, after slowing down in some moderator, may be 
detected by conventional methods at a time much later than the end of the 
X-ray pulse. 

In a previous note (Reid and McNeill 1953) it was observed that a 
background activity was induced in a Nal (Tl) crystal which was being 
used to detect x-rays associated with the *°°Pb (yn) 7°’*Pb reaction. A 
study of this activity has now shown that it is caused by the interaction of 
fast photoneutrons with the iodine of the crystal. This interaction 
produces #281, which decays by 8 emission with a half-life of 25 minutes. 
As this activity is actually produced within the crystal, the detection 
efficiency for the resulting B-rays is nearly 100%, and the Nal crystal may 
therefore be used as a reasonably efficient fast neutron detector. 


+ Communicated by the Author. 
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The present paper describes some measurements of relative cross sections 
for the (yn) process in different elements, using this fast neutron detection 
technique. In principle this technique resembles those previously 
quoted in that the measurements are of the integrated effects of the beam 
pulse and are carried out well after the end of the beam pulse. However, 
as the Nal crystal is activated by fast neutrons, the photoneutrons do not 
need to be thermalized prior to detection. 


§ 2, PRELIMINARY EXPERIMENTS 


During the initial experiments, different targets were bombarded by a 
bremsstrahlung x-ray beam of maximum energy 22 Mey. A 1-0 inch 
cube Nal crystal, mounted on an EMI 5311 photomultiplier, was placed 
under the targets, but outside of the collimated beam. The radioactivity 
produced in a standard copper foil was in all cases used as a beam intensity 
monitor. With all the targets used, a radioactivity was induced in the 
Nal crystal, and the half-life of this activity was found to be independent 
of the target material. The average of all results gave the half-life as 
25:4+0-2 minutes. Moreover, it was found that for a given target the 
strength of the induced activity was proportional to the x-ray beam 
intensity. These facts show that the activity is due to *I produced by 
the reaction 1*7I (ny) ?8I, and, therefore, that the intensity of the 
activity may be used as a measure of the neutron flux through the crystal. 

A considerable activity was induced in the Nal crystal when there was 
no target in the beam. This was considered to be due to photoneutrons 
produced at the internal target and do-nut walls of the synchrotron, 
and at the walls of the beam room. ‘To reduce this beam background the 
detecting system was screened by placing it within a large cubical paraffin 
box. The box had an external dimension of 78 em, and had inside it a 
chamber lined with cadmium and iron. Inside this cavity were placed 
the various targets and the crystal, photomultiplier, and head amplifier. 
The centre of the detection cavity was 122cm from the target of the 
synchrotron. Slots were cut in the box so that the x-ray beam could 
pass through it. 

This screening arrangement was used in all the work subsequently 
described. 

§ 3. ENERGY OF THE PHOTONEUTRONS 

Some knowledge of the energy of the neutrons reacting with the iodine 
nuclei is essential for the comparison of neutron yields from different 
targets. As the (ny) cross section in iodine has been measured and found 
to be markedly dependent on neutron energy, the absolute determination 
of this cross section in a particular experiment would enable the energy of 
the incident neutrons to be found. Photoneutrons are known to be emitted 
with a whole spectrum of energies, and therefore in the case of photoemission 
the experimental determination of the (ny) cross section will lead to a 
mean cross section defined by | N.o,de/ | N.de, where NV, are the numbers of, 


ru 
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and o, the (ny) cross sections for, neutrons of energy «, the integral being 
taken over all available neutron energies. ‘To this mean cross section will 
_ correspond a mean interaction neutron energy. 

To determine this cross section it is necessary to bombard the Nal 
crystal with a known flux of neutrons. As the 22 mev Cu (yn) yield is well 
established (Byerley and Stephens 1951), this known flux was produced by 
irradiating with the x-ray beam a copper target placed at a known distance 
above the crystal. The absolute strength of the beam was measured with a 
Victoreen meter, and the copper monitor activity normally used for 
measuring the beam intensity was calibrated in terms of this Victoreen 
meter (Wright 1954). 

In a typical experiment, after 15 minutes irradiation of the copper 
target, the beam was switched off and the decay of the induced activity 
followed for 30 minutes, readings being taken every 5 minutes. The 
rate of decay at switch-off was determined using the measured half-life 
of 8J, and thus the rate of production of !8I found. The average of 
many runs was taken, and allowance was made for the fact that with only 
15 minutes irradiation the iodine activity would not be saturated. 
A fuller discussion of the experimental technique is given in §4. From the 
rate of production of I, the known flux of neutrons, and the knowledge 
of the number of iodine nuclei present in the Nal crystal, the cross section 
of the !2’I (ny) 1281 reaction was calculated. 

This cross section was found to be 1-6--0-5 x 10-*° cm?, corresponding to 
a mean interaction energy of 0-979 mev (The Science and Engineering 
of Nuclear Power, Vol. I, 1947). 

Byerley and Stephens (1951) measured the energy distribution of Cu (yn) 
neutrons, and from their results an average neutron energy of 2-0 Mev is 
indicated. The mean interaction energy will be less than this, as the 
low energy neutrons will give a disproportionately large yield. An 
estimate has been made of the mean interaction energy using Byerley 
and Stephens’ results and the known 12’I (ny) cross sections, and this gives 
that the mean interaction energy equals 1-2 Mev, in reasonable agreement 
with the experimental result quoted above. 

The fact that resonance neutrons, with their much higher cross section, 
play no part in the interaction was checked by surrounding, the Nal 
crystal with a layer of iodine during the course of an irradiation, and 
removing this layer before counting of the induced activity was commenced, 
The iodine layer was sufficiently thick to absorb 80% of any 40 ev 
resonance neutrons, and at least 16° of any neutrons below this resonance 
energy. To within the experimental error of 5%, this iodine layer had ne 
effect, indicating that in fact only fast neutrons were effective in the 
interaction with the iodine of the crystal. 

These facts show that the neutrons reacting in the crystal are coming 
straight from the target, and have not been thermalized in the paraffin 
sereening prior to capture. Apparently the neutrons have an opportunity 
of interaction when first they pass through the crystal, but if once they 
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enter the paraffin, they will be slowed down there and the cadmium sheet 
will prevent them re-entering the detection chamber. 


§ 4. COMPARISON OF PHOTONEUTRON YIELDS 


As the neutrons interacting with the iodine of the crystal have come 
directly from the target without losing energy by collision, the distribution 
of energies of the photoneutrons for a given maximum X-ray energy will 
depend only on the shape of the (yn) cross section curve of the target 
material. If the shapes of this curve for different materials are the same, 
it is possible to say that the ratio of the induced iodine activity for two 
materials is equal to the ratio of their (yn) yields, for although the !?’I (ny) 
cross section varies with neutron energy, the relative proportions of neu- 
trons in each energy range will be the same for both materials. Yield 
curves have been found for many nuclides, and the great majority of these 
show the same ‘ giant resonance’ shape. In particular Cu and Pb have 
similar characteristics (Montalbetti, Katz and Goldemberg 1953), and it is 
reasonable to assume that the other two elements investigated, cadmium 
and mercury, have the same general characteristics, especially as the 
multi-isotopic characters of these two elements will smooth out any 
peculiarity which might be present in any one nuclide. 

It is therefore possible to find the ratios of the 22 Mey (yn) yields by 
comparing the iodine activity induced when different targets are irradiated. 
Experiments were performed to compare the yields of the four elements 
quoted. In all cases thick targets were required to give adequate yields, 
the thicknesses being chosen so that equal beam absorption occurred in all 
targets. 

Normally the crystal was left a sufficiently long time between runs so 
that it was quite inactive before the next run was started. If this was not 
possible, allowance was made for the residual activity. As in all cases 
background had to be taken into account, the normal procedure was to 
subtract ‘cosmic’ background, then normalize to a standard beam 
strength using the copper monitor activity, and then to subtract the 
normalized ‘no-target’ background. The result was taken as the 
normalized yield due to the target. 

The * no-target ’ background was in all cases an appreciable fraction of 
the total yields. In the case of the copper target, for instance, the ‘ no- 
target ’ yield was almost a third of the total yield obtained with the 
copper target in position. Higher yields were obtained from the other 
targets, however, and therefore the ‘ no-target ’ yield was relatively less 
important. The causes of the ‘ no-target ’ yield were probably the stray 
neutrons which can enter the detection chamber via the beam channels, 
and the (yn) reactions that may take place in the detection chamber 
screening itself, 

The comparative yields obtained from the four target materials are 
shown in the table in two ways : in column 2 is given the normalized yield 
from each target, and in column 3 the 22 mev yields per mole are given 
relative to the yield from copper. The statistical errors are indicated 
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after each result. For comparison, it may be noted that the ‘ no-target ’ 
yield was 143 counts in the first 30 minutes after switch-off per 1000 Cu 
monitor counts. This background has, of course, been subtracted from 
~ the total ‘ with-target ’ yields before inclusion in the table. 


———SS__—E ee ee eee 


Target Counts in 30 minutes 22 Mev yield/mol/r | Yield/mol/r 
element per 1000 monitor counts relative to copper x 10-8 | 
Cu 288 +15 A 1-0 3.2. 

Cd 647 +28 4-140:3 13 
Hg 661 +26 9-5 +0-9 os 30 
Pb 470 +17 8-4 0-5 27 


§ 5. Discussion 


The ance yields of lead and copper have been measured by 
Montalbetti, Katz and Goldemberg (1953). From their graphs the 22 mev 
yields of these two elements are seen to be in the ratio of 9:1, a figure 
which is in satisfactory agreement with the present result of 8-4: 1. In 
the cases of cadmium and mercury, direct comparison with previous work 
is not possible as no data on these elements is available. 

Photoneutron yields at 22 mev are known for many elements of widely 
differing atomic mass (Nathans and Halpern 1954, Montalbetti, Katz and 
Goldemberg 1953), and using these results a graph can be drawn of yield 
against atomic mass. It is found that a smooth curve may be drawn 
through the experimental points, and the results of the present investi- 
gation may therefore be checked by plotting them on this same graph. 
Before this is done a photoneutron yield must be assumed for one of the 
four elements studied, since the present experiment gives only relative 
yields. The present results for Pb have a smaller statistical error than do 
those for Cu, the other element for which independent figures are available, 
and therefore Pb has been taken as the standard. The value of 27 10° 
neutrons/mol/7, taken from the graphs of Montalbetti, Katz and 
Goldemberg (1953), has been accepted, and using this figure the yields of 
the other elements were calculated. These values are given in column 4 
of the table. When these yields are plotted on a Yield~Atomic Mass graph, 
as suggested above, it is in fact found that they fit very satisfactorily 
with the results of other authors. 
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ABSTRACT 


When a vertical or inclined nozzle is fixed in an open tank with its 
outlet flush with the surrounding water, the jet from the nozzle is irregu- 
larly disturbed by waves due to the fall of the jet. Small-scale experiments 
in a circular basin showed that, if the nozzle were properly sited and 
trained, and supplied at the correct pressure, the water in the basin 
could be thrown into one of its natural states of stationary oscillation. 
The jet then pulsated in an almost steady manner that was much more 
interesting to contemplate than an irregular or monotonously uniform 
jet. A simple, approximate theory is put forward to explain the main- 
tenance of the self-excited oscillations in the basin. 


$1. INTRODUCTION 


BINNIE AND HooKines (1948) showed that, when water is supplied 
tangentially to a circular tank and allowed to escape in a whirlpool 
at the centre of the bottom, the flat free surface in the tank is sometimes 
disturbed by oscillations. These form a system of waves that travels 
unchanged round the centre, and they appear to be due to resonance 
between the natural periods of the water in the tank and of disturbances 
in the whirlpool. But the nature of the instability of a whirlpool under 
gravity is obscure, and the arrangement cannot be used to excite a 
predetermined natural period of the tank. The present paper describes 
another way of producing self-excited oscillations ina tank. Experiments 
were carried out at the Engineering Laboratory, Cambridge, with a small 
nozzle fixed with its outlet level with the surrounding water. It was 
found that, if the trajectory of the jet were properly arranged. the fall 
of the jet influenced the submergence of the nozzle in such a way that 
many of the natural stationary periods of the tank could be excited. 

The kinds of free oscillation that are possible in a circular basin must 
first be examined, and an account of the most important is given in § 2. 
It seems unlikely that the action of the rising jet on the wave over the 
nozzle outlet can have any exciting effect. For the relative velocity 
of the jet and the wave is greater as the crest of the wave subsides than as 
it rises; hence, if the nozzle is not submerged for the whole of every 
cycle, the average force on the wave over a eycle opposes the motion. 


* Communicated by the Author. 
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Aun explanation of the oscillations will be sought by investigating the 
trajectory of the jets, and a simple theory is developed in §3. The 
experiments are described in § 4, and when they are analysed, the main- 
tenance of the waves can be accounted for. Presumably they were set 
going by some chance combination of the random disturbances which 
were always seen to occur. The work has an application to the design 
of pulsating fountains, and it was in fact suggested by an examination 
of the fountain of this type that may be seen at Villeneuve, Switzerland. 


§2. Tor THEORY OF FREE STATIONARY OSCILLATIONS 
IN A CrrcuLAR BASIN 


For a circular basin having uniform depth and vertical walls the’ 
complete theory of small oscillations was worked out by Rayleigh (1876), 
a summary of whose analysis was given by Lamb (1932). For motion’ 
possessing circular symmetry the waves have annular crests and hollows,’ 
and the shapes of the’ first four modes are displayed in fig. 1, which is” 


Fig. 1 


0 0.25 ans 0.75 l 


Profiles of symmetrical waves. 


drawn to an arbitrary vertical scale on a base of R=r/a, where r is the 
radius at any point and a is the radius of the basin. The radius and depth 
of the experimental tank were 18 and 9 in., and for these values the upper 
part of table 1 shows the periods + of the waves of fig. 1. 


Table 1. Periods of Symmetrical Oscillations and of Asymmetrical 
Oscillations with one Diametral Node. g=32:2 ft./sec® 


Symmetrical 
mode 1 2 3 4 
T Sec © 0-708 0-513 0-425 0-372 


Asymmetrical 
mode i i ili 
T SeC 1-173 0-590 0-464 
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Experiments were carried out not only on symmetrical oscillations 
but also on the species for which there are one diametral node and 
6, 1, 2,....cireular nodes. Cross sections of the first three waves, 
taken through the centre at right angles to the diametral node, are given 
in fig. 2 and the periods in the lower part of table 1. Only in the funda- 
mental mode (i) was the period appreciably influenced by the depth. 
In this mode the water sways from side to side in a simple manner. 
The second mode (ii) is more complicated, and the contoured plan given 
by Lamb on the lower half of page 288 shows the important feature of 
a mound at R—0-345 accompanied by a hollow on the other side of the 
diametral node. 


Fig. 2 


0 0.25 ys 0.75 


Profiles of unsymmetrical waves. 


No exact solution is known for a circular basin of variable depth, 
but by means of Rayleigh’s principle an approximation may be found to 
the period of the slowest mode in a hemispherical basin. As an approxi- 
mate type we assume that the free surface remains plane. Lamb (1932) 
shows that, for a liquid of density p, a semicircular element of radius } 
and thickness dx possesses kinetic energy 7’ given by 


a : 
T= (- — 3) posBedi. 4 2 sete meee ch) 
a  & 
and potential energy V given by 
Vetn pbk dae 2. eye (2-2) 


when the plane surface is inclined at an angle B to the horizontal. We 
will employ these expressions in the manner indicated in fig. 3 which shows 
a plan of a hemisphere of radius a. It follows from (2.1) that the kinetic 
energy of the hemispherical mass of liquid is 


2 7 5 49, n/2 5 § dd—0-2602 a6? , (2.3) 
2(-—< cos =—0'2602pa"6*, 
: € 3) as I, 
and from (2.2) that the potential energy 1s 
71/2 : d 
gpa | cos! § dd=0-3827gpa4B?. . . . . (2.4) 
0 


Z2 
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We assume that Boc cos ot’, and the method yields 
__ 0-B827 g 


2 t 2.5 
o= 0-2602 a’ (eh 
hence the period 7 is given by 
27 a\ 2 
r= =o (<) Oe hee 
g g 


It is known that this estimate is slightly too small. To check (2.6). 
a hemisphere of radius about 52 in. was placed on a table and filled 
with water almost to the brim. When the table was momentarily 
shaken, it was found that the mode of oscillation contemplated could 
in fact occur and that its period was 0-60sec. The value calculated 
from (2.6) is 0-62 sec, showing reasonable agreement with the experiment. 


Fig. 3 


Plan of hemispherical basin. 


For a circular basin of radius and depth both equal to a, the gravest 
period is 7=4-738(a/g)"?, which is smaller than (2.6) because the kinetic 
energy is less for a given displacement, while the potential energy is the 
same. But the effect of rounding off the sharp corner in the circular 
basin is not large. Many fountain basins have a cross section of shallow 
elliptical form, intermediate between rectangular and hemispherical. 
An approximation to the gravest period can be made by means of 
Rayleigh’s theory (1876) with the depth taken to be the maximum 
depth. This estimate will be somewhat too small because the kinetic 
energy is greater than when the depth is uniform. 


$3. THe THEORY OF THE TRAJECTORY OF THE 
PARTIALLY OBSTRUCTED JET 


In order to discover how the impact of the falling jet stimulates the 
oscillations we must know when and where a jet from a submerged nozzle 
will impinge. We consider a nozzle delivering, when unsubmerged, 
a jet of velocity V, at an angle « to the horizontal and fixed with its 
outlet exactly at the mean level in the basin. Time ¢’ will be measured 
from the instant when submergence of the nozzle by the rise of the 
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neighbouring surface has just begun. Submergence lasts until t'=7/2, 
and it will be assumed that the minimum jet velocity, denoted by V,, 
occurs at t’=7/4. Nothing seems to be known about the velocity and 
diameter of the jet cast into the air from a submerged nozzle, even when 
the surrounding water is quiescent. Here we will suppose that refraction 
is negligible and that the jet velocity v varies in the sinusoidal manner 


v=Vi{1- (1 pi) sin 2a Orn Ye ee bev cond) 
a eS 


At i’=0 the surrounding surface is moving with its maximum velocity, 
and (3.1) is arranged to give the greatest rate of reduction at that instant. 
We will now follow the trajectory of a particle leaving the nozzle at time ¢’ 
on the supposition that its path through the surrounding water is neg- 
ligibly short compared with its path through the air. Air resistance 
and interference in the air between different parts of the jet being ignored, 
the required relations are obtainable by very elementary methods. 
If the maximum and minimum heights attained are denoted by H, 
and /7, and the maximum and minimum ranges by #, and #,, 


(V, sin a)? 


== 3045) 
A= 2g (3. 2) 
Vee a 
ofp Sino De eae Oe Tie ee Sas) 
gJ 
and the ratio V,/V, can be found from either 
i (FE 1/2 ; 
Va co(22| 2 0 Saar te ere 
V; H, 
V. RN 1/2 : 
— = aA . . . . . . . . ono 
= Vy (F) oy 
Corresponding to a jet velocity v, the range is 
2 
Gi sila ee 
J 
and the time of flight is 
v Q 7 
2-sin o. Me ee 8. BAD) 
Y 
We now introduce dimensionless quantities by writing 
F' i 
aN ae ‘ L— —, SPS, sn ye ee 
Y ge Fi, and t a (3.8) 


Then from (3.1), (3.3) and (3.6) the non-dimensional range is 
( loa. lg <a 
1 fe . 1 =2) sini eee se = (8.9) 
V5 J 
and with the origin of time mentioned above, the non-dimensional 
time of striking the surface is found from (3.2) and (3.7) to be 


2 /2H,\12 ( r) — \ 
=t+ —{— 1— {| 1— =} sin 2at;, pe a ted) 
nats (NY i? 


/ o 


. 
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the factor of the term in twisted brackets being the non-dimensional 
time of flight of the unobstructed jet. The eqns. (3.9) and (3.10) hold 
good for 0<t<z/2, and for the remainder of the period we have merely 
to take V,/V,as unity. When H,, 7 and V/V, are known, the procedure 
is to insert in (3.9) and (3.10) a series of equidistant values of ¢, and the 
results, when plotted together, show the nature of the fall of the jet over 
a complete period. For a vertical jet we put sin 2=1 and omit all 
mention of Z#. 

Occasionally it was noticed that the jet was unobstructed for appre- 
ciably less than half a period. To meet these conditions we will suppose 
that submergence lasts for 37/4, beginning at t/=—7/8 and ending at 
t’=57/8. Again we assume a sinusoidal variation in velocity, so that 
(3.1) is replaced by ; 


Ake a 5 
r=v,f1—(1— 72) sin F(Z +5). £0 SENSE EY 


Equations (3.2)—-(3.8) are still applicable, and (3.9) and (3.10) become 
f Va\.. -4er 1\)2 
ida ce) es a ie Se ligne 
R 2 (1 7) sin Ga :)} ; (3.12) 
2 (2H,\v2 V\ 5. oh 1\) 
—=t-+ -—{— —{l——} = —}>. 5. (pjyl ts: 
and t. += / ; ) {1 ( r) sin 3 (1+ =)} (3.13) 


§4. DESCRIPTION OF THE EXPERIMENTS AND COMPARISON WITH THEORY 


A symmetrical arrangement of three jets, one of which is shown in 
fig. 4, was fitted at the centre of a circular tank of diameter 36 in. and 
depth 9in. The nozzles were made of glass tubing drawn out to about 
0-15in. diameter. It was soon discovered that a single nozzle was 
sufficient to excite oscillations. Thereafter, usually only one was used, 
but the others served as gauges indicating whether the oscillation was 
symmetrical or not. For operation in other positions in the tank, the 
nozzle was mounted on a clamp and stand. The water escaped from 
the tank near the periphery through a vertical pipe which terminated in 
a small circular weir on a screw thread; thus by rotation of the weir 
the water level could be finely adjusted. 

In experiments with an inclined jet the usual method was to set the 
nozzle just clear of the surrounding surface. After the maximum height 
H, of the trajectory had been measured, the water level was raised by 
means of the weir until disturbance of the jet began; at intervals of 
some minutes the level was raised further in small steps, ample time 
being allowed for pulsations to make themselves apparent. If steady 
oscillations set in, they were timed by watching the nozzle, observing 
the top of the trajectory, or listening to the rhythmic ‘ flop ’ that’ often 
occurred, due to the simultaneous impact of a large part of the jet. 
The surface in the tank was violently disturbed by the falling jet, which 
threw up random waves of small length, and it was impossible to discern 


ee 
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“any regular oscillation at the periphery, still less the position of nodes. 
“The general commotion made the action of the submerged nozzle irregular 
: and the accurate observation of H, and #, a matter of extreme difficulty 
‘when the trajectory was long. No attempt was made to determine 
‘the wave amplitude, which was probably of the order of 4in. near 
‘the nozzle. When once the correct setting of H, and of the circular 
‘ weir had been discovered, the pulsations set in quickly after the apparatus 
had been started from rest, the maximum time required being not more 
than a couple of minutes. On no occasion were the oscillations stimulated 
_artificially. 


Fig. 4 


Elevation of nozzle and support. 


In the vertical tests, the nozzle was set slightly out of the plumb in 
order to reduce interference between different parts of the jet. Never- 
theless, the water that rose to the minimum height was often struck 
there and dispersed by the free jet following up behind. It proved 
possible to excite symmetrical and asymmetrical oscillations with both 
inclined and vertical jets, and the experimental results can be conveniently 
divided into the four groups shown in tables 2 and 3. In all the tests 
there was no ambiguity about the mode, for the observed periods were 
close to one of the theoretical values set out in table 1. 

In the first group—oscillations with circular symmetry excited by an 
‘inclined jet—the nozzle was held nearly at the centre and arranged to 
‘strike (when unobstructed) close to the periphery. tn the ao 
mode (1), the nozzle was seen to be submerged for about half the oe 
‘and so (3.9) and (3.10) apply. They have been used with the a is 
‘in table 2 to produce fig. 5, in which # is plotted upwerds on - ase 
‘of t,.. The-diagram begins at (= eee ea R= 1, and eaot im oe 
girection of the arrow. For 0<t<1, the points on the curve were 
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calculated at intervals of ¢ equal to 1/20; and since the time intervals 
are the same, approximately equal volumes of water impinged between 
each pair of adjacent points. Figure 1 indicates that the waves at 
the centre and at the periphery were in anti-phase, and the node is marked 
on the left-hand vertical axis in fig. 5. It was noticed that a concentrated 
fall of low-velocity water took place at the minimum radius, preceded 
by a quick fall from the periphery inwards ; the subsequent fall outwards 


Table 2. Experiments with Symmetrical Modes 


Mode “ Hy H, A; ue z (3 
sec in. in. in. i ir \ tg 
Inclinejd jet 
I 0-714 SS a= 7 0-62 1-20 
v4 0-520 8 _- 8 0-67 0-79 
3 0-432 I3 — 8 0-67 1-21 
4 0-380 1h _ 4) 0-71 1-57 
Vertica)l jet 
1 0-709 61 24 — 0-63 1-59 
2 0-508 33 13 — 0-63 1-63 
3 0-432 22 9 — 0-64 1-57 
4 0-380 17 7 — 0-64 1-57 


0 0 
5 | t. 1.5 2 


Fall of jet from inclined nozzle at centre ; circular symmetry, one node 
and three nodes. 


along the radius seemed very light. The diagram explains the first 
and third of these observations but not the second; the trend of the 
initial part of the curve, instead of being to the left, should have been 
slightly to the right. Thus the cut-off of the jet by the rising wave was 
not as sharp as postulated in the theory. The minimum radius was rather 
variable, but whatever the exact timing of the above events may have 
been, they occurred across the node not far from ¢,=1 when the surface 


EE 
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was flat, hence their effect on the oscillation could not have been large. 
The main stimulus seems to have come from the unobstructed jet which 
struck in the interval 1-7<t,<2-2. At t,=1-75 the wave at the periphery 
-Was at a crest, hence a powerful force was exerted on the descending 
wave. Points in fig. 5, for which ¢, differs by unity, are simultaneous, 
therefore the low-velocity concentration struck while the high-velocity 
fall was in progress, and it is not surprising that a ‘ flop’ was heard. 

The second mode (2) was remarkably regular and easy to obtain. 
Conditions were now different because the waves at the centre and at 
the periphery were in phase. The corresponding #—t, curve is given 
in fig. 6 with the nodes marked on the left. The obstructed jet was seen 


Fig. 6 


T T T = T T i 


05 5 2 
: t, 


Fall of jet from inclined nozzle at centre ; circular symmetry, two nodes. 


to fall in much the same manner as before, and now the initial part of 
the curve is nearly vertical, in better agreement with observation. 
At t,=0-75 the region between the two nodes was at its maximum elevation 
and the motion was assisted by the concentrated fall of water in that 
area at that instant. At ¢,—1-25 the wave at that periphery was at 
a crest, and again the unobstructed jet fell on a subsiding wave. It 
should be added that with the nozzle fixed near the periphery and trained 
on the centre a short attempt to produce the gravest mode failed when 
the jet was arranged in accordance with the first line of table 2. Instead, 
the second mode appeared at H,—36in., #,=9 in., to which the corres- 
ponding time of flight is 1-63. Thus the unobstructed jet was in ube 
air for almost an additional period compared, with the value in table 2. 
This experience of searching for one mode and finding another was not 
o The third mode (3) is seen from the final column of table 2 to a 
dynamically identical with the first, and fig. 5, to which the poe Oo 
the nodes have been added on the right, will serve to explain the ies 
The fourth mode (4) was difficult to obtain, the setting of the ie aa : : 
being delicate. In this test only, the nozzle was clearly pepagnEt kd 
more than half the period, and the curve in fig. 7 has been crawn wi 
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the aid of (3.12) and (3.13). At t,=1-25 the region between the two 
middle nodes was at its maximum elevation, thus the motion was assisted 
by the concentrated low-velocity fall. At the periphery the wave crest 
occurred at #,=2-25, and again the unobstructed jet was helpful. [t 
was not found possible to animate this mode with a jet travelling in the 
ait for only about half a period. 


Fig. 7 


05 


0 0.5 | ti 1.5 2 25 


Fall of jet from inclined nozzle at centre: circular symmetry, four 
nodes. Above: fall of jet from vertical nozzle at centre. 


The corresponding experiments with an almost vertical jet placed 
at the centre are summarized in the lower part of table 2; dynamically 
they are all the same as the fourth mode above, which also was principally 
excited by impact at a point in phase with the wave at the nozzle. For 
the slowest mode (1), one jet was used, but for the three higher modes 
the feebler jets were often momentarily upset by interference, and the 
usual timing over 50 cycles was difficult. Therefore, for these three 
experiments only, three jets were employed to facilitate timing. The 
A—t, curve now degenerates into the straight line with the numbered 
points plotted at the top of fig. 7. The wave at the centre was at a crest 
at t,=1:25, and the concentrated low-velocity fall occurred also at 
that instant. The unobstructed jet fell on both sides of the next crest 
at t,=2-25. We recognize that water falling on a rising wave that has 
passed its mean position does not entirely oppose the motion ; although 
the impulsive pressure is in the contrary direction, yet weight added to 
an elevation promotes the oscillation. Thus assistance was provided by 
the unobstructed jet. Experiments with an almost vertical jet placed 
close to the periphery were not pursued, because interference caused 
much water to be thrown out of the tank. 

The experimental results with the asymmetrical modes are set out in 
table 3. The test with the slowest mode (i) was spectacular, and the 
summit of the trajectory varied gracefully between heights about 102 
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and 52in. Since the jet was nearly vertical, there was much interference - 
but when the water from the submerged nozzle avoided collision, it fell 
close to the centre of the basin and had no influence on the oscillation. 
The second mode (ii) was produced in four ways with the nozzle placed 
in the mound position. The first three were with the nozzle trained on an 
antiphase point, and they are dynamically identical with the first and 
third entries in table 2, which also were carried out with impingement 
on an antiphase point. Now the modes and the fall of the obstructed 
jet were widely different in these five experiments, thus it appears that 
the main stimulus to the oscillations came from the properly timed 
impact of the unobstructed jet, the feebler obstructed jet having an effect 
only of minor importance. The fourth and the last lines in table 3 refer 
to experiments with impact on a point in phase, and the times of flight 
differed by half a period from the first three tests. The fifth test was 
with anti-phase conditions, the time of flight being a whole period 
greater than before. 


Table 3. Experiments with Asymmetrical Modes 


Mode * 


sec ED edits 10 Vir 


A, BR, BR V 2 (7% 


1/2 
; ) Nozzle at} Trained on 


Inclined jet 
e200" We lor 1-36 18 | 0-71 1:22 _ |periphery| other end of 

same diameter 

u 0-600 DG I IA (ey) aril 1-19 mound | hollow 

ii 0-595 25 |12 54 0-65 1-21 mound | nearer end of 

same diameter 


ii 0-597 8 94 8 | 0-58 0-68 mound | further end of 
same diameter 
ii 0-469 52 alee oe) 4- 2-22 mound | further end of 


same diameter 


Vertical jet 
ii 0-600 


OU 
bo 


H,=21 | 0-64 1-74 mound == 


When a pulsating jet is required in a fountain basin, the following 
procedure should. be adopted. First, the mode and the method of 
obtaining it must be decided, and then an estimate of the period is made, 
as explained in §2. The unobstructed jet is set at the value of i A 
calculated from the dimensionless time of flight given in table 2 or 3, 
and the basin level is slowly raised until pulsations set in. 
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§ 1. InTRODUCTION 


A cryostat with a hermetically sealed vacuum calorimeter has recently 
been developed in this laboratory for investigating the specific heats of 
solids at low temperatures. In this apparatus the use of an exchange 
gas and the operation of high vacuum equipment necessary with calori- 
meters of conventional design are avoided by the introduction of a 
mechanical arrangement which can make or break a thermal contact of 
the specimen with the bath as necessary.+ With slight modifications in 
technique the cryostat, which is now being used for measurements in the 
liquid helium range, could also be used for calorimetry up to 50°K. Its 
performance has been tested by some preliminary measurements on the 
specific heat of tin in the range 1-5—4-2°K. 


§ 2. DescripTION OF APPARATUS 


Figure 1 shows the contact making mechanism and the calorimeter 
chamber C. Inside C the specimen 8, is suspended by fine eureka wires 
of negligible heat conductance from the metal disc D, maintained at the 
helium bath temperature by the copper wire heat-leak L. It can be 
moved up and down without breaking the vacuum by means of the 
assembly at the top consisting of the metal bellows B, the spring S., 
and the nut N. The disc D, shields the specimen from any possible 
heating effects due to radiation or * hot molecules’ of gas from the top 
of the cryostat. The disc D, which is also maintained at liquid helium 
temperature, presses against the specimen when it is in the lowest position 
to give the thermal contact with the bath. 

The temperature of the specimen is measured with the carbon resistance 
thermometer T, (a } watt resistor of 12 ohm room temperature resistance, 
manufactured by Allen Bradley). This is cemented inside a hole drilled 
in the specimen. Its resistance temperature curve is determined by 
measurements with the specimen in contact with the calorimeter bottom, 
the bath temperature being found from the helium vapour pressure 


* Communicated by D. Shoenberg. 

+ While this work was in progress, we heard from Dr. D. Shoenberg, that a 
calorimeter based on the same principle had been developed at the Clarendon 
Laboratory, Oxford (Webb, F. J., Thesis, Oxford 1953). After submitting for 
publication we also heard of the work of Westrum, et a/. (1953, J. Chem. Phys., 
21, 419), in which a mechanical contact is used for cooling specimens in a 
different type of calorimeter. 
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tables published by Van Dijk and Shoenberg (1949). It exhibits a 
variation from 120 to 1300 ohms in the range 4-2-1-5°x. The resistance 
measuring bridge could detect changes equivalent to temperature changes 
of 10-4°K at 4°K with a measuring current of 10 microamps. This 
sensitivity increases at lower temperatures and is about seven times 
higher at 1-5°xK. 

The electrical leads to the thermometer and heating coil H, which are 
in thermal contact with the bath inside the metal capillary T,, are sealed 


Fig. 1 Fig. 2 
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off with wax at W; in an earlier design the leads were brought out 
through glass—metal seals in the calorimeter, but these often gave rise 
to leaks and were abandoned. At all other joints the vacuum system 1s 
sealed hermetically by soldering. On filling the cryostat with liquid 
helium the air inside the calorimeter is frozen out and a vacuum better 
than any obtainable with diffusion pumps is produced. The Breese 
is kept raised during this process to prevent it being frozen to the calori- 
meter bottom by frozen air. Later it is cooled. to the bath temperature 
by making contact with the bottom of the calorimeter. 


340 On A New Vacuum Calorimeter for Low Temperatures 


The specific heats at various temperatures are measured by isolating 
the specimen from the bath, passing a known electrical power for a 
definite time through the heating coil and measuring the rise in tempera- 
ture 47. In these measurements the specimen is always kept within 
0-04°x of the bath temperature so that the heat lost during the measuring 
interval along the leads and suspension wires is negligible (less than 
0-25°, of the total heat). When this temperature difference is much 
larger a correction has to be applied for this loss. 

The quantity of heat dissipated in the specimen is calculated from the 
known resistance of the coil, the heating time and the current. The 
change 47' in temperature, which is usually of the order of 0-02°K, could 
be measured to 1° accuracy. The accuracy of determination of specific 
heats is estimated to be better than 3°, down to 1-5°K. 


§ 3. RESULTS 

Figure 2 (a) is a plot of the variation of the specific heat of tin in the 
liquid helium range. It is seen that the scatter of points in the curve is 
less than 3°, except near the superconducting transition point where 
very small intervals (0-003°K) of temperature had to be used. The 
magnitude of the jump in the atomic heat at the superconducting 
transition point 3-71°K, is 0-0025 cal/°K. This agrees well with the 
value 0-0026 cal/°K calculated from Rutgers’s formula using the known 
value of dH,/dT at the transition point (H, is the critical magnetic field). 
Figure 2(b) shows the variation of C/T? with 7 observed by us. The 
results obtained by Keesom and Van Laer (1938) are also reproduced 
in fig. 2(c) for comparison. It can be seen that the two curves are similar 
in character, and there is good agreement between the absolute magnitudes 
of C/T. These results show that the new calorimeter compares well in 
performance with those of more conventional design. 
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XL. The Emission of Alpha Particles from Argon under Bombardment 
with Neutrons 


By E. H. Bretuamy and F. C. Fuack* 
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ABSTRACT 
Groups of particles have been observed from the bombardment of 
argon with neutrons using a high pressure gridded ionization chamber. 
The chief reaction is #°A(n, «)?’S with a measured Q-value of —2-5-+-0-1 
Mev, and transitions were observed to excited levels in 37S at 1-3-40-05, 
2-2+0-1, and 2-7-+-0-1mev. The cross section for the ground state transi- 
tion is approximately 30 «barns. 


§ 1. EXPERIMENTAL 


A GRIDDED ionization chamber with a counting volume of 1-5 litres con- 
taining argon at a pressure of 30 atmospheres was bombarded with 
neutrons from the *H(d,n)*He reaction. The resultant pulses after 
amplification were analysed on an 80 channel kicksorter (Hutchinson 
and Scarrott 1951). The result of such an experiment carried out at 90° 
to the 250 kev deuteron beam is shown in the figure. In addition to 
well defined peaks (A and B), evidence exists of other peaks as indicated 
by the arrows. 

Owing to slight alterations in overall gain and kicksorter bias, from day 
to day, which altered the position of the peaks by about 2%, the data 
from experiments which were made at intervals of several days cannot 
be immediately added on one graph, but particle groups corresponding 
to positions C and D shown in the figure are present in each of several pulse 
distributions. _The presence of a group at position E is also strongly 
indicated, although the data are not conclusive. 

The higher energy end of the distribution can be compounded from peaks 
of full width at half-maximum of 6% at energies of 11-6, 10-3, 9-4, 8-9 
and 8-1 Mev with relative intensities of 1-0, 1:45, 2-2, 1-5 and 1:0 respect- 
ively, together with a steeply rising lower energy component which has 
zero intensity at 9-4 Mev, but an intensity at 8-1 Mev which is more than 
twice that of the peak at this position. 

Polonium «-particles from a weak source were used to calibrate the 
energy scale and provide the resolution expected for a running time of 
several hours. The resolution obtained with a strong source when 


* Communicated by Professor P. I. Dee. 
+ Now at University College of the South West, Exeter. 
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amplifier and kicksorter drifts were negligible was 3%. An experiment 
with neutrons having a mean energy of 14-8 Mev gave the same Q values 
for the peaks A and B, but owing to the much larger spread in neutron 
energy (--0-3 Mev), the resolution was much poorer than in the figure, 
and the peaks C and D were not resolved. 

The neutron spectrum was examined by observing the recoil protons 
in a plastic phosphor. There was a small contamination of 2-5 Mev 
neutrons from the d(d, n)t reaction. 
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§ 2. INTERPRETATION OF RESULTS 


There can be little doubt that peak A corresponds to the reaction 
“A(n, «)®*S in which the sulphur nucleus is left in the ground state. 
The reactions (n, «) and (n, p) in 9%A (0-33°% abundant) are expected to 
produce peaks at H=16-1 and 14:2 Mev, of which there is little evidence. 
A is present only to 0-07°% and corresponding peaks are expected at 
13-9 Mev for the «-reaction which does not appear, and at 10 Mey where 
it would be hard to detect. (n,d) reactions would give peaks at 7-9 and 
8-1 Mev respectively for these two nuclei, but owing to the low abundances, 


excessive cross sections would be required for these reactions to give a 
detectable contribution to the peaks found. 


Pretty gp Cage 
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This leaves, then, only reactions in “A to account for the remaining 
peaks. The Q values for the reactions 38A(n, p)36Cl and 3A(n, p)3Cl 
would indicate that the Q value for the reaction 40A(n, p)*°Cl is certainly 
less than —4 Mey and is likely to be about —8 mev, so that all the assign - 
ments can be made to the reaction #°A(n, «)37Cl. The resulting Q-values 
and excited levels in °’Cl are shown in the table. The cross section for 
the ground state transition was calculated from the integrated deuteron 


current incident on the target and is unlikely to be correct to better than 
50%. 


Peak Intensity E peak (Mev) Q(mev) | Levels in 3’Cl (mev) 
A 1-0 11-6 —2-5+0-1 ground state 
B 1-45 10-3 —3°8 1-3-+.0-05 
C 2 9-4 —4:7 2-2+0-1 
D 1-15 8-9 —5:-2 2-7+0-1 
(#) 1 8-1 —6:0 3-5+0-2 


The energy limits for charged particles from the (n, d) and (n, pn) 
reactions in4°A are marked in the figure. The Q-value for the *°A(y, p)89Cl 
reaction as determined by Wilkinson and Carver (1951) has been used to 
determine these limits. The steps in the curve may be associated with 
these reactions. 

The Q-value for the ground state reaction together with the mass- 
spectrometer determinations of the *°A and 3’Cl masses can be used to 
compute the end-point energy of the B-ray (10%) to the ground state of 
87Cl from the 5 min 2’S decay. A value of 4-9-++0-2 Mev is obtained which 
is not in strong disagreement with the value of 4:3-10-3 Mev obtained by 
an Aluminium absorption determination (Bleuler and Zunti 1946). This 
Q-value also confirms the assignment (Endt and Kluyver 1954) of the 
a-particles observed by Graves and Coon (1946) with 2-5 Mev neutrons 
to the 3*A(n, «)?3S reaction. 
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XLII. Spallation of Elements in the Mass Range 51-75 


By 8. G. Rupstam 


The Gustaf Werner Institute for Nuclear Chemistry, University 
of Uppsala, Uppsala, Sweden* 


[Received November 17, 1954] 


SUMMARY 


Empirical formulae for the yields of the products obtained in spallation 
reactions with high energy particles are constructed. These formulae 
reproduce most of the experimentally determined cross sections within 
a factor of 2, and they simplify the comparison between different spallation 
investigations. 


§ 1. [yTRODUCTION 


Tue spallation of medium weight elements has been fairly well studied. 
The mass region best investigated is probably 51-75. In this region 
the following spallation studies are found in the literature: vanadium, 
manganese, and cobalt + 187 Mev protons (Rudstam 1953), iron + 340 
Mev protons (Rudstam, Stevenson, and Folger 1952), cobalt + 240 Mev 
protons (Wagner and Wiig 1952), cobalt + 370 Mev protons (Belmont 
and Miller 1954), copper -+ 60 Mev protons (Carleson 1954), copper +190 
Mev deuterons, 340 Mev protons and 380 Mev alpha particles (Batzel, 
Miller, and Seaborg 1951, Miller, Thompson and Cunningham 1948), 
copper + 400 Mev neutrons (Marquez 1952), copper + 2-2 Bev protons 
(Friedlander, Miller, Wolfgang, Hudis, and Baker 1954), zinc + 340 Mev 
protons (Worthington 1952), and arsenic + 190 Mev deuterons (Hopkins 
and Cunningham 1948). It is seen that the conditions of the investigations 
differ widely. It is of great interest to compare these different groups of 
spallation reactions,, but this is not easily done if the spallation yields 
are given merely as a set of tabulated cross sections. 

However, for vanadium bombarded with 187 Mev protons, the cross 
sections can be expressed by a formula of the following type (Rudstam 


i= 4 
hee o(Ay Z,)=exp(P)Aj—-Q--R(Zie SA) ae 


where o(A,, Z,) is the formation cross section of the nuclide with the mass 
number A, and the atomic number Z,, and P, Q, R, and S are constants. 

If the spallation groups mentioned before could be characterized by an 
equation of this type, this would greatly facilitate the comparison between 
them. The purpose of the present study was to investigate the application 


* Communicated by Professor The Svedberg, F.R.S. 

+ By a group of spallation reactions, or a spallation group, is meant all the 
spallation reactions associated with the bombardment of a given element (X) 
with particles of a given kind (a) and energy (Z). Notation: G(X; a; Z). 
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of the formula (1) to the spallation of elements in the mass range 51-75 
and then to use the results for a comparison between these spallation 
groups. 

§ 2. DETERMINATION OF THE Cross SECTION FoRMULA 


The cross section formula (1) for a group of spallation reactions can be 
determined from the measured cross sections in the following way. The 
formula is written : 

log, o(A,, Z;)=P.A,—Q—R(Z,—S.A,). Page) 

An arbitrary value is chosen for the parameter S, and the method of 
least squares is then used to determine the best values of the parameters 
P,Q, and & from the known data of log, o(A,, Z,), A;, and (Z,—S.A,)2. 
The root-mean-square error « of the logarithm of the cross sections is 
calculated, and different values of S are tried until the minimum of « is 
found. The corresponding values of P, Q, R, and S give the best formula 
of the type (1) for the cross sections, and the value of «,,,, shows how 
successful the formula is, because the average ratio betweenexperimental] 
and calculated cross sections is given by the factor e**. 

The experimental cross sections used for the determination of the 
cross section formulae should be the independent cross sections for 
primary formed spallation products. However, the independent cross 
section of a nucleus is generally much lower than the independent cross 
section of its daughter nucleus for the actual spallation reactions. For 
this reason the measured cross sections can be regarded as independent. 

Since the cross section formula cannot be expected to hold over all 
mass numbers, the mass region used in the calculations was arbitrarily 
limited to between 30 and 2 mass numbers below the lowest mass number 
of the target isotopes. 


§ 3. RESULTS OF THE CALCULATIONS 

The parameters P, Q, R, and S, and the factor e* are tabulated in table 1 
for twelve different groups of spallation reactions. 

The cross section formula fits remarkably well as is seen from the 
values of e€. In most cases this factor is less than 2, which means that 
the cross section formula reproduces the experimental cross sections 
within a factor of 2 on the average. This is satisfactory since the cross 
sections often vary by a factor of 10? or more, and moreover the experi- 
mental cross sections may sometimes be inaccurate to a factor of 2. 
Thus one may conclude that the cross section formula (1) is a good 
approximation to the cross section distribution for spallation of elements 
in the actual mass region. For a given spallation group, the cross section 
distribution can then be roughly characterized as follows : . 

(1) The yields* of the spallation products vary regularly with their 
atomic number and mass number. ‘This is indicated by the mere 
possibility of expressing all the cross sections by a simple formula. 


* The yields are proportional to the cross sections. 
ZAZ 
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(2) The yield-mass curve, that is the total isobaric yield o(A) versus 
the mass number, increases exponentially with the mass number. This 
is seen by the following argument. 


o(A)=Zo(A, Z)=exp [P. A—Q]2 exp [—R(Z—S . A)?] 
Z Z 


~exp[P.A—Q] | é exp [—R(Z—S . A)? dZ. 


The limits of integration, Z, and Z,, are the smallest and the largest 
obtainable value of Z. If R is not too small this integral will not differ 
much from the integral with the limits — co and +00. One then obtains 


o(A)~exp[P.A— ay exp[—R(Z—S.A)*] dZ 
7 —(n/R)2-exp[P.A—Q]. . (3) 
Table 1 


Spallation group 


G(V ; p; 187 mev) 


G(Mn ; p; 187 Mev) 


G(Cu ; d; 190 mev) 


G(As ; d; 190 mev) 


G(Fe ; p ; 340 Mev) 


G(Co ; p; 370 Mev) 


G(Zn ; p ; 340 Mev) 


G(Cu ; p ; 340 Mev) 


G(Cu ; p ; 60 Mev) : 
+0-119 


G(Cu ; p ; 2200 Mev) 0-072 
+0-015 


G(Cu ; n; 400 Mev) 0-350 
+-0-067 


G(Cu ; « ; 380 Mev) 0-236 
+ 0-036 
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The parameter P is thus a measure of the steepness of this exponential 
increase in the yield—mass curve. 

(3) The charge distribution curve, that is, the independent yield versus 
- the atomic number for isobars, is a gaussian function of the atomic 
number with a precision index (=R1/2) which is independent of the mass 
number. In this connection it is of interest to note that the assumption 
of a charge distribution curve independent of mass number is consistent 
with evidence from the study of yields of fragments in the fission of 
uranium with thermal neutrons (Pappas 1953). 

(4) The most probable primary charge, that is, the ‘ atomic number ’, 
generally non-integral, giving maximum yield in the charge distribution 
curve, is a linear function of the mass number, passing through the 
origin and having a slope equal to S in a Z versus A plot. 

It may be remarked that the constant @ merely fixes the scale of the 
cross sections and is of minor interest here, since the measured cross sections 
are only relative in most cases. 

The measured and calculated cross sections for the spallation groups 
discussed in this paper can be compared in the tables 2-13. 


Table 2. G(V; p; 187 Mev). Cross Sections in millibarn. Experimental 
data: Rudstam (1953) 


: oO oO . oO oO 
Nuclide measured | calculated Nuelide measured | calculated 

22 Nai 0-030 0-020 258) 3:6 4-7 
"Mo 0-028 0-032 47Ca O11 0:08 
3181 0-18 0-16 43S¢e 4-0 2-5 
SA B- 0-34 0:27 44Sc 17 | 12 
33? 0:32 0-34 46Sc 20 DAS 
35S 0-61 0-74 47S¢ 10 12 
BIC} 0-031 0-056 *8Sc 4:8 2-4 
38C] 0-74. 0:97 457i 6-1 3°8 
220] 0:30 eld a sok A 
42 : ise ) 
Fe Ae 1:8 48Cr 0:23 0:59 


Table 3. G(Mn; p; 187 mev). Relative Cross Sections. Experimental 
data: Rudstam (1953) 


| “4 
a o o 
Nuclide Ae eee Pe suinted Nuclide measured | calculated 
d 5 N a4 * 0-31 
3181 0-027 0:022 45Ca, 0-25 0: ‘ 
ss 0:034 0-047 47Ca, 0-007 | 0-007 | 
ral 5@ 0:39 0:40 ad Wi 1:00 | 0:94 | 
43K 0-21 0-14 49Cr 9-9 | 2-3 | 
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Table 4. G(Cu; d; 190 Mev). Relative Cross Sections. Experimental 
data: Batzel, Miller and Seaborg (1951) 


+ oO oO . oO oO 
Nuclide measured | calculated Nuclide measured | calculated 
32P 0:0005 0-00018 56VIn 0:10 0-096 
38] 0-00032 0-00072 o2Be 0-003 0:0024. 
39] 0:00007 0:00039 53He 0-047 0-022 
4517 0-002 0:0029 59Fe 0-045 0-083 
48V 0-035 0-016 55Co 0-022 0-025 
43Cr 0-01 0-0071 81Co 0-094 0-192 
51Cr 0:03 0-061 57Ni 0-054 0-037 
5IMn 0-029 0-011 60Cu 0-20 0:30 
52MIn 0-087 0-043 81Cu 1-00 0:98 


Table 5. G(As; d; 190 Mev). Relative Cross Sections. Experimental 
data : Hopkins and Cunningham (1948) 


. oO oO . oO Oo 
Nuclide measured | calculated Nuclide measured | calculated 
52MIn 0:0002 0:00065 61Cu 0-1 0-018 
56Mn 0-002 0:0025 64Cu 0-1 0-057 
59He 0-005 0:0040 67Cu 0-02 0-013 
55Co 0-003 0-0007 697n 0-12 0-032 
81Co 0-003 0-0093 6Ga 0-1 0-13 
57Ni 0:0002 0-0011 8Ga 0-2 0-26 
S5Ni 0-001 0:0054- 71As 0:3 0-85 

66Ni 0-002 0-0018 72As 1-00 1-1 
Cu 0-06 0:0068 


Table 6. G(Fe; p; 340 Mev). Cross Sections in millibarn. Experimental 
data: Rudstam, Stevenson, and Folger (1952) 


: Oo o : oO oO 

Nuclide measured | calculated Nuclide measured | calculated 
318i 0-12 0-07 47S¢ 0-84. 0:98 
2? 0-044 0-15 48S 0:45 0:20 
35S 0:23 0:23 4577} 3:7 1:9 
84] 0-11 0-07 ATV 2-4. 2°7 
38C] 0-17 0-21 48V 10:3 7 aft 
89C] 0-045 0-059 48Cr 0:80 0-73 
SAR 0:7 0:75 49Cr 4:8 4:0 
43K 0:4 0:25 51Mn 4-0 5:7 
45Ca, 0:56 0-50 2Mn 13 18 
47Ca 0-007 0-010 52He 0:68 1-32 
46Sc Bie 2-4 
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Table 7. G(Co; p; 370 Mey). Cross Sections in millibarn. Experimental 
data: Belmont and Miller (1954) 


. oO oO . oO Co 
Nuclide measured | calculated Nuclide measured | calculated 

32/—P 0:30 0-21 49V7 31-1 iil 
ob) 0:03 0:20 49Cr 4-1] 2:5 
390] 0-50 0:14 PLCr 27:5 16 
38k 0:31 0:13 5IMn 4-1] 3-2 
Ae 0-85 122 52MIn 19-7 12 

a3 KK 0:50 0:52 54In 25-4 31 
o(Ca; 0-66 1:00 56Min 3°8 8:3 
47Ca 0:06 0-04 52He 0:28 0:61 
46S 2-2 3:7 53He 1-7 4-3 
4577} 3°5 1-4 55e 37:0 39 
47 2-1 1:9 55Co 5:2 5:6 
a8 0-6 5:8 56Co 


Table 8. G(Zn; p; 340 Mev). Cross Sections in millibarn. Experimental 
data : Worthington (1952) 


a o oO : o (or 
Nuclide measured | calculated Nuclide measured | calculated 
32P 0-02 0-063 Mn 2 1-1 
ee oe [ae ees 
42 ANS * n ‘ 2 
& ee 0-12 52He 0-22 0-17 
45Ca, 0-22 0:24 Fe 0:67 0-61 
47Ca, 0:0055 0:0054 55Co 1:5 1-9 
Bese | oh | Sea ae 
. . O 

a Ss 2-0 56Ni 0-090 0:25 
48Cr 0-20 0-12 St 1-7 ae 
49Cy 0-8 61Cu 33 5 


Table 9. G(Cu; p; 340 Mev). Relative Cross Sections. Experimental 


data: Batzel, Miller, and Seaborg (1951) 
{2 ee eT 


ee re oO oO 
. Nuclide Be ited Perini Nuclide measured | calculated 
31 0:17 
2A 0:0053 0:0019 52MIn 0:3 
34] 0:0003 0:0006 54Min 0:53 te 
38C] 0:0025 0:0038 56Mn 0-11 0: 
39] 0-0004 0:0013 52He 0-008 0)-0065 
45Ca 0:007 0-011 53He 0:073 0-056 
47Ca 0:0004 0:00033 5b We 0-49 0-63 
46Sc 0:075 0-047 59He on ogee 
as 0:033 0:023 Co -] oul 
se 0-044 0-015 **Co oA ve 
48V 0-032 0-078 200 2-5 | ee 
49Cr 0-041 0-030 $1Co O21 | he 
51Cr 0:31 0:25 57Ni 0-079 BU 
51Mn 0-071 0-041 610u 1:00 2-4 


eee 
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Table 10. G(Cu; p; 60 Mev). Relative Cross Sections. Experimental 
data: Carleson (1954) 


. oO oO : oO oO 
Nuclide measured | calculated Nuclide measured | calculated 
54)Vin 0-0012 0-0014 58C0 0-25 0:07 
56MIn 0:00091 0:00075 61Co 0-013 0-018 
59He 0:0018 0-0019 57Ni 0-014 0-006 
56Co 0-0022 0-0086 60Cu 0:32 0:33 
57Co 0-031 0-031 61Cu 1-00 1-51 


Table 11. G(Cu: p; 2200 Mev). Cross Sections in millibarn. Experimental 
data: Friedlander et al. (1954) 


A oO oO . Co oO 
Nuclide measured | calculated Nuelide measured | calculated 

32p 6-4 9-5 45Cr 0-5 0-3 
33p 0-9 23 49Cr 2-2 Lei 
358 0:95 2-8 510r 19 9-5 
340] 0-89 0-58 51Vin 1-4 1-6 
38] 0-67 2-2 Mn 7:0 5-2 
390] 0-40 0-79 54\In 12 12 
421 3°7 3°5 56Vin 2:3 3°3 
43K 1:3 1:5 52He 0-19 0-24 
45Ca 1:1 2-1 53Re 1:3 1:5 
47Ca 0-10 0:09 55 We 12 11 
46Se 6:5 5-6 59Re 0:6 1-2 
47S¢ 3:5 2:8 55Co 1:5 1-4 
48Sc 2-4 0:8 81Co 3°6 1:7 
5T7 2°5 1-8. 57Ni 0-6 1:2 
47VV 2:7 1:7 60Cu 1:8 5:2 
48V 9-6 5-0 61Cu 6:5 14 
49VV 97 8-5 


Table 12. G(Cu; n; 400 Mev). Cross Sections in millibarn. Experimental 
data: Marquez (1952) 


. Oo Oo . Oo oO 
Nuclide measured | calculated Nuclide measured | calculated 

Tj 0-0078 0-071 He 2-44. 2:3 
490r 0:35 0:20 2 Oo) 0-415 0:87 
5IMn 0:76 0:39 58Co 62 24 

52MIn 4-68 1:5 81Co 3:78 5:8 
56Min 2°81 2:8 57Ni 0-54 1-4 
52he 0:134. 0-067 81Cu 15-4 47 


53He 1-24 0-54 
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Table 13. G(Cu; «; 380 Mev). Relative Cross Sections. Experimental 
data: Batzel, Miller, and Seaborg (1951) 


. Co oO . (on Oo 
Nuclide measured | calculated Nuclide measured | calculated 
38] 0:0022 0:0014 56Mn 0:05 0-050 
49Cr 0-006 0-012 59He 0-034. 0-040 
51Cr 0-016 0:049 5INi 0:066 0-049 
51MIn 0:025 0-018 60Cu 0-12 0-21 
52Min 0-08 0-044 61Cu 1:00 0:44 


§ 4. COMPARISON BETWEEN SPALLATION GROUPS 


The cross section distributions of different spallation groups can now 
be compared by studying the dependence of the parameters P, R, and S 
on the composition of the target and on the kind and energy of the bom- 
barding particles. It must be remarked in advance, however, that it is 
not possible to obtain more than general trends from the available data. 

The figs. 1-8 indicate how the parameters of the cross section equation 
vary with the energy of the bombarding particles and with the mass 
number Ay or the atomic number Zp). A, is defined as the mass number 
of the nucleus formed by adding the bombarding particle to the target 


Fig. 1 
&£=190 Mey 
Pp 
0.4 
& 
i d 
mle ¢ 
0.2 
0.1 
0 
50 60 70 80 4 


nucleus, in other words the ‘ compound ’ nucleus. Zy is the corresponding 
atomic number. Aj, defined in this way, applies only to monoisotopic 
targets. For other target elements an average of A, has to be used. 
The kind of bombarding particle is denoted by the letters p, d, n, or « 
appearing above the corresponding point in the figures. aap onl 
The dependence of the parameter P on A, can be studied in figs. 


1 and 2 for the energies 190 mev (fig. 1) and 340-400 Mey (fig. 2). It is 
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seen from the figures that P is almost independent of A, in the mass region 
investigated. It should also be noted that P does not depend on the kind 
of bombarding particle, with the possible exception of neutrons. Figure 3 
shows P as a function of the energy. The parameter P decreases with 
increasing energy. This is to be expected. The steepness of the yield— 
mass curve should decrease with increasing energy of the bombarding 
particles. 
Fig. 2 


© = 340-400 Mev 


; na 
0.4 | 
0.3 
iesi aes 7 
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The parameter R is plotted against A, for proton bombardments in 
fig. 4* and for bombardments with other projectiles in fig. 5. The energies 
of the bombardments (in Mev) are denoted in the figures. 

Figure 4 shows that R is quite insensitive to variation of A, and of the 
energy. Thus the charge distribution curve has almost the same width 


Se eae eee 
* Three points at Aj>—64-5 (spallation of copper) have been spread out in the 
figure. 
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in the spallation of copper with 60 mev and with 2200 mev protons. 
One might expect R to depend on the isotopic composition of the target, 
but this does not appear from the figure. 

Figure 5 indicates that the parameter R is somewhat smaller for 
deuteron and alpha particle bombardments than for proton bombardment. 


Fig. 4 
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It is not very profitable to compare the values of the Caen S 
for different spallation groups. Instead, one may study the agate 
of the factor U=S—(Z,/A,), which is a measure of = preference or 
neutron emission. This factor will be zero if the most probable see 
between emitted protons and neutrons (where both free and bounce 
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nucleons are considered) is the same as the ratio between protons and 
neutrons in the compound nucleus. If U is positive, neutron emission 
is more probable than proton emission. The factor U is plotted versus 
Z, for 190 mev in fig. 6 and for 340-400 mev in fig. 7 and versus the energy 
of the bombarding particles for Z)=30 in fig. 8. 
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Figure 6 indicates an increase of U with increasing Z). This increase 
is to be expected since the potential barrier for charged particles rises 
with increasing Zp. ; 

Figure 8 shows that U decreases with increasing energy. This trend 
is not surprising. The potential barrier, which favours the emission 
of neutrons, is relatively more important when the energy of the bombard- 
ing particles is low. 


Fig. 8 
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By integrating the cross section equation it is possible to estimate the 
total cross section, o(A,; A,), for formation of products with mass 
numbers between A, and A,. Thus one obtains : 


A, co 
o(Ay; Ay~| dA| exp [PA—Q—R(Z—S. Ay] dZ 


-Q 1/2 
= 5—gur [exp (PA»)—exp (PA,)]. eens) 


§ 5. Usk oF THE Cross SECTION FoRMULA 


It is of practical interest to be able to estimate cross sections of spalla- 
tion products as a function of the mass number (and the atomic number) 
of the target and of the kind and energy of the bombarding particles. 
A systematic study of spallation reactions, preferably with ee 
targets, and determination of cross section distributions might make ae 4) 
an estimation possible since the parameters of the cross section foray a 
(1) could then be evaluated as functions of the variables mentioned above. 
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XLUI. CORRESPONDENCE 


Radioactive 69As and As 


By F. D. 8S. Buremenr and E. G. Provr 
Atomic Energy Research Establishment, Harwell 


[Received January 13, 1955] 


THE lightest radioarsenic isotope previously known is As. This letter 
describes a new neutron deficient isotope ®As. 

Samples of 20 mg of ‘Specpure’ germanium dioxide were irradiated 
with protons at a series of energies from 8-60 Mev. Radioarsenic was 
then isolated by dissolving the germanium dioxide in 3ml. of 27N 
hydrofluoric acid, adding 5 ml. of 6N ammonium hydroxide and preci- 
pitating the arsenic as sulphide by passing in hydrogen sulphide. After 
a second precipitation with inactive germanium hold-back carrier, the 
arsenious sulphide was isolated for counting. 

At proton energies of 8-20 Mev the radioarsenic showed half-lives 
of 52 minutes (As), 26 hours (7As and 76As) and 17-6 days (74As). 
Between 20 and 35 Mev a new 15 minute activity appeared in appreciable 
yield (®°As), due to a (p, 2n) reaction on “Ge. Between 40 and 60 Mev 
the initial half-life was ~10 minutes, suggesting the formation of a 
mixture of ®°As with a shorter lived activity, of half-life ~7 minutes, 
possibly ®*As produced by a (p, 3n) reaction on Ge. 

The mass assignment of the 15 minute activity to ®°As was made by 
milking off and identifying the radioactive germanium daughter. 
Radioarsenic isolated from germanium dioxide irradiated with 30 Mev 
protons was prepared as a solution in hydrochloric acid containing inactive 
germanium carrier, and the germanium was then quantitatively isolated 
by distillation for two minutes in a current of chlorine. Successive 
separations of germanium were carried out at regular time intervals, 
fresh germanium carrier being added on each occasion. Germanium 
sulphide sources for counting were prepared from each distillate. All 
the germanium sources showed a half-life of 40 hours, due to “Ge, and 
the yields for successive separations decreased at a rate indicating a 
15 minute half-life for the parent, °°As. 

The ®As emitted positrons, the maximum energy of which, as deter- 
mined by absorption in aluminium, was 2-9 Mev. The abeorption curve 
did not show any evidence of lower energy components. The gamma 
radiation was examined with a Nal(Tl) scintillation spectrometer and. 
a single channel pulse analyser. In order to avoid the interference 
caused by the backscatter peak of the annihilation. radiation, the source 
was mounted on a very thin aluminium foil in a magnetic field sufficient 
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to prevent any positrons from reaching the scintillator crystal. A single 
gamma-ray with an energy of 0-23 Mev was found. In both the positron 
and gamma-ray energy determinations corrections were made for the 
effects of the longer lived activities. The gamma-ray spectrum of 52 
minute 7As was also examined and gamma-rays of energies 1-06, 1-74 
and 2-04 Mev were found in agreement with the results of Aten and 
Verkerk (1952). Additional gamma-rays of energies 0-18 and 1:36 Mev 
were also found. 
REFERENCE 
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The Elastic Scattering of Protons by Helium 3 


By D. R. SwEETMAN 
Department of Physics, The University of Birmingham* 


[Received January 29, 1955] 


Tue angular distribution for the elastic scattering of protons by helium 3 
is of interest for comparison with the theoretical results given by Swan 
(1953). 

By operation of the Birmingham 60 in. cyclotron at }? the magnetic 
field used for ?D+ acceleration and using helium 2°% enriched in *He 
fed into the ion source, an internal beam of *He ions suitable for trans- 
mutation experiments has been obtained (Fremlin 1952). Unfortunately, 
the unsuitable shimming at the low magnetic field results in the loss of 
most of the particles before they reach the extraction radius. Thus 
the present external 7He beam is too small for scattering purposes and 
an attempt has therefore been made to use the internal *He beam with 
a specially designed scattering chamber. 

As space was very restricted a nuclear plate was used as a detector, 
the method being similar to that used by Chadwick, May, Pickavance 
and Powell (1944) except that solid targets were used. The scatterer 
was mounted at 45° to the incident beam which was defined by a stop 
system ; the natural collimation of the beam was found to be sufficient. 
The beam was received in a Faraday cup which was connected to a 
condenser, the voltage across which was measured with a Farmer 
electrometer voltmeter. Scattered particles entered the nuclear plate 
at a glancing angle of between 5 and 10 degrees. As the position of the 
beam hitting the scatterer was known accurately, relative cross sections 
could be determined with an accuracy essentially limited by the number 
of tracks that were measured in the plate. The beam was located and the 
cyclotron adjusted to resonance using a sensitive current meter connected 


to a mechanically operated shutter which was raised for the exposures. 


* Communicated by J. H. Fremlin. 
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Trouble was experienced at first due to fogging of the plates by intense 
X-rays emitted by the dee feelers and dee edges. It was found, however, 
that thick layers of copper on the top and bottom of the chamber and 
a lead cap on the end facing the dee feelers were sufficient to reduce the 
fog background to negligible proportions. 

Since it was difficult to determine absolute cross sections directly 
(due to the large vertical spread of the beam at the scatterer position) 
cross sections were normalized using a plate exposed to particles 
elastically scattered from gold, the cross sections for which were calculated 
from the Rutherford relationship. 


EXTREME CURVES 
GIVEN BY SWAN 


Cross section (barns/steradian) 


RUTHERFORD CURVE 


0 30 60 90 120 150 180 
C—M angle (degrees) 


Angular distribution for p-*He elastic scattering. C—M energy 3-72 Mev. 


The apparatus has been used to scatter 14-9 Mev *He ions from protons 
(C-M energy=3-72 Mev) the latter being conveniently obtainable in 
the form of polythene foil. A proton group of distinct range was obtained 
at all angles to the incident beam up to 70° when the range of the seat- 
tered protons became so small that some were lost in the scatterer. 
In order to establish the nature of the background to be subtracted from 
the proton group a second plate was exposed using a carbon scatterer. 
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The resulting angular distribution when converted to the C-M system 
for the scattering of protons by He is shown in the figure. The probable 
errors shown have been calculated on a purely statistical basis from the 
number of tracks measured at each angle. The method of normalization 
employed limits the absolute accuracy of the cross sections to about 10%. 

A preliminary result has been obtained for the higher C-M energy of 
6-45 Mev when the minimum appears to have moved to about 104° 
C-M angle. 

The angular distribution given agrees with the general trend of 
the experimental results at lower energies given by Famularo, Brown, 
Holmgren and Stratton (1953) which also show a pronounced minimum 
near 90° and a steep rise towards 180°. 

The distribution shown is not in agreement with the theoretical curves 
given by Swan (1953) who predicts no pronounced rise towards 180° 
for this energy. However, he suggests* that the agreement may be 
considerably improved by not using approximate Gaussian ground-state 
wave functions. 
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XLIV. REVIEWS OF BOOKS 


Optics : Lectures on Theoretical Physics, Vol. IV. By ARNOLD SOMMERFELD. 
Translated from the German by Otto Laporte and Peter A. Moldauer. 
[Pp. vii+-383.] (Academic Press Inc., New York.) Price $6.80. 


Tuts book forms Volume IV of the author’s series with the general title 
Lectures on Theoretical Physics. Tt includes a thorough treatment of the theory 
of diffraction regarded as a boundary value problem in electro-magnetic 
theory, and this is the most valuable part of the book. The rationalized 
M.K.S. system of electrical units is used. Unfortunately, the sections on 
elementary diffraction theory are not so satisfactory. They are in places 
inaccurate and suffer from the complete omission of any vector diagrams, 
which are so valuable in understanding diffraction phenomena. On page 220 
it is stated that the intensity due to the first Fresnel zone for a circular aperture 
is equal to that of the unobstructed wave (instead of twice that of the un- 
obstructed wave). Figure 64 is also incorrect. On page 215 the increase in 
strength of radio signals at the point on the Earth antipodal to the transmitting 
antenna is attributed to diffraction round the spherical Earth, but in practice 
this is due to the trapping and focussing effect of the ionosphere. 

The treatment of anomalous dispersion is good, and includes an outline of 
the wave-mechanical theory. Other interesting topics, not usually found in 
standard textbooks on optics, include a quantum treatment of Doppler effect, 
the theory of Cerencov radiation, and a thorough discussion of wave propagation 
in dispersive media. The section on white light gives only a qualitative 
discussion and seems inadequate. 

‘The methods of Fourier transforms and contour integration are frequently 
used, but not explained, reference being made to Volume VI of the series for 
a treatment of these mathematical topics. 

Because of the omission of most of the experimental material, the book 
cannot be recommended as a standard textbook for students of experimental 
physics, though it should be useful as a reference book. Students of theoretical 
physics should find it valuable and stimulating. Brite: 


Radio-Isotope Conference 1954, Vol. I, Medical and Physiological Applications. 
Edited by J. E. Jounston, R. A. Farres and R. J. Mitterr. [Pp. xi+418.] 
(Butterworths Scientific Publications.) Price 65s. 


Tuts volume is a report of the Second Isotope Conference held in Oxford 
during July, 1954, the first one having been held in July 1951. There are 
46 papers under the following headings :—Therapy, Diagnosis, Animal Physiology 
and Pathology, Biochemistry, and Plant Nutrition and Allied Subjects. The 
whole Conference, of which the above proceedings report a part, consisted of 
70 papers chosen from 280 offered. The editors, in choosing the papers to 
be read, gave preference to new ideas and to new methods and stress the 
importance of new results. The present volume testifies admirably to their 
choice. Pee, . 
Whereas in 1951 the emphasis was still on the measuring and assaying of 
radioactive isotopes, it has now definitely moved over, in both clinical work 
and research work, to their uses. Nevertheless, techniques of measurement 
are still improving and from this Conference one can see, for eet Pom, 
during the past three years, the crystal scintillation counter has replace Lt he 
Geiger counter in this field. The uses of radio-iodine, ae ore, gee 
radio-phosphorus still form a large part of the Conterence but phe field s 
obviously expanding rapidly ; however, uses for such isotopes as *Y, °°Fe anc 
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206Bi, which had not appeared three years ago, are also described. Tracer 
work in metabolic studies is, of course, rapidly expanding; biochemical 
compounds labelled with C and 3H are now common tools. This again is 
reflected in the number of papers in this field. 

The proceedings of the Conference demonstrate, if any were needed, how 
radioactive isotopes have developed into being one of the established working 
tools for both the medical specialist and biochemist. This Conference gives 
an excellent view of the present situation and points the way for further 
advances. 

The volume is well bound, clearly printed and singularly free from mistakes, 
though at 65s. the price seems rather high. Both the editors and publishers are 
to be congratulated on the short time which elapsed between the Conference 
and the date of publication. Colas. 


BOOK NOTICES 


Luminescence with Particular Reference to Inorganic Phosphors. A symposium 
held by the Electronics Group of the Institute of Physics at the Cavendish 
Laboratory, Cambridge from the 7 to 10 April, 1954. British Journal 
of Applied Physics, Supplement No. 4. Price 25s. 


Petrology for Students. An introduction to the study of rocks under the 
microscope. By ALFRED HarKER. [Pp. viii+283.] (Cambridge 
University Press.) (8th edition, revised.) Price 18s. 


Mathematics of Engineering Systems (Linear and Non-Linear). By DErex F. 
Lawpen. [Pp. 380.] (Methuen & Co. Ltd.) Price 30s. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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D. SHAW and B. J. MASON | Phil. Mag. Ser. 7; Vol: 46, Pl. 3. 


(¢) (f) (7) 

Examples of crystals grown on the metal plate. (a) Hexagona | plates ( —15°e), 
(b) prisms (—6°c), (c) prisms (—22°c), (d) prisms (—33°C), (e), (f), (9), 
prisms with pyramidal faces (—22°c). 
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Oscillograph traces of 8, and C, counter signals. The scale is in units of LO~® sec. 


(a) An S-particle, delay 36 x 10~ sec. 

(6) A fast meson, delay 1 x 10~* sec. 

(c) An electron shower of high density ; many particles enter the counter Cg. 
The delay between the 8, and C, signals is measured between points 
at which the trace is deflected through a fixed distance of about half a 
trace width. 
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